Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,902 G→C 5.6% E226Q (GAG→CAG)  yhdZ → amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,486,9020GC 5.6% 598.4 / 18.6 179E226Q (GAG→CAG) yhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (71/98);  new base C (5/5);  total (76/103)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 187 total aligned reads displayed.

GACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCT  >  NC_002947/1486754‑1487051
                                                                                                                                                    |                                                                                                                                                     
gACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCga                                                                                                                                                      <  1:5073313/150‑1 (MQ=255)
   cTCGGCGCTCGCTCCGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGAt                                                                                                                                                   >  1:1106476/1‑150 (MQ=255)
   cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGACCCGGGTTTTCTTCATGGACAAGGGCGAGAt                                                                                                                                                   >  2:1453626/1‑150 (MQ=255)
   cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAt                                                                                                                                                   >  2:1062397/1‑150 (MQ=255)
   cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAt                                                                                                                                                   >  1:4951289/1‑150 (MQ=255)
   cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAt                                                                                                                                                   >  2:4441456/1‑150 (MQ=255)
   cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAt                                                                                                                                                   >  1:3610107/1‑150 (MQ=255)
   cTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAt                                                                                                                                                   >  2:4703236/1‑150 (MQ=255)
    tCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGAtt                                                                                                                                                  >  1:3771405/1‑150 (MQ=255)
     cGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtg                                                                                                                                                 >  1:2063013/1‑150 (MQ=255)
       gCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgtg                                                                                                                                               <  2:1147073/150‑1 (MQ=255)
       gCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATtgtg                                                                                                                                               <  1:3250621/150‑1 (MQ=255)
        cACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTgg                                                                                                                                              >  2:1557867/1‑150 (MQ=255)
         aCTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGTCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGa                                                                                                                                             <  1:3511387/150‑1 (MQ=255)
         aCTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGa                                                                                                                                             <  2:2811412/150‑1 (MQ=255)
         aCTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGa                                                                                                                                             <  2:1260519/150‑1 (MQ=255)
         aCTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGa                                                                                                                                             >  1:3040003/1‑150 (MQ=255)
         aCTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCAGGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGa                                                                                                                                             <  2:4840464/150‑1 (MQ=255)
               tccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTTGAACAgaa                                                                                                                                       >  2:3655207/2‑148 (MQ=255)
               tccGGAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAAGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTCGAACAgaa                                                                                                                                       <  1:4121435/149‑3 (MQ=255)
               cccGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGc                                                                                                                                       <  2:2247960/150‑1 (MQ=255)
                ccGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAGATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGcc                                                                                                                                      >  1:4915750/1‑150 (MQ=255)
                ccGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGcc                                                                                                                                      >  1:312435/1‑150 (MQ=255)
                ccGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGcc                                                                                                                                      >  1:3655398/1‑150 (MQ=255)
                ccGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGCGGAGCAGGcc                                                                                                                                      <  1:2962428/150‑1 (MQ=255)
                   gAAATGGTCAAGGAAGTGCTGGATACGATGGTGAGTCTTGCCGATGAGGGCATGACCATGATCTGTGTGACGCACGAAATGGGCTTTGCGCGTCAGGTTGCCAATCGTGTGATCTTCATGGATCAGGGCCAGATCGTTGAACAgaattcg                                                                                                                                   >  2:4226764/1‑144 (MQ=255)
                     aaTGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGcc                                                                                                                                 >  2:2028710/1‑150 (MQ=255)
                       tGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCgg                                                                                                                               <  2:247630/150‑1 (MQ=255)
                        ggTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGa                                                                                                                              <  1:5085756/150‑1 (MQ=255)
                         gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                                             >  2:4108354/1‑150 (MQ=255)
                         gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                                             >  2:661613/1‑150 (MQ=255)
                         gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                                             >  1:1046000/1‑150 (MQ=255)
                         gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                                             <  2:1815329/150‑1 (MQ=255)
                         gTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGAt                                                                                                                             <  1:2333352/150‑1 (MQ=255)
                          tCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATg                                                                                                                            >  2:1410538/1‑150 (MQ=255)
                          tCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATg                                                                                                                            >  2:1721907/1‑150 (MQ=255)
                          tCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATg                                                                                                                            >  2:1127132/1‑150 (MQ=255)
                            aaGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGAc                                                                                                                          >  2:4638467/1‑150 (MQ=255)
                              ggAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCTACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGActt                                                                                                                        >  2:5052772/1‑150 (MQ=255)
                              ggAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGc                                                                                                                                           >  2:2910155/1‑131 (MQ=255)
                              ggAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGActt                                                                                                                        >  2:4849567/1‑150 (MQ=255)
                              ggAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGActt                                                                                                                        >  1:2626778/1‑150 (MQ=255)
                                 aGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACttcttc                                                                                                                    <  2:1157049/150‑1 (MQ=255)
                                  gTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACttcttc                                                                                                                    <  2:1121285/150‑1 (MQ=255)
                                  gTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACttcttc                                                                                                                    <  1:3704529/150‑1 (MQ=255)
                                  gTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACttcttc                                                                                                                    <  2:1082486/150‑1 (MQ=255)
                                     cTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGAc                                                                                                                 >  1:3446633/1‑150 (MQ=255)
                                      tCGATACCAAGGTCGGCCTGGCCGAGGATGGCATGACCAAGCTCTGCGTGACCCACGAAATGGGCTATGCCCGTACCGTGGCGACCCGGGTTATCTTCATCGACAAGGGCGAGATTGTGGAGCAGGCAGCGCCGGATGACTTCTTCGAcc                                                                                                                >  2:3075017/1‑150 (MQ=255)
                                               tggtggGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCa                                                                                                       <  1:3905260/150‑1 (MQ=255)
                                               tggtggGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCa                                                                                                       <  1:1812241/150‑1 (MQ=255)
                                                        tGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGa                                                                                              <  1:3258465/150‑1 (MQ=255)
                                                        tGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGa                                                                                              <  2:206391/150‑1 (MQ=255)
                                                         ggCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGAc                                                                                             <  2:4125700/150‑1 (MQ=255)
                                                             gAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAg                                                                                         <  2:1046245/150‑1 (MQ=255)
                                                             gAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAg                                                                                         <  2:2418695/150‑1 (MQ=255)
                                                             gAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAg                                                                                         >  1:4785520/1‑150 (MQ=255)
                                                                   ggCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTc                                                                                   <  2:1940862/150‑1 (MQ=255)
                                                                   ggCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTc                                                                                   <  1:623120/150‑1 (MQ=255)
                                                                   ggCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTc                                                                                   <  1:1225325/150‑1 (MQ=255)
                                                                   ggCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTc                                                                                   <  1:1794590/150‑1 (MQ=255)
                                                                   ggCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTc                                                                                   <  1:1968194/150‑1 (MQ=255)
                                                                   ggCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTc                                                                                   >  2:1430122/1‑150 (MQ=255)
                                                                   ggCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTc                                                                                   <  1:1012265/150‑1 (MQ=255)
                                                                       tGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCa                                                                               <  1:1389500/150‑1 (MQ=255)
                                                                        gACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAg                                                                              <  1:3406254/150‑1 (MQ=255)
                                                                         aCCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGc                                                                             >  2:3110349/1‑150 (MQ=255)
                                                                         aCCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGc                                                                             <  1:4979532/150‑1 (MQ=255)
                                                                          ccATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGcc                                                                            >  1:999606/1‑150 (MQ=255)
                                                                             tGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGa                                                                         >  1:2579997/1‑150 (MQ=255)
                                                                               ctctGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATc                                                                       <  2:4773458/150‑1 (MQ=255)
                                                                                tctGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATcc                                                                      >  2:1335079/1‑150 (MQ=255)
                                                                                      tGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCAtt                                                                <  2:4152792/150‑1 (MQ=255)
                                                                                          ccACGAAATGGGCTTTGCCCGTACGGTGGCGAAGCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACGAAGTTGTTCCTCAGCCAGATCCTGCATTGATg                                                            >  1:4069706/1‑150 (MQ=255)
                                                                                          ccACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATg                                                            >  2:1609910/1‑150 (MQ=255)
                                                                                              gAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGt                                                        >  1:2737987/1‑150 (MQ=255)
                                                                                                 aTGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTg                                                     <  2:3489864/150‑1 (MQ=255)
                                                                                                 aTGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTg                                                     <  2:3490093/150‑1 (MQ=255)
                                                                                                   gggCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTggg                                                   <  1:3177221/150‑1 (MQ=255)
                                                                                                     gCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCt                                                 <  1:2617463/150‑1 (MQ=255)
                                                                                                           cccGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCgg                                           <  1:3663666/150‑1 (MQ=255)
                                                                                                             cGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCa                                         <  1:3042396/150‑1 (MQ=255)
                                                                                                                  cGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCT                                    <  1:4943495‑M1/150‑1 (MQ=255)
                                                                                                                     ggCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCCGTGCAGCCCaaac                                 <  1:3277524‑M1/150‑5 (MQ=255)
                                                                                                                      gCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCCGCGCAGCCCaaaca                                <  2:2153457‑M1/150‑6 (MQ=255)
                                                                                                                       cGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGTCTTGTTCGTCCGGCCGCGCAGTGGCCGGTCCAAGTGGTTCCGCAGCCAGATCCTGCTGTGAGGGGGTTGGGGCTGCGGGGCGGCCGTTggcg                               >  2:5023778/1‑150 (MQ=255)
                                                                                                                        gAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGG                              <  1:2244988‑M1/150‑1 (MQ=255)
                                                                                                                        gAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGG                              <  1:4438253‑M1/150‑1 (MQ=255)
                                                                                                                            cGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCAC                          >  1:4631286‑M1/1‑150 (MQ=255)
                                                                                                                              ggTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGC                        <  1:947687‑M1/150‑1 (MQ=255)
                                                                                                                                     ttCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTgcgc                                             >  1:1023272/1‑122 (MQ=255)
                                                                                                                                     ttCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCC                 <  2:3478812‑M1/150‑1 (MQ=255)
                                                                                                                                      tCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCC                >  2:2732335‑M1/1‑150 (MQ=255)
                                                                                                                                        aTGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCAC              <  2:3189806‑M1/150‑1 (MQ=255)
                                                                                                                                           gACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGG           <  1:4419040‑M1/150‑1 (MQ=255)
                                                                                                                                            aCAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGa          <  1:1609910‑M1/150‑2 (MQ=255)
                                                                                                                                              aaGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt        <  2:312435/150‑1 (MQ=255)
                                                                                                                                               aGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGaatt       >  2:2366915‑M1/1‑146 (MQ=255)
                                                                                                                                                gggCGAGATTGTGGAGCAGGCGGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCCAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAg      <  1:2882573/150‑1 (MQ=255)
                                                                                                                                                   cGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAgcgc   <  1:4264785/150‑1 (MQ=255)
                                                                                                                                                    gagaTTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCt  <  2:3040003/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACCTCGGCACTGGACCCGGAAATGGTCAAGGAAGTGCTCGATACCATGGTGGGCCTGGCCGAGGATGGCATGACCATGCTCTGCGTGACCCACGAAATGGGCTTTGCCCGTACCGTGGCGAACCGGGTTATCTTCATGGACAAGGGCGAGATTGTGGAGCAGGCCGCGCCGGATGACTTCTTCGACCGGCCGCGCAGTGACCGGACCAAGTTGTTCCTCAGCCAGATCCTGCATTGATGCTGTTTGGGCTGCGCGGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTAGCGCT  >  NC_002947/1486754‑1487051

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