Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,074 T→G 11.8% intergenic (+35/+117) yfhM → / ← PP_0571 lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947668,0740TG11.8% 596.7 / 8.5 235intergenic (+35/+117)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (53/144);  new base G (31/5);  total (86/149)
Fisher's exact test for biased strand distribution p-value = 2.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

GCAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCC  >  NC_002947/667926‑668223
                                                                                                                                                     |                                                                                                                                                      
gCAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTT                                                                                                                                                       <  1:5680868‑M1/150‑1 (MQ=255)
gCAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTT                                                                                                                                                       <  2:3409068‑M1/150‑1 (MQ=255)
 cAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGTGAGCGGCCTTG                                                                                                                                                      <  2:7004958‑M1/150‑1 (MQ=255)
 cAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTG                                                                                                                                                      <  1:2367563‑M1/150‑1 (MQ=255)
  aGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGC                                                                                                                                                     >  2:1534699‑M1/1‑150 (MQ=255)
    tGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCTTTGCGg                                                                                                                                                   <  2:7751653‑M1/150‑2 (MQ=255)
       ccccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGC                                                                                                                                                >  2:5680352‑M1/1‑150 (MQ=255)
         ccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGA                                                                                                                                              <  2:4075356‑M1/150‑1 (MQ=255)
          cGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCGGGCGTCGCGgt                                                                                                                                             >  1:6961185‑M1/1‑148 (MQ=255)
            ggCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAg                                                                                                                                           <  1:6866641/150‑1 (MQ=255)
            ggCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAg                                                                                                                                           <  1:2086792/150‑1 (MQ=255)
            ggCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAg                                                                                                                                           <  2:7277032/150‑1 (MQ=255)
                ccTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGcc                                                                                                                                       <  2:5360277/150‑1 (MQ=255)
                 cTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCg                                                                                                                                      >  1:3363275/1‑150 (MQ=255)
                  tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGc                                                                                                                                     >  2:4464813/1‑150 (MQ=255)
                  tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGGGATGGGGCCgg                                                                                                                                     >  1:7251569/1‑149 (MQ=255)
                       cgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑Gc                                                                                                                               <  2:8787190/150‑1 (MQ=255)
                       cgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑Gc                                                                                                                               <  2:5750114/150‑1 (MQ=255)
                        gTGCCGCCGCCGCAGGTCGAATCGAGGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGACAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCg                                                                                                                              <  2:7659876/150‑1 (MQ=255)
                            cgccgccgcAGGTCGACTCGAGGTACCGGCCGAACCTGCAGGCCGTGGGCGCTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGccc                                                                                                                          <  2:5708047/150‑1 (MQ=255)
                            cgccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGccc                                                                                                                          <  2:4759606/150‑1 (MQ=255)
                             gccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTT‑GCGGccct                                                                                                                         <  2:832920/150‑2 (MQ=255)
                              ccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAA‑GCGGCCCCg                                                                                                                        >  2:2152362/1‑150 (MQ=255)
                                gccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTT‑GCGGccctatc                                                                                                                      <  1:3843292/150‑5 (MQ=255)
                                gccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGAGAGGGCCGCTTT‑GCGGccctatc                                                                                                                      <  2:1206498/150‑5 (MQ=255)
                                   gcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAA‑GCGGCCCCGACCAt                                                                                                                   >  2:6372640/1‑150 (MQ=255)
                                     aGGTCGAACCGATGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGCATGGGGAGTTGACGGTGAAAGCCCGATTGAGACCCGCTGTACGCGG‑AATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTc                                                                                                                 <  1:7083518/150‑1 (MQ=38)
                                       gTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGAAAGGGCAGATCACGGTGAAAGCCCGCTAGACACCCGCTCTACGCGC‑AAAGTCTGGGAGCGGCCTTGCGTCCCGATAGGCCCGCAAA‑GCGGCCCCGACAATTCaa                                                                                                               <  2:8458729/150‑1 (MQ=255)
                                          gactCGATGTACCGGCCGACCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAgcg                                                                                                            <  2:7104983/147‑1 (MQ=255)
                                          gAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGGCGCGCGAGGGCCGCACA‑GCGGCCCCCACCCATCAAgcg                                                                                                            >  2:3100447/1‑150 (MQ=255)
                                          gAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCCAA‑GCGGCCCCGACAATTCAAgcg                                                                                                            >  2:8521170/1‑150 (MQ=255)
                                           aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCGTGCGGGGCGGTTGGGCCCCCAA‑CCGGCCCCCCCCAATCAAgccg                                                                                                           >  1:7469873/1‑148 (MQ=255)
                                           aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCGGGGCGGCGCGGTGGGGCGGCCGA‑GCGGCGCCGACCAttacaccgg                                                                                                           >  1:2503843/1‑143 (MQ=255)
                                           aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAgcgg                                                                                                           >  2:8205307/1‑150 (MQ=255)
                                            aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGTTGGGAAAGGGTATATGAGGGTTAAAGCCCGCTTGAGACCCGCTCTACGGGG‑AAAGTGTGGTAGGGGCCTTGCGTCGCTATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAgcggc                                                                                                          <  1:6162445/150‑1 (MQ=37)
                                            aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAC‑GCGGCCCCGACAATTCAAgcggc                                                                                                          <  1:6523997/150‑1 (MQ=255)
                                            aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCGCAAA‑GCGGCCCCGACCATTCAAgcggc                                                                                                          >  2:848578/1‑150 (MQ=255)
                                                ttGTACCGGCCGAACCTGCAGGCGGTGGGCGTTGGGCATGGGTATATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                      <  1:1351867/149‑1 (MQ=38)
                                                  gTACCGGCCGAACCTGCAGGCGGGGGGCGAGGGGCCTGGGGATATGCGGGTGACCGCCCGCTAGCGACCCGCTCTCCGCGG‑AAAGTGGGGGCGGGGCCTTGCGTCGCGATCGGGCCGAAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                    <  1:6611502/150‑1 (MQ=255)
                                                  gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                    <  2:654187/150‑1 (MQ=255)
                                                   gaCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTg                                                                                                   <  2:5262133/149‑1 (MQ=255)
                                                    aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCGTGCGGCGCGATAGGGCCGCCCA‑GCGGCCCCGACCAATCAAGCGGCGAAGCGGa                                                                                                  >  2:5875168/1‑150 (MQ=255)
                                                    aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCGGa                                                                                                  >  2:452238/1‑150 (MQ=255)
                                                     ccGGCCGAACCTGCAGGCCGTGGGCGTTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGATACCCGCTCTACGCGGAAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGa                                                                                                  <  2:2815495/150‑1 (MQ=255)
                                                     ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCCAA‑GCGGCCCCCACAATTCAAGCGGCGAAGCGGaa                                                                                                 >  2:2207322/1‑150 (MQ=255)
                                                      cGGCCGAACCTGCTGGCCGTGGGCGATGGGCAAGGGGAGTTGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGT‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                <  1:5537959/150‑1 (MQ=255)
                                                           gacccTGCAGGCCTTGGGCTTTGGGAAAGGGTTGATGACGGTTAAAGCCCTCTAGAGACCCGCTCTCCGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCTATTGGGCCCAAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                           <  1:7714657/147‑1 (MQ=37)
                                                           gAACCTGCAGGCCGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                           <  1:8803497/150‑1 (MQ=255)
                                                            aaCCTGCAGGCGGTGGGCGTTGGGAATGGGTATTTGACGTTTAAAGCCCGCTTGAGACCCGCTCTACGCGT‑AAAGTGTGGGAGCGGCGTTGCGTCGCGATTGGGCCGAAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                          <  1:6182187/150‑1 (MQ=255)
                                                            aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                          <  2:1913650/150‑1 (MQ=255)
                                                              ccTGCTGGCGGTGGGCGAGGGGCATGGGGAGTTGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGGGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgggg                                                                                        <  1:4480142/150‑1 (MQ=37)
                                                                   aGGCCGTGGGCGATGGGCAAGGGGAGATGACTGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCCGCGCGATAGGGCCGCAAA‑GCCGCACCGACAAATCAAGCGCCCAACCAGACAGACTGGGgcagca                                                                                   >  2:7452607/1‑149 (MQ=255)
                                                                   aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGGGCGGCGCGCTTGGGCCGCCAA‑GCGGCCCCGCCCCGTGACGCGGCGAAGCGGAAAGCCTGGGgcggct                                                                                   >  1:3617728/1‑150 (MQ=255)
                                                                     gCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGGTAGGGCCGCCAA‑GCGGCCCCGACAAATCAAGCGGCGCAGCGGACCGCCTGGggcggcggg                                                                                 >  2:8225191/1‑147 (MQ=38)
                                                                      ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCGGGCGGGGCGGTAGGGCCGCCCC‑GCGGCCCCGCCCAAAAAAGCCGCGCAGCGGCGAGCCGGGggcggcgggg                                                                                >  1:2715387/1‑146 (MQ=255)
                                                                      ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGGCGCGGTAGGGCCGCAAA‑GCGGCCCCGCCAAATCAAGCGGCGCAGCGGTAAGCCGGGGGCTGCggtg                                                                                >  2:2363085/1‑150 (MQ=255)
                                                                        gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                              <  1:8618906/150‑1 (MQ=255)
                                                                        gTGGGCGATGGGCAAGGGGAGATGAAGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGAGTGGGAGCGGCCATGCCGCGCCATAGCGCCCCCCA‑GCCGCCCCCACAAATCAACCGGCGCACCGGCAACCCTGGGGCGGCtgtgaa                                                                              >  2:3805193/1‑148 (MQ=255)
                                                                          gggggttGGGCATGGGGAGATGCCGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGT‑AATGTGTGGGTGCGGCCTTGCGTCGCTATTGGGCCCCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                            <  1:3690017/144‑1 (MQ=255)
                                                                          gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggc                                                                             <  2:3372959‑M1/150‑66 (MQ=255)
                                                                          gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggc                                                                             <  2:4264194‑M1/150‑66 (MQ=255)
                                                                          gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCAG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                            <  1:3557094/150‑1 (MQ=255)
                                                                           ggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggcc                                                                            <  1:4078085‑M1/150‑67 (MQ=255)
                                                                             ggtggggCACGGGGTGATGACGTTGAAAGCCCGCTTGAGACCCGCTCTACGCGG‑AACGTGTGGGCGCGGCCTTGCGTCGCGATCGGGCCGAAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                         <  1:2990447/146‑1 (MQ=255)
                                                                             cgtggggCTAGGGGCGAGTAGGGTGAACGCCCGCTCCCGCCCCGCTCTACGCGT‑AACGTGTGGGTGCGGCCTTGCGTCGCGATTGGGCCGCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                         <  1:275128/146‑1 (MQ=255)
                                                                             cgcggggCACGGGGAGATGCCGGTGAACGCCCGCTCGAGCCCCGCTCTACGCGG‑AACGGGGGGGAGCGGCCTGGCGTCGCGATCGGGCCGCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                         <  1:2540664/146‑1 (MQ=255)
                                                                              gATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCACAGCAGcccca                                                                        <  2:226726/150‑3 (MQ=255)
                                                                                tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatc                                                                       <  2:8751868‑M1/150‑72 (MQ=255)
                                                                                    cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGGCAGGCGTGGCGATAGGGCCGCCAA‑GCGGCCCCGACCATACAAGCGGCCAAGCGGTACGCCAGGGGCGGCGGTGCCGCCCAGGCGCGa                                                                  >  1:1027646/1‑150 (MQ=255)
                                                                                    cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                  >  2:8052571‑M2/1‑144 (MQ=255)
                                                                                         gttgaTGCCGTTCAAAGCCCGATATATCCCCGATATACGCGG‑AATGTGTGGTTGGGTCGTTGCTTCGCGATAGGGCCCCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCaa                                                             <  1:1608485/147‑1 (MQ=255)
                                                                                          gagaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaag                                                             <  1:8346003‑M1/150‑82 (MQ=255)
                                                                                            gttgaGGGTCAACGCCCGCTATACACCGGCTCTACGGGT‑AACGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAC‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggc                                                          >  1:1751967‑M2/3‑136 (MQ=255)
                                                                                            gaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggc                                                           <  1:6001027‑M1/150‑84 (MQ=255)
                                                                                              tGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCGGGGGCTGCTGTGCAGCCCAGtcgcgacgcaaggccg                                                        <  2:4509690‑M2/150‑17 (MQ=255)
                                                                                                 cGGTGAACGCCCGGTCGCGCCCCGCTCTCGGGGG‑ATGGTGTGGGCGGGGCCTGGGGGCGCGTTGGGGCCGCAAC‑GCGGCCCCGACAATTCACGCGGCGAAGCTGAAAGCCTGGGGGTGCTGGGCAGCCGTGTCGCGACGCACGGCCGCTc                                                     <  2:8381461/150‑1 (MQ=255)
                                                                                                  ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCCAGCGGCGAAGCTGTAAGCCGGGGGCGGCTTTGCAGCCCTGtcgcgacgccaggccgcgcc                                                    <  2:5209047‑M2/150‑21 (MQ=255)
                                                                                                      aaaGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccaca                                                <  2:7855435‑M2/150‑25 (MQ=255)
                                                                                                        aGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggcc                                                                            >  1:5246947‑M1/1‑55 (MQ=255)
                                                                                                               ttgagcCCCGCTCTCCGGGG‑ACTGTGTGGGCGCGGCCTGGCGCCGCTATTGGGCCCCAAC‑GCGGCCCCGACAATTCACGGGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGATCCCACAAAGACcgcg                                       <  2:1704375/148‑1 (MQ=255)
                                                                                                                      ccGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTTCGGCGGGATAGGGCCGCCCA‑GCGGCCCCCCCCATGCAAGCGGCGAAGCTGAAATCCGGGGGCTGCGGGGGGGCCCTGGCGCGACGCAAGGCCGCGCCCCCCAAGCCCgcgcccgtca                                >  2:5182499/1‑144 (MQ=255)
                                                                                                                              cgcgG‑AAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAA‑GCCGCCCCGACAAATCAAGCGGCCCAGCTGAAAGCCGGGGGCTGCGGGGCAGCCCAGTCGCGACGCAAGGCCGCGCCCACCAAGACCGCGCAGATCAGTGGATCa                        >  2:3982047/1‑150 (MQ=37)
                                                                                                                                cgG‑AACGTGGGGGTGCGGCCTTGCGTCGCGATCGGGCCGCAAC‑GCGGCCCCGACAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCAcc                      <  2:7935680/150‑1 (MQ=37)
                                                                                                                                      ttggtgggaGCGGCCTTGCGTCGCGActgggctgcccagcagccccaggctttcagcttcgccgctttaattgtcggggccgctgtgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggct                  <  2:2039686‑M1/147‑125 (MQ=255)
                                                                                                                                       gtgtGGGAGCGGGCTTGCGTCGCGCTAGGGCCGCAAA‑GCGGCCCCGACAATTAAAGCGGCGAAGCTGAAAGCCAGGGGCTGCTGTGGCGCCCTGTCGCGACGCAAGGCCGCTCCCCCAAAGACCGCGCACATCAGTGGATCACCCAGCtc                >  2:4888678/1‑150 (MQ=37)
                                                                                                                                        tgtgGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAA‑GCGGCCCCCACAAATCAAGCGGCCGAGCTGAAAGCCTGGGGCTGCTGTGCCGCCCATTAGCCTCCCCAGGCCGCGCCCCCAAAGACCCCGCCGGTCACTGGATCACCCAGCtca               >  2:2491074/1‑150 (MQ=21)
                                                                                                                                          tgGGCGGGGCCTTGGGTCGCTTTCGGGCCGCAAC‑GGGGCCCCGACAATTCACGCGGCCAAGCTGAAAGCCTGGGGCTGCTTTGCAGCCCTTTGGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCtcatc             <  2:1359164/150‑1 (MQ=25)
                                                                                                                                            ggAGCGGCCTTGCGCCGCCAccgggctgcacagcagccccaggcttacagctgcgccgctatatttgtcgggcccgcttggcggccctatcgcgacgcacggccgctcccacactttccgcgtagagcgggtctctagcgggctttcacc            <  2:2567195‑M1/150‑131 (MQ=255)
                                                                                                                                             gtggggCGTTGCGTCGCGAT‑GGGCCCCAAACGCGCCCCCGACAATTCACGCGGCGATGCTTACAGCGTGGGGGTGGTGGGCCGCCCTGCCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAcc          <  2:4634528/146‑1 (MQ=255)
                                                                                                                                             gcgCGGCCTTGCGCCGCTATCGGGCCGAAAC‑GCGGCCCCGACAATTAACGCGGCGACGCTGAACGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGAAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAcc          <  2:4346061/148‑1 (MQ=37)
                                                                                                                                              cgcggcgTTGCGTCGCGTACGGCCCGAACC‑CCGGCCCCGACAATCCCACCGCGGCCGCTGACACCCGGGGCCGGCGGTCCCGCCCGGTCCCCACGCACGCCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAccc         <  2:8284486/149‑1 (MQ=255)
                                                                                                                                               gCGGCCTTGGGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgatttgtcggggccgctttgcggccctatcgcgacgcacggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtc         <  2:2588624‑M1/150‑134 (MQ=255)
                                                                                                                                               gCGGCCTGGCGTCGCGActgggctgcaccgcagcccctggctttcagcttcgccgctttatttgtcggggccgcttggcggccctatcgcgacgcacggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtc         <  2:2588190‑M1/150‑134 (MQ=255)
                                                                                                                                                  gCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAg     >  1:948292/1‑150 (MQ=255)
                                                                                                                                                   ccTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggcccgatcgcgacgcaaggcccctcccacactttccgcggagagcgggtctctagcgggctttcaccggcctca     >  1:946702‑M1/1‑13 (MQ=255)
                                                                                                                                                   ccTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctgtgcggccctatcgcgacgcaaggccgctcccacactgtccgcgtagagcgggtctctagcgggctttcaccgtcatct     >  2:3507025‑M1/1‑13 (MQ=255)
                                                                                                                                                    cTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctc    >  2:2995081‑M1/1‑12 (MQ=255)
                                                                                                                                                    cTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctc    >  1:1399461‑M1/1‑12 (MQ=255)
                                                                                                                                                     ttGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacacattccgcgtatagcgggtctctagcgggctttaaccgtcaacacc   >  1:8363029‑M1/1‑11 (MQ=255)
                                                                                                                                                     ttGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc  >  1:8481763/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGG‑AAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAA‑GCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCC  >  NC_002947/667926‑668223

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: ‑