Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,442,037 +C 6.5% coding (1456/1896 nt) thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,442,0371.→C 6.5% 69.2 ‑1.7 ‑0.0232coding (1456/1896 nt)thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (10/5):  ref base (101/116):  total (111/121)
Fisher's exact test strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.91e-01

Only 100 of 254 total aligned reads displayed.

GCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAT  >  NC_007779/3441937‑3442137
                                                                                                     |                                                                                                    
    tgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgca                                                                                                   <  1:9392035/95‑5 (MQ=255)
    tGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTg                                                                                                  <  1:4169675/101‑1 (MQ=255)
     gtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcag                                                                                                  <  1:1677956/96‑6 (MQ=255)
        aCGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGttt                                                                                              <  1:13003281/101‑1 (MQ=255)
        aCGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGttt                                                                                              <  1:2375757/101‑1 (MQ=255)
         ggAAGCTCCGTTCTCCACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctg                                                                                              <  1:13826390/100‑10 (MQ=255)
         cGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                                                                             >  1:7865544/1‑101 (MQ=255)
         cGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                                                                             >  1:8142549/1‑101 (MQ=255)
         cGAAGCGCCGTTTTACACTCTGAGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                                                                             <  1:8357490/101‑1 (MQ=255)
            aGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATCTTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCt                                                                                          <  1:10635875/101‑1 (MQ=255)
             gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTg                                                                                         <  1:7928163/101‑1 (MQ=255)
             gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTg                                                                                         <  1:6644996/101‑1 (MQ=255)
              cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                                                                        <  1:13981008/101‑1 (MQ=255)
               gcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcg                                                                                       >  1:1564095/1‑101 (MQ=255)
                ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGACTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                                                                      >  1:5010457/1‑101 (MQ=255)
                 cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                                     <  1:5863339/101‑1 (MQ=255)
                   ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                                                   <  1:12648711/101‑1 (MQ=255)
                   ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                                                   <  1:8485158/101‑1 (MQ=255)
                   ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                                                   <  1:1861594/101‑1 (MQ=255)
                       acacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGctc                                                                               <  1:6225141/101‑1 (MQ=255)
                          ctctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGt                                                                            <  1:11330085/101‑1 (MQ=255)
                           tctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGtt                                                                           >  1:12627926/1‑101 (MQ=255)
                           tctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGtt                                                                           >  1:10591265/1‑101 (MQ=255)
                            ctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTa                                                                          >  1:7493228/1‑101 (MQ=255)
                             tGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTAc                                                                         >  1:5262355/1‑101 (MQ=255)
                               gggCCGCTAACTACCGATATTGCGCCGGGCTATGACCTCTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGt                                                                       <  1:6971747/101‑1 (MQ=255)
                                 gCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTaa                                                                     >  1:10810026/1‑101 (MQ=255)
                                   cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                                                                   <  1:13748013/101‑1 (MQ=255)
                                   cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                                                                   <  1:2892372/101‑1 (MQ=255)
                                   cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                                                                   >  1:3729186/1‑101 (MQ=255)
                                    gCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGc                                                                  >  1:2063631/1‑101 (MQ=255)
                                      tggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGccg                                                                >  1:9125128/6‑100 (MQ=255)
                                      tggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGccg                                                                >  1:5614397/6‑100 (MQ=255)
                                      tAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCa                                                                <  1:3473687/101‑1 (MQ=255)
                                       ggcaaCCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGccga                                                               >  1:13590803/5‑99 (MQ=255)
                                       aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                                               >  1:1750371/1‑101 (MQ=255)
                                        aCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaa                                                              >  1:12376845/1‑101 (MQ=255)
                                         cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa                                                             <  1:11964530/101‑1 (MQ=255)
                                         cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa                                                             >  1:11847926/1‑101 (MQ=255)
                                          tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                                            >  1:13568891/1‑101 (MQ=255)
                                           aCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAaga                                                           <  1:3944732/101‑1 (MQ=255)
                                             cGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc                                                         >  1:11014837/1‑101 (MQ=255)
                                              gATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCa                                                        >  1:6026108/1‑101 (MQ=255)
                                                tatTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATc                                                      <  1:14722439/101‑1 (MQ=255)
                                                catcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCacc                                                      <  1:4630784/97‑3 (MQ=255)
                                                 atcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCACCt                                                     <  1:13169344/98‑1 (MQ=255)
                                                 atTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCCTCt                                                     >  1:4866426/1‑101 (MQ=255)
                                                   cgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCACCTgg                                                   >  1:734922/2‑101 (MQ=255)
                                                       ccGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCt                                               >  1:5516136/1‑101 (MQ=255)
                                                       ccGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCt                                               <  1:12010572/101‑1 (MQ=255)
                                                         gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                                             >  1:14976058/1‑101 (MQ=255)
                                                         gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                                             >  1:6617377/1‑101 (MQ=255)
                                                         gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                                             >  1:8437070/1‑101 (MQ=255)
                                                         gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                                             >  1:56414/1‑101 (MQ=255)
                                                           gCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                                           <  1:5127414/101‑1 (MQ=255)
                                                              aTGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAAt                                        <  1:253170/101‑1 (MQ=255)
                                                              aTGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAAt                                        >  1:4908967/1‑101 (MQ=255)
                                                               tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                       <  1:13323856/101‑1 (MQ=255)
                                                               tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                       <  1:12876002/101‑1 (MQ=255)
                                                                gACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa                                      >  1:5842962/1‑101 (MQ=255)
                                                                 aCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaaa                                     >  1:9729591/1‑101 (MQ=255)
                                                                 aCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaaa                                     <  1:1002995/101‑1 (MQ=255)
                                                                 aCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaaa                                     >  1:12210740/1‑101 (MQ=255)
                                                                  ccACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaag                                    <  1:4360723/101‑1 (MQ=255)
                                                                  ccACGTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaag                                    <  1:11548447/101‑1 (MQ=255)
                                                                   cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga                                   <  1:108163/101‑1 (MQ=255)
                                                                   cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga                                   <  1:7548212/101‑1 (MQ=255)
                                                                    aCTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                  >  1:5710175/1‑101 (MQ=255)
                                                                      ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGCTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                                <  1:2245020/101‑1 (MQ=255)
                                                                        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATg                              >  1:9041191/1‑101 (MQ=255)
                                                                        cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATg                              <  1:11155375/101‑1 (MQ=255)
                                                                         aCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGt                             <  1:1461488/101‑1 (MQ=255)
                                                                          cGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt                            >  1:9796102/1‑101 (MQ=255)
                                                                           gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                           >  1:10517716/1‑101 (MQ=255)
                                                                           gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                           <  1:6773012/101‑1 (MQ=255)
                                                                           gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                           <  1:4690823/101‑1 (MQ=255)
                                                                            tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCCGGGTCTGCCCAATAAAGAAGATGTTaa                          <  1:2198725/101‑1 (MQ=255)
                                                                             cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg                         <  1:3265555/101‑1 (MQ=255)
                                                                             cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg                         <  1:1927262/101‑1 (MQ=255)
                                                                                 gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTTCCCAATAAAGAAGATGTTAAGCAgg                     <  1:6310736/101‑1 (MQ=255)
                                                                                 gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgg                     <  1:11646233/101‑1 (MQ=255)
                                                                                  aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAggg                    >  1:10181271/1‑101 (MQ=255)
                                                                                  aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAggg                    <  1:3324230/101‑1 (MQ=255)
                                                                                     ggTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCt                 <  1:1406818/101‑1 (MQ=255)
                                                                                     ggTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCt                 <  1:3079141/101‑1 (MQ=255)
                                                                                        gcggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTAt              <  1:11288074/101‑1 (MQ=255)
                                                                                          ggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCa            <  1:7697594/101‑1 (MQ=255)
                                                                                          ggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCa            >  1:14363442/1‑101 (MQ=255)
                                                                                             gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCt         >  1:5357718/1‑101 (MQ=255)
                                                                                             gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCt         >  1:6886943/1‑101 (MQ=255)
                                                                                             gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCt         <  1:3606435/101‑1 (MQ=255)
                                                                                              atgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCta        <  1:12256168/101‑1 (MQ=255)
                                                                                              atgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCta        <  1:12843742/101‑1 (MQ=255)
                                                                                               tgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtat       <  1:1906920/101‑1 (MQ=255)
                                                                                                  ttGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAg    >  1:10952763/1‑101 (MQ=255)
                                                                                                    ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAt  >  1:7258629/1‑101 (MQ=255)
                                                                                                    ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAt  <  1:7066420/101‑1 (MQ=255)
                                                                                                    ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAt  >  1:10489345/1‑101 (MQ=255)
                                                                                                    ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAt  <  1:4807159/101‑1 (MQ=255)
                                                                                                    ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAt  <  1:14498200/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
GCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAT  >  NC_007779/3441937‑3442137

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.