Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,441,988 T→C 8.1% I469I (ATT→ATC thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,441,9880T→C 8.1% 40.5 0.0 ‑0.0233I469I (ATT→ATCthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (10/9):  ref base (107/107):  total (117/116)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.79e-01

Only 100 of 240 total aligned reads displayed.

GCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGG  >  NC_007779/3441888‑3442088
                                                                                                    |                                                                                                     
  cACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgc                                                                                                     <  1:6959059/101‑1 (MQ=255)
  cACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgc                                                                                                     <  1:5320485/101‑1 (MQ=255)
   aCGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgcg                                                                                                    >  1:75990/1‑101 (MQ=255)
    cGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgcgc                                                                                                   <  1:12533351/101‑1 (MQ=255)
    cGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgcgc                                                                                                   >  1:4860511/1‑101 (MQ=255)
    cGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgcgc                                                                                                   <  1:2116037/101‑1 (MQ=255)
    cGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTgcgc                                                                                                   >  1:14724886/1‑101 (MQ=255)
     gTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGcc                                                                                                  >  1:2765903/1‑101 (MQ=255)
      tGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCg                                                                                                 <  1:7744574/101‑1 (MQ=255)
       gCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCgg                                                                                                <  1:7274882/101‑1 (MQ=255)
          gATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCt                                                                                             >  1:14551610/1‑101 (MQ=255)
          gATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCt                                                                                             <  1:12043100/101‑1 (MQ=255)
           aTGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTa                                                                                            >  1:15061203/1‑101 (MQ=255)
                   gatCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACt                                                                                    <  1:10155205/101‑1 (MQ=255)
                    atCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACtt                                                                                   >  1:1922054/1‑101 (MQ=255)
                     tCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTc                                                                                  >  1:10184220/1‑101 (MQ=255)
                      ccgccgCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCa                                                                                 <  1:5454182/101‑1 (MQ=255)
                       cgccgcAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCAc                                                                                <  1:11681633/101‑1 (MQ=255)
                       cgccgcAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCAc                                                                                <  1:892874/101‑1 (MQ=255)
                          cgcAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTc                                                                             <  1:4243631/101‑1 (MQ=255)
                           gcAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCg                                                                            <  1:14424976/101‑1 (MQ=255)
                            cAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCgg                                                                           >  1:7904226/1‑101 (MQ=255)
                              atatGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCgggg                                                                         <  1:4911933/101‑1 (MQ=255)
                                atGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGAt                                                                       >  1:6829258/1‑101 (MQ=255)
                                  gACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTg                                                                     <  1:2536331/101‑1 (MQ=255)
                                  gACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTg                                                                     >  1:12217883/1‑101 (MQ=255)
                                      gaggagTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgc                                                                 >  1:3126944/1‑101 (MQ=255)
                                       aggagTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcg                                                                <  1:1167061/101‑1 (MQ=255)
                                        ggagTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTgcgg                                                               >  1:3012187/1‑101 (MQ=255)
                                             tAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatg                                                          >  1:6706586/1‑101 (MQ=255)
                                              agagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatga                                                         >  1:6540223/1‑101 (MQ=255)
                                               gagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatgat                                                        <  1:12949499/101‑1 (MQ=255)
                                                agCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGAtt                                                       <  1:25023/101‑1 (MQ=255)
                                                 gCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTg                                                      >  1:9034072/1‑101 (MQ=255)
                                                 gCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTg                                                      <  1:4579751/101‑1 (MQ=255)
                                                  cACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTgg                                                    >  1:5237648/1‑101 (MQ=255)
                                                     tgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgca                                                  <  1:9392035/95‑5 (MQ=255)
                                                     tgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgca                                                  >  1:6233546/7‑97 (MQ=255)
                                                      gtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcag                                                 <  1:1677956/96‑6 (MQ=255)
                                                        cACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGtt                                              >  1:2363087/1‑101 (MQ=255)
                                                         aggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggct                                              <  1:8481722/99‑9 (MQ=255)
                                                         aggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggct                                              >  1:4779581/3‑93 (MQ=255)
                                                         aCGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGttt                                             <  1:4962299/101‑1 (MQ=255)
                                                         aCGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGttt                                             <  1:13003281/101‑1 (MQ=255)
                                                          ggAAGCTCCGTTCTCCACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctg                                             <  1:13826390/100‑10 (MQ=255)
                                                          cGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                            >  1:11297447/1‑101 (MQ=255)
                                                          cGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                            <  1:7550469/101‑1 (MQ=255)
                                                          cGAAGCGCCGTTTTACACTCTGAGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                            <  1:8357490/101‑1 (MQ=255)
                                                           gAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctgg                                            >  1:3669272/1‑91 (MQ=255)
                                                           gAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctgg                                            >  1:1324261/1‑91 (MQ=255)
                                                             aGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATCTTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCt                                         <  1:10635875/101‑1 (MQ=255)
                                                             aGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCt                                         <  1:10834969/101‑1 (MQ=255)
                                                              gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTg                                        <  1:7928163/101‑1 (MQ=255)
                                                              gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTg                                        <  1:6644996/101‑1 (MQ=255)
                                                              gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTg                                        <  1:10385468/101‑1 (MQ=255)
                                                               cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                       <  1:13981008/101‑1 (MQ=255)
                                                                gcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcg                                      >  1:13777605/1‑101 (MQ=255)
                                                                gcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcg                                      >  1:1564095/1‑101 (MQ=255)
                                                                gcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcg                                      <  1:10785763/101‑1 (MQ=255)
                                                                 ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                     >  1:3609233/1‑101 (MQ=255)
                                                                 ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                     >  1:11485049/1‑101 (MQ=255)
                                                                  cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                    <  1:5863339/101‑1 (MQ=255)
                                                                    ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                  <  1:12648711/101‑1 (MQ=255)
                                                                    ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                  <  1:327250/101‑1 (MQ=255)
                                                                      ttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGc                                >  1:12703338/1‑101 (MQ=255)
                                                                      ttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGc                                <  1:10818777/101‑1 (MQ=255)
                                                                          acTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTg                            <  1:14358079/101‑1 (MQ=255)
                                                                           ctctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGt                           <  1:13437887/101‑1 (MQ=255)
                                                                            tctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGtt                          >  1:12627926/1‑101 (MQ=255)
                                                                             ctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTa                         >  1:7493228/1‑101 (MQ=255)
                                                                             ctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTa                         >  1:14632412/1‑101 (MQ=255)
                                                                             ctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTa                         >  1:1906241/1‑101 (MQ=255)
                                                                               ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                       >  1:12454123/1‑101 (MQ=255)
                                                                                gggCCGCTAACTACCGATATTGCGCCGGGCTATGACCTCTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGt                      <  1:6971747/101‑1 (MQ=255)
                                                                                 ggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTa                     >  1:11194587/1‑101 (MQ=255)
                                                                                   ccGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAAc                   >  1:14088818/1‑101 (MQ=255)
                                                                                   ccGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAAc                   >  1:7215197/1‑101 (MQ=255)
                                                                                    cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                  <  1:13748013/101‑1 (MQ=255)
                                                                                     gCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGc                 >  1:4338385/1‑101 (MQ=255)
                                                                                       tggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGccg               >  1:5614397/6‑100 (MQ=255)
                                                                                       tggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGccg               >  1:9125128/6‑100 (MQ=255)
                                                                                        ggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGccga              >  1:13590803/5‑99 (MQ=255)
                                                                                        aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa              >  1:1750371/1‑101 (MQ=255)
                                                                                         gcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCgaa             >  1:5639288/4‑101 (MQ=255)
                                                                                          caaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCgaag            >  1:8111134/3‑100 (MQ=255)
                                                                                          cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa            <  1:963155/101‑1 (MQ=255)
                                                                                          cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa            <  1:744363/101‑1 (MQ=255)
                                                                                          cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa            >  1:11847926/1‑101 (MQ=255)
                                                                                          cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa            <  1:11964530/101‑1 (MQ=255)
                                                                                           tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag           >  1:12262093/1‑101 (MQ=255)
                                                                                           tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag           >  1:13568891/1‑101 (MQ=255)
                                                                                            aCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAaga          <  1:3944732/101‑1 (MQ=255)
                                                                                             ccGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAggag         <  1:2572304/101‑1 (MQ=255)
                                                                                              cGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc        >  1:11014837/1‑101 (MQ=255)
                                                                                              cGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc        >  1:6358191/1‑101 (MQ=255)
                                                                                                atatTGCGCTGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCAt      <  1:6205552/101‑1 (MQ=255)
                                                                                                atatTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCAt      <  1:14909659/101‑1 (MQ=255)
                                                                                                 tatTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATc     <  1:240090/101‑1 (MQ=255)
                                                                                                   ttGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTg   <  1:3440937/101‑1 (MQ=255)
                                                                                                    tGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTgg  >  1:13151382/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
GCCACGTGCCGATGCAGATGATCCGCCGCAATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGG  >  NC_007779/3441888‑3442088

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.