Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,095 A→T 5.8% intergenic (+56/+96) yfhM → / ← PP_0571 lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947668,0950AT 5.8% 604.5 / 21.4 215intergenic (+56/+96)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (76/115);  new base T (5/7);  total (89/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-02

Only 100 of 217 total aligned reads displayed.

GCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTC  >  NC_002947/667948‑668243
                                                                                                                                                    |                                                                                                                                                     
  gTGCCGCCGCCGCAGGTCGAATCGAGGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGACAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCg                                                                                                                                                  <  2:7659876/150‑1 (MQ=255)
     ccgccgccgCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGcc                                                                                                                                               <  2:2514227/150‑1 (MQ=255)
       gccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCGGGCGGCGCGGTCGGGCCGCAAA‑GCGGcccc                                                                                                                                             >  1:8809914/1‑150 (MQ=255)
       gccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGcccc                                                                                                                                             <  2:3820994/150‑1 (MQ=255)
        ccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTT‑GCGGcccta                                                                                                                                            <  1:386299/150‑3 (MQ=255)
         cgccgcAGGTCGAATCGTTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTATAGACCCGCTCTACGCGT‑AAAGTGTGGGTGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGa                                                                                                                                           <  1:6893730/150‑1 (MQ=39)
         cgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTT‑GCGGccctat                                                                                                                                           <  2:1281174/150‑4 (MQ=255)
         cgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTT‑GCGGccctat                                                                                                                                           <  2:2481614/150‑4 (MQ=255)
             gcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAA‑GCGGCCCCGACCAt                                                                                                                                       >  2:6372640/1‑150 (MQ=255)
               aGGTCGAACCGATGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGCATGGGGAGTTGACGGTGAAAGCCCGATTGAGACCCGCTGTACGCGG‑AATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTc                                                                                                                                     <  1:7083518/150‑1 (MQ=38)
                   cGAATGGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGT‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAgc                                                                                                                                 <  2:2280580/150‑1 (MQ=255)
                   cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAgc                                                                                                                                 <  2:5978779/150‑1 (MQ=255)
                    gactCGATGTACCGGCCGACCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAgcg                                                                                                                                <  2:7104983/147‑1 (MQ=255)
                     aaTCGATGTACCGGCCGAACCTGCTGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAGGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGAAAA‑GCGGCCCCGACAATTCAAgcgg                                                                                                                               <  2:2244850/150‑1 (MQ=255)
                     aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTTGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAgcgg                                                                                                                               <  1:8702533/150‑1 (MQ=255)
                     aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCGGGGCGGCGCGGTGGGGCGGCCGA‑GCGGCGCCGACCAttacaccgg                                                                                                                               >  1:2503843/1‑143 (MQ=255)
                     aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCCTCGCGCTATGGCCGCCAA‑GCGGCCCCGACAAAAAAAgcgg                                                                                                                               >  2:3774189/1‑150 (MQ=255)
                      aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCTTGCGTGGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTAAAgcggc                                                                                                                              >  2:4065137/1‑150 (MQ=255)
                      aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGGGAGGGCCGCAAA‑GCGGCCCCGACAAATAAAgcggc                                                                                                                              >  2:848416/1‑150 (MQ=255)
                          aTGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                                          <  1:7896678/150‑1 (MQ=255)
                          aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCCAA‑GCGGCCCCGACAAGTCCAGCGGCGAAg                                                                                                                          >  2:8142396/1‑150 (MQ=255)
                            gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                                        <  2:654187/150‑1 (MQ=255)
                            gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                                        <  2:7892158/150‑1 (MQ=255)
                            gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                                        <  2:7318515/150‑1 (MQ=255)
                             tACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGCAGCTg                                                                                                                       >  2:3343192/1‑150 (MQ=255)
                             gaCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTg                                                                                                                       <  2:5262133/149‑1 (MQ=255)
                              aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCGGa                                                                                                                      >  2:452238/1‑150 (MQ=255)
                              aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCCCAAA‑GCCGCCCCGACAATACAATCGgcggagcgga                                                                                                                      >  2:83580/1‑150 (MQ=255)
                               ccGGCCGAACCTGCAGGCCGTGGGCGTTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGATACCCGCTCTACGCGGAAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGa                                                                                                                      <  2:2815495/150‑1 (MQ=255)
                                cGGCCGAACCTGCTGGCCGTGGGCGATGGGCAAGGGGAGTTGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGT‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                                    <  1:5537959/150‑1 (MQ=255)
                                cGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                                    <  2:7815068/150‑1 (MQ=255)
                                 ggCCGAACCTGCTGGCCGTGGGCGTTGGGCATGGGGAGATGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGT‑AATGTGTGGGTGCGGCCTTGCGTCGCTTTAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAg                                                                                                                   <  1:5183488/150‑1 (MQ=255)
                                  gCCGACCCTGCAGGCGGGGGGGGTTGGGCATGGGGAGTTGACGGTGAAAGCCCGCTCGCGACCCGCTCTAGGGGG‑AATGTGGGGGTGGGGCCTTGCGGCGGTATCGGGCCGCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGc                                                                                                                  <  1:2989842/150‑1 (MQ=255)
                                     gAACCTGCAGGCCGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                                               <  1:6535772/150‑1 (MQ=255)
                                     gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGGTAGGGCCGCCAA‑GCGGCCCCGACAATTGAAGCGGCGTCGCGGAAAGccgg                                                                                                               >  1:1370219/1‑148 (MQ=38)
                                      aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                                              <  2:1913650/150‑1 (MQ=255)
                                        ccggccggCCGTGGGCGTTGGGCATGGGGAGAGGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGG‑AATGTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAAC‑GGGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgggg                                                                                                            <  1:1497149/144‑1 (MQ=37)
                                        ccTGCTGGCGGTGGGCGAGGGGCATGGGGAGTTGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGGGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgggg                                                                                                            <  1:4480142/150‑1 (MQ=37)
                                        ccTGCTGGCCGTGGGCGTGGGGCATGGGGTGCTGACGGTGAAAGCCCGCTAGAGCCCCGCTCTAGGGGG‑AATGTGTGGGCGCGGCCTTGCGTCGCGATTGGGCCGCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgggg                                                                                                            <  1:965880/150‑1 (MQ=255)
                                           gCAGGCCGTGGGCGATGGGCCAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCGCCAA‑GCGGCCCCGACCATTCAAGCGGCGGCGCGGGAGGCCGGGGgcgg                                                                                                         >  2:5404525/1‑148 (MQ=255)
                                           gCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCCAA‑GCGGCCCCCAAAAAACAAGCGGCGAACCGGAAAGCCGGGGgctg                                                                                                         >  2:1782315/1‑150 (MQ=38)
                                             aGGCGGTGGGCGAGGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                                       <  2:2185762/150‑1 (MQ=255)
                                             aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGGGCGGCGCGCTTGGGCCGCCAA‑GCGGCCCCGCCCCGTGACGCGGCGAAGCGGAAAGCCTGGGgcggct                                                                                                       >  1:3617728/1‑150 (MQ=255)
                                               gCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGGTATGGGCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCGGGGGCGGCggt                                                                                                     >  2:8224934/1‑150 (MQ=255)
                                                ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCGGGCGGGGCGGTAGGGCCGCCCC‑GCGGCCCCGCCCAAAAAAGCCGCGCAGCGGCGAGCCGGGggcggcgggg                                                                                                    >  1:2715387/1‑146 (MQ=255)
                                                 cGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGc                                                                                                   <  2:573335/150‑1 (MQ=255)
                                                  gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                  <  1:8618762/150‑1 (MQ=255)
                                                  gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCCACAAATCAAGCCGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                  >  2:8515696/1‑150 (MQ=255)
                                                    gggggttGGGCATGGGGAGATGCCGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGT‑AATGTGTGGGTGCGGCCTTGCGTCGCTATTGGGCCCCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                <  1:3690017/144‑1 (MQ=255)
                                                    gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCAG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                <  1:3557094/150‑1 (MQ=255)
                                                     ggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                               <  2:7058323/150‑1 (MQ=255)
                                                     ggCGAGGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                               <  2:4800484/150‑1 (MQ=255)
                                                       cgtggggCTAGGGGCGAGTAGGGTGAACGCCCGCTCCCGCCCCGCTCTACGCGT‑AACGTGTGGGTGCGGCCTTGCGTCGCGATTGGGCCGCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                             <  1:275128/146‑1 (MQ=255)
                                                       cgcggggCACGGGGAGATGCCGGTGAACGCCCGCTCGAGCCCCGCTCTACGCGG‑AACGGGGGGGAGCGGCCTGGCGTCGCGATCGGGCCGCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                             <  1:2540664/146‑1 (MQ=255)
                                                       cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGT‑AAAGTGTGGGAGCGGCCTTGCGTCGCTATAGGCCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                             <  1:116934/150‑1 (MQ=255)
                                                       cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGCAAGGGCCGCAAA‑GCGGCCCCGACAAAAAAACCGGCGAAGCCGAAACCCAGGCGCTGCTGAGCAGacca                                                                                             >  1:7275898/1‑149 (MQ=37)
                                                           gggCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCGTGCGGCGGGTTAGGGCCGCAAA‑GCGGCCCCGTCAAAGGAAGCGGCGAAGGGGAAAGCCGGGGGCTGCTGTGGAGcccaggcg                                                                                         >  1:6464384/1‑145 (MQ=255)
                                                              cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCAGGCGTGGCGATAGGGCCGCCAA‑GCGGCCCCGACCATACAAGCGGCCAAGCGGTACGCCAGGGGCGGCGGTGCCGCCCAGGCGCGa                                                                                      >  1:1027646/1‑150 (MQ=255)
                                                              cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                                      >  2:8052571‑M2/1‑144 (MQ=255)
                                                               aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCCTCGCGATAGGGCCGCCAA‑GCGGCCCCCACAAATCAAGCGGCGAACCGGAAAGCCTGGGGCTGCAGTGCAGCCCGGtcgcgac                                                                                     <  2:4440916‑M2/150‑8 (MQ=255)
                                                                    gagaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaag                                                                                >  2:2995081‑M2/1‑138 (MQ=255)
                                                                     tgaggccggtgaaAGCCCGCTAGAGACCCGCTCTCCGCGG‑AAAGTGTGGGAGGGGCCTTGCGTCGCGATCGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                               >  1:946702‑M2/7‑137 (MQ=255)
                                                                      gttgaGGGTCAACGCCCGCTATACACCGGCTCTACGGGT‑AACGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggc                                                                              >  1:1751967‑M2/3‑136 (MQ=255)
                                                                         gACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAA‑GCGGCCCCGACAAATAAAGCGGCGAAGCTGAAAGCCAGGGGCTGCGGTGCAGCCCAGtcgcgacgccaggccgc                                                                           <  2:2588190‑M2/150‑18 (MQ=255)
                                                                            ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGCGCGGCCGTGCGGCGCGCTAGGGCCGCACA‑GCGGCCCCGACACGGCCCGCGGCGCAGCGGCAAGCCGGGGGCGGCGGGGCCGCCCTGGCGCGGCGCCCGGCCGCGcc                                                                        >  2:5098460/1‑150 (MQ=255)
                                                                            ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCCAGCGGCGAAGCTGTAAGCCGGGGGCGGCTTTGCAGCCCTGtcgcgacgccaggccgcgcc                                                                        <  2:5209047‑M2/150‑21 (MQ=255)
                                                                                  aGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCACA‑GCGGCCCCGACAATTAAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGGGCAGCCCAGtcgcgacgcaaggccgctcccaccaa                                                                  <  2:2039686‑M2/150‑27 (MQ=255)
                                                                                    cccGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAA‑GCGGCCCCGACAAATCAAGCGGCGAAGCGGAAAGCCAGGGGCTGCTGGGCAGCCCGGTCGCGCCGCAAGGCCGCGCCCACCAAGa                                                                >  2:8624150/1‑150 (MQ=37)
                                                                                        cTAGAGACCCGCTCTCCGGGT‑AATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAC‑GCGGCCCCGAAAATTCACGCGGCGAAGCTGAACGCCTGGGGCTGCTGTGCAGCCCTGTGGCGCCGCAAGGCCGCTCCCACAAAGACcgc                                                            <  2:5831604/150‑1 (MQ=25)
                                                                                           gagaCCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCCACAAAAAAAGCGGCGAAGCGGGAAGCCAGGGGCTGCGGTGCAGCCCAGTCGCGACGCAAGGCCGCGCCCACAAACACcgcccc                                                         >  2:2923815/1‑147 (MQ=25)
                                                                                             caCCCGCTCTCCGCGG‑ACTGTGGGGGCGCGGCCTTGCGTCGCGATGGGGCCGAAAC‑GCGGCCCCGACACTTCACGCGGCGACGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGTAGa                                                       <  2:8242546/149‑1 (MQ=25)
                                                                                              aCCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAA‑GCGGCCCCGACAAATCAAGCGGCGGAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCAGTCGCGACGCAAGGCCGCTCCCCCAAAGACCgcgccgct                                                      >  2:3398330/1‑146 (MQ=37)
                                                                                                        cgcgG‑AAAGTGTGGGAGCGGCCTTGCCTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAAATCAACCGGCCAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCCAGGCCCCTCCCACCAAGACCGCGCCGATCAGTGGATCa                                            >  2:8727765/1‑150 (MQ=37)
                                                                                                          cgG‑AACGTGGGGGTGCGGCCTTGCGTCGCGATCGGGCCGCAAC‑GCGGCCCCGACAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCAcc                                          <  2:7935680/150‑1 (MQ=37)
                                                                                                                 ttgtgGGTGCGCCCTTGCGAGCCTATTGGCCCCCAAC‑GCGCCCCCCACAATACAAGCGGCGAAGCTTAAAGCCTGGGGCTGCTGTGCCGCCCTGTCGCCACGAACGCCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCtc                                    <  2:8293367/149‑1 (MQ=255)
                                                                                                                    tgGGCGGGGCCTTGGGTCGCTTTCGGGCCGCAAC‑GGGGCCCCGACAATTCACGCGGCCAAGCTGAAAGCCTGGGGCTGCTTTGCAGCCCTTTGGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCtcatc                                 <  2:1359164/150‑1 (MQ=25)
                                                                                                                       gtggggCGTTGCGTCGCGAT‑GGGCCCCAAACGCGCCCCCGACAATTCACGCGGCGATGCTTACAGCGTGGGGGTGGTGGGCCGCCCTGCCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAcc                              <  2:4634528/146‑1 (MQ=255)
                                                                                                                       gcgCGGCCTTGCGCCGCTATCGGGCCGAAAC‑GCGGCCCCGACAATTAACGCGGCGACGCTGAACGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGAAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAcc                              <  2:4346061/148‑1 (MQ=37)
                                                                                                                       gAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                               <  2:5696415/85‑1 (MQ=255)
                                                                                                                       gAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                               >  1:5696415/1‑85 (MQ=255)
                                                                                                                        tgCGGCCTTGGGCCGCGATTGGGCCGCAAC‑GCGGCCCCTAAAATACACGCGGCGAAGCTCAAAGCCTGGGGCTGCTGTGCCCCCCTGTCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAccc                             <  2:8284482/149‑1 (MQ=255)
                                                                                                                        aGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAccc                             >  2:8412343/1‑150 (MQ=255)
                                                                                                                            gCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatc                         <  1:6001027‑M2/150‑68 (MQ=255)
                                                                                                                               ttGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc                      >  1:8481763/1‑150 (MQ=255)
                                                                                                                                 gCGTCGCTATAGGCCCGCAAA‑GCGCCCCCCAAAATACAACCGGCTAAGCTGAAAGCCTGGGCCTGCTGTGCACCCCTGTCCCGACGAAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCa                    <  2:4925876/150‑1 (MQ=255)
                                                                                                                                   gTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGc                  >  1:5825507/1‑150 (MQ=255)
                                                                                                                                    tCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttg                 <  1:4329158‑M2/150‑76 (MQ=255)
                                                                                                                                         aTAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagc            >  1:3770202/1‑150 (MQ=255)
                                                                                                                                           ggggCCGCTTT‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca          >  1:6206114/2‑150 (MQ=255)
                                                                                                                                           ggggCCGCTTT‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca          >  2:7413528/2‑150 (MQ=255)
                                                                                                                                           ggggCCGCTTT‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca          >  2:1695873/2‑150 (MQ=255)
                                                                                                                                           aGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca          >  1:8225887/1‑150 (MQ=255)
                                                                                                                                             ggCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgcccatcgcc        <  1:6164513‑M2/150‑85 (MQ=255)
                                                                                                                                             ggCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggct                                     >  1:5246947‑M2/1‑66 (MQ=255)
                                                                                                                                             ggCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcacacactttccgcgtagagcgggtctctagagggctttcaccgtcatctccccttgcccatcgcc        <  2:8784541‑M2/150‑85 (MQ=255)
                                                                                                                                              gCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgcccatcgccc       <  2:4264194‑M2/150‑86 (MQ=255)
                                                                                                                                              gCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgat       >  1:2122464/1‑150 (MQ=255)
                                                                                                                                              gCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGAGCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgat       >  1:3224090/1‑150 (MQ=255)
                                                                                                                                               ccGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatg      >  1:7526909/1‑150 (MQ=255)
                                                                                                                                                   aaa‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTc  >  1:2964616/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTC  >  NC_002947/667948‑668243

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