Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,867,789:1 +C 16.8% intergenic (‑28/+24) thyA ← / ← lgt thymidylate synthase/prolipoprotein diacylglyceryl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,867,7891.C16.8% 304.5 / 70.0 101intergenic (‑28/+24)thyA/lgtthymidylate synthase/prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base . (39/45);  new base C (9/8);  total (48/53)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 103 total aligned reads displayed.

CCCTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTT  >  NC_002947/5867640‑5867937
                                                                                                                                                      |                                                                                                                                                    
cccTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAgggg                                                                                                                                                       >  1:2313290/1‑150 (MQ=255)
  cTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGa                                                                                                                                                    >  1:2168541/1‑150 (MQ=255)
     cTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCTCTGCAACGTGCCGTTATTAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCg                                                                                                                                                 <  1:463966/150‑1 (MQ=255)
          ggTCGAAACGCAGCATCGCGACCGGCAGGCTGATAGTGCGGATGCCGGTGCGGTTTCACTGCAACGTGCCGTTATCAATGACGTTGCGGATCAAGTCCAGATACTTTTTCATGGATAACAAGTAAAAAAAACGGCAGGGGGATCGCCCCt                                                                                                                                            <  2:1032192/150‑1 (MQ=255)
           gTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTg                                                                                                                                           >  2:1033740/1‑150 (MQ=255)
             cGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGcc                                                                                                                                         <  2:169534/150‑1 (MQ=255)
             cGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGcc                                                                                                                                         <  2:2313290/150‑1 (MQ=255)
             cGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGcc                                                                                                                                         <  2:2175148/150‑1 (MQ=255)
               aaaCGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGt                                                                                                                                       >  1:1487569/1‑150 (MQ=255)
                 aCGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACTTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTgg                                                                                                                                     <  1:1033740/150‑1 (MQ=255)
                  cGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGCCGtttt                                                                                                                                    <  1:256595/150‑4 (MQ=255)
                       caTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGCCGttttt                                                                                                                                   >  2:86224/1‑142 (MQ=255)
                       caTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGCCGttttt                                                                                                                                   <  1:86224/146‑5 (MQ=255)
                        aTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAAt                                                                                                                              <  1:2235477/150‑1 (MQ=255)
                             gcCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGc                                                                                                                         <  1:129690/150‑1 (MQ=255)
                               ccGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTAACTCTAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCgg                                                                                                                       <  2:2372722/150‑1 (MQ=255)
                                   cAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGttt                                                                                                                   >  1:1798004/1‑150 (MQ=255)
                                      gCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGc                                                                                                                >  1:1212694/1‑150 (MQ=255)
                                       cTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGcc                                                                                                               <  2:1212694/150‑1 (MQ=255)
                                          aTGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTa                                                                                                            <  1:274593/150‑1 (MQ=255)
                                               gCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCtt                                                                                                       <  2:2206121/150‑1 (MQ=255)
                                                    tGCGGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGTCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGAt                                                                                                  <  2:458783/150‑1 (MQ=255)
                                                         gTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCTGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             <  2:1487569/150‑1 (MQ=255)
                                                          tGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAg                                                                                            <  1:26197/150‑1 (MQ=255)
                                                           gCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGc                                                                                           <  2:2168541/150‑1 (MQ=255)
                                                              gTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaaacaaa                                                                                        >  2:11447/1‑145 (MQ=255)
                                                              gTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCCTACCTGCAGCAAACACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAACCcca                                                                                        >  1:91846/1‑150 (MQ=255)
                                                                ttCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGt                                                                                                                                       >  1:1518918/1‑101 (MQ=255)
                                                                ttCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccacc                                                                                      >  1:1827137/1‑150 (MQ=255)
                                                                 tCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccacca                                                                                     <  2:1213759/150‑1 (MQ=255)
                                                                   ccTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGa                                                                                   >  2:2039107/1‑150 (MQ=255)
                                                                    cTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGa                                                                                   >  1:2058089/1‑149 (MQ=255)
                                                                    cTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGa                                                                                   <  2:2058089/149‑1 (MQ=255)
                                                                     tGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATc                                                                                 <  2:1488547/150‑1 (MQ=255)
                                                                     tGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATc                                                                                 >  1:2397403/1‑150 (MQ=255)
                                                                      gCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGt                                                                                                                                       >  1:1102006/1‑95 (MQ=255)
                                                                      gCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGt                                                                                                                                       <  2:1102006/95‑1 (MQ=255)
                                                                             gCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGcccc                                                                                                                                             <  2:2029547/82‑1 (MQ=255)
                                                                             gCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGcccc                                                                                                                                             >  1:2029547/1‑82 (MQ=255)
                                                                             gCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                         <  2:1827137/150‑1 (MQ=255)
                                                                              ccGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGAAAAAATGGCAGGGGGATCGTCCTTTCCGTGGGAATAACCATTCGAGTTTCTCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAgg                                                                        <  2:337070/150‑1 (MQ=255)
                                                                                 ttATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAg                                                                                                           <  2:588109/112‑1 (MQ=255)
                                                                                 ttATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAg                                                                                                           >  1:588109/1‑112 (MQ=255)
                                                                                      aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTgg                                                                                                                                     >  2:220724/1‑81 (MQ=255)
                                                                                      aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTgg                                                                                                                                     <  1:220724/81‑1 (MQ=255)
                                                                                      aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGAc                                                                <  1:1516923/150‑1 (MQ=255)
                                                                                      aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGAc                                                                >  1:1589859/1‑150 (MQ=255)
                                                                                       aTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCg                                                                                                                 >  2:1235654/1‑100 (MQ=255)
                                                                                       aTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCg                                                                                                                 <  1:1235654/100‑1 (MQ=255)
                                                                                        tGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGGAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTa                                                              >  1:211836/1‑150 (MQ=255)
                                                                                            gTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcat                                                          <  1:2039107/150‑1 (MQ=255)
                                                                                                  gACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAAt                                                    >  1:389249/1‑150 (MQ=255)
                                                                                                        gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGa                                                                          >  1:108041/1‑122 (MQ=255)
                                                                                                        gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGa                                                                          <  2:108041/122‑1 (MQ=255)
                                                                                                        gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAg                                              >  1:1934733/1‑150 (MQ=255)
                                                                                                         tCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGc                                             >  1:690324/1‑150 (MQ=255)
                                                                                                            aGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCt                                                                                                                                            <  1:1305323/52‑1 (MQ=255)
                                                                                                            aGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCt                                                                                                                                            >  2:1305323/1‑52 (MQ=255)
                                                                                                             gATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAg                                                                                            >  2:886372/1‑99 (MQ=255)
                                                                                                             gATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAg                                                                                            <  1:886372/99‑1 (MQ=255)
                                                                                                              aTACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTg                                        <  2:389249/150‑1 (MQ=255)
                                                                                                                     ttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                                 <  2:690324/150‑1 (MQ=255)
                                                                                                                     ttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                                 <  2:2397403/150‑1 (MQ=255)
                                                                                                                         tGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGccc                                                                                                                                              <  2:2134108/37‑1 (MQ=255)
                                                                                                                         tGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGccc                                                                                                                                              >  1:2134108/1‑37 (MQ=255)
                                                                                                                         tGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTg                                                                                                                                      <  1:1893081/45‑1 (MQ=255)
                                                                                                                         tGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTg                                                                                                                                      >  2:1893081/1‑45 (MQ=255)
                                                                                                                         tGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACATCAACTGGCCTtgagtga                             <  1:11447/150‑1 (MQ=255)
                                                                                                                             ttACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             <  2:478753/82‑1 (MQ=255)
                                                                                                                             ttACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             >  1:478753/1‑82 (MQ=255)
                                                                                                                               aCCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGagccagc                       >  2:2396814/1‑150 (MQ=255)
                                                                                                                                 cTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGNGATTAATCANGCGGATTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGCTCAAGACGAGGCCACCGACTAGCATCGTCATGCACAGCAACTGGCCTTGAGTGAGCccgccg                     >  1:1001545/1‑150 (MQ=255)
                                                                                                                                 cTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATAGGAATGCACAGCAACTGGCGTTGAGTGAGCCAGCCg                     >  1:351929/1‑150 (MQ=255)
                                                                                                                                      gCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGc                                                                                           <  1:612486/75‑1 (MQ=255)
                                                                                                                                      gCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGc                                                                                           >  2:612486/1‑75 (MQ=255)
                                                                                                                                      gCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGa                                                                                   <  1:603387/83‑1 (MQ=255)
                                                                                                                                      gCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGa                                                                                   >  2:603387/1‑83 (MQ=255)
                                                                                                                                       cAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAcgcg               >  2:2254290/1‑150 (MQ=255)
                                                                                                                                        tcccaCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa              >  1:1199815/5‑150 (MQ=255)
                                                                                                                                        aaaaaCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa              <  1:766107/150‑1 (MQ=255)
                                                                                                                                         cccaCGGCAGGGGCGATCCCCCTGCCGTGGGCTTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCACAACAGCTCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGAt             >  1:1857014/4‑150 (MQ=255)
                                                                                                                                         aaaaCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                        <  1:198768/107‑1 (MQ=255)
                                                                                                                                         aaaaCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                        >  2:198768/1‑107 (MQ=255)
                                                                                                                                          ccaCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagtga                             <  1:1026781/131‑1 (MQ=255)
                                                                                                                                          ccaCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagtga                             >  2:1026781/3‑133 (MQ=255)
                                                                                                                                          ccaCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATg            >  2:1215484/3‑150 (MQ=255)
                                                                                                                                           caCGGCAGGGGCGATCGCCCTGCCGTGGGATTAATCAGGAGGGAATGTTCGTAGGCTTGCGATTATAAGCCCCCCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGt           >  2:695768/2‑150 (MQ=255)
                                                                                                                                           caCGGCAGGGGCGATCCCCCTTCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCAACAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGGGCCAGACGAACGCGATGt           >  2:720751/2‑150 (MQ=255)
                                                                                                                                           caCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGa                                                                                   >  2:1210857/2‑78 (MQ=255)
                                                                                                                                           caCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGa                                                                                   <  1:1210857/77‑1 (MQ=255)
                                                                                                                                           caCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGt           >  2:2040059/2‑150 (MQ=255)
                                                                                                                                            aCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGt           >  2:1764368/1‑149 (MQ=255)
                                                                                                                                            aCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGt           <  1:1764368/149‑1 (MQ=255)
                                                                                                                                              ggCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAc        <  1:2396814/150‑1 (MQ=255)
                                                                                                                                               gCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                                 <  2:122747/124‑1 (MQ=255)
                                                                                                                                               gCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                                 >  1:122747/1‑124 (MQ=255)
                                                                                                                                                 aGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGaa                                                     >  2:1630932/1‑102 (MQ=255)
                                                                                                                                                 aGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGaa                                                     <  1:1630932/102‑1 (MQ=255)
                                                                                                                                                 aGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGa     <  2:1798004/150‑1 (MQ=255)
                                                                                                                                                    gggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGtt  >  2:1128239/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCCTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTT  >  NC_002947/5867640‑5867937

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: