Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,074 T→G 17.5% 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,0740TG17.5% 409.9 / 13.3 179intergenic (+35/+117)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (53/89);  new base G (34/3);  total (87/92)
Fisher's exact test for biased strand distribution p-value = 1.12e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

GCAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCC  >  NC_002947/667926‑668223
                                                                                                                                                    |                                                                                                                                                     
gCAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTT                                                                                                                                                      <  1:1681286‑M1/150‑1 (MQ=255)
gCAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTT                                                                                                                                                      <  1:7639359‑M1/150‑1 (MQ=255)
 cAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTG                                                                                                                                                     >  2:7883735‑M1/1‑150 (MQ=255)
 cAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTG                                                                                                                                                     >  1:7304277‑M1/1‑150 (MQ=255)
  aGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGC                                                                                                                                                    <  1:579465‑M1/150‑1 (MQ=255)
  aGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGC                                                                                                                                                    >  2:8558889‑M1/1‑150 (MQ=255)
   gTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCG                                                                                                                                                   <  2:5577796‑M1/150‑1 (MQ=255)
   gTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCG                                                                                                                                                   <  1:6012844‑M1/150‑1 (MQ=255)
      aCCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCG                                                                                                                                                >  2:8791434‑M1/1‑150 (MQ=255)
       ccccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGC                                                                                                                                               >  2:732446‑M1/1‑150 (MQ=255)
       ccccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGC                                                                                                                                               >  2:4739293‑M1/1‑150 (MQ=255)
        cccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTGGCGGCGCG                                                                                                                                              >  1:8034679‑M1/1‑150 (MQ=255)
         ccGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGA                                                                                                                                             <  1:4744418‑M1/150‑1 (MQ=255)
                ccTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGcc                                                                                                                                      <  1:1202639/150‑1 (MQ=255)
                  tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCCAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGc                                                                                                                                    >  2:8268493/1‑150 (MQ=255)
                  tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGTCGCGATAGGGgcgc                                                                                                                                    >  1:5070485/1‑150 (MQ=255)
                  tACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGc                                                                                                                                    >  2:6328302/1‑150 (MQ=255)
                   aCCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCGAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCa                                                                                                                                   <  2:581751/150‑1 (MQ=255)
                     cgcgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGTCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaaa                                                                                                                                 <  2:5437179/150‑1 (MQ=255)
                     cgcgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaaa                                                                                                                                 >  2:3617274/1‑150 (MQ=255)
                       cgTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGc                                                                                                                               <  2:1660648/150‑1 (MQ=255)
                        gTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGGGCGGTAGGGCCGCAAAGCg                                                                                                                              >  1:5226861/1‑150 (MQ=255)
                         tGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCgg                                                                                                                             <  2:7254911/150‑1 (MQ=255)
                          gccgccgccGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATATGGCCGCAAAGCGGc                                                                                                                            <  2:7751839/150‑1 (MQ=255)
                            cgccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGACAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGccc                                                                                                                          <  2:8426199/150‑1 (MQ=255)
                             gccgccgcAGGTCGAATCGATGTACCGGCCGAGCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcccc                                                                                                                         <  2:3461046/150‑1 (MQ=255)
                             gccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcccc                                                                                                                         <  1:6555228/150‑1 (MQ=255)
                             gccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcccc                                                                                                                         <  1:4644013/150‑1 (MQ=255)
                             gccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcccc                                                                                                                         <  2:8870396/150‑1 (MQ=255)
                               cgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGa                                                                                                                       <  1:2949786/150‑1 (MQ=255)
                                       gTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                            >  1:242657/1‑137 (MQ=255)
                                         cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGTTGGGCAAGGGGAGTTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgc                                                                                                             <  1:9016602/150‑1 (MQ=255)
                                         cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgc                                                                                                             <  2:4308065/150‑1 (MQ=255)
                                            aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcggc                                                                                                          <  2:1772045/150‑1 (MQ=255)
                                            aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGAAAGGGGAGATGACGGTGAAAGCCCGCTAGACACCCCCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcggc                                                                                                          <  2:2398949/150‑1 (MQ=38)
                                            aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGGGGCGGTGGGGCGGGCAAGCGGCCCCGCCCAGTCAAtcggg                                                                                                          >  1:6936301/1‑149 (MQ=255)
                                               gATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGaa                                                                                                       <  2:5933028/150‑1 (MQ=255)
                                                aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCCTTAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAg                                                                                                      >  2:1179929/1‑150 (MQ=255)
                                                aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAg                                                                                                      >  2:4035483/1‑150 (MQ=255)
                                                    aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGa                                                                                                  <  1:5686648/150‑1 (MQ=255)
                                                    aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACACATCAAGCGGCGAAGCTGa                                                                                                  >  2:8807668/1‑150 (MQ=255)
                                                     ccGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCAAAGCGGCCCCGACAATTAAAGCGGCGACGCGGaa                                                                                                 >  1:8244246/1‑150 (MQ=39)
                                                      cGGCCGAACCTGCTGGCGGTGGGCGTGGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGGGTAATGTGTGGGCGCGGCCTTGCGTCGCGATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                <  1:4833021/150‑1 (MQ=37)
                                                      cGGCCGAACCTGCAGGCGGTGGGCGATGGGAATGGGGTGATTACGGTTAAAGCCCGCTCGAGACCCGCTCTACGCGGAAAGTGTGGGAGGGGCCTTGGGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                <  1:3151362/150‑1 (MQ=255)
                                                      cGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                <  2:3911750/150‑1 (MQ=255)
                                                      cGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCCAAGCGGCCCCCACAATTCAAGCGGCGACGCTGgaa                                                                                                >  2:6912481/1‑150 (MQ=255)
                                                           gAACCTGCAGGCCGTGGGCGATGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAACGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTg                                                                                           <  1:547065/150‑1 (MQ=255)
                                                           gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGTCGCGATAGGGCCGCCCGGCGGCCCCGCCAAATCAAGCGGCGAACCTGAAAGccag                                                                                           >  2:1006307/1‑148 (MQ=255)
                                                           gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCTGAAAGCCTg                                                                                           >  2:6539462/1‑150 (MQ=255)
                                                            aaCCTGCTGGCCGTGGGCGATGGGCACGGGGAGTTGACGGTGAACCCCCGCTCGCGACCCGCTCTACGCGGAATGTGTGGGAGGGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                          <  1:3301644/150‑1 (MQ=255)
                                                            aaCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgg                                                                                          <  2:3246878/150‑1 (MQ=255)
                                                             aCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGTCGGTGAAAGCCCGCTATAGACCCTCTCTGCGCGGTCTTTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagaagccccatgctttctgcttagccgattgaattgtcggt                                                                                         <  2:5099012‑M1/150‑53 (MQ=255)
                                                                 gctggCCGTGGGCTTTGGGCATGGGGATATGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGGAACGTGTGGGAGCGGCCTTGGGTCGCGATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctg                                                                                     <  1:4179694/147‑1 (MQ=38)
                                                                     gcggggggCGTTGGGCATGGGTAGTTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAATGTGTGGGAGCGGCCTTGCGTCGCTATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgt                                                                                 <  1:1104616/145‑1 (MQ=37)
                                                                     gcggggggCGTTGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAACGGGTGGGCGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgt                                                                                 <  1:1104555/145‑1 (MQ=255)
                                                                      ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                <  2:7129599/150‑1 (MQ=255)
                                                                      ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                <  2:1369845/150‑1 (MQ=255)
                                                                      ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                <  2:8211475/150‑1 (MQ=255)
                                                                       cGTGGGAGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCAACGCGGAAAGTGTGCGAGAGGCCTATCGACGAGATAGGGCCGCAAAGAGCCCACGACAATTCAAGCGGCCAAACTGAAATCATGGGGCTGCagttc                                                                               >  2:6807064/1‑148 (MQ=255)
                                                                        gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                              <  2:2128318/150‑1 (MQ=255)
                                                                        gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                              <  1:4511945/150‑1 (MQ=255)
                                                                        gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                              <  2:158963/150‑1 (MQ=255)
                                                                          gggCGTGGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGCGCGGCCTTGCGAGGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                            <  1:893730/150‑1 (MQ=38)
                                                                             cGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCTATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                         <  1:1917380/150‑1 (MQ=255)
                                                                             cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCCTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                         >  2:2449896/1‑150 (MQ=255)
                                                                               ttgGGCATGGGGTTATGACGGTGAAAGCCCGCTAGAGACCCGCTCTCCGCGGAACGTGTGGGAGCGGCGTTGCGTGGCTATCGGGCCCCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGt                                                                       >  1:3871822‑M2/2‑149 (MQ=255)
                                                                               aTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGACAGGGTGGGAGCGGCCTTGCGTCGCGTTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCACAGCAGccccag                                                                       <  2:1915952/150‑4 (MQ=38)
                                                                                tGGGCATGGGTATATGACGTTTAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCATTGCGTCGCGATAGGCCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtc                                                                      >  1:7883735‑M2/1‑148 (MQ=255)
                                                                                tGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtc                                                                      >  2:2643254‑M2/1‑148 (MQ=255)
                                                                                    cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcga                                                                  <  1:2398213‑M1/150‑76 (MQ=255)
                                                                                    cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcga                                                                  <  2:5224943‑M1/150‑76 (MQ=255)
                                                                                    cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATTGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                  >  1:4602684‑M2/1‑144 (MQ=255)
                                                                                         ggAGATGACGGTGAACGCCCGCTCGAGCCCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaa                                                             >  1:7270976‑M2/1‑139 (MQ=255)
                                                                                         ggAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaa                                                             >  2:1870458‑M2/1‑139 (MQ=255)
                                                                                           agaTGCCGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGACAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACACTTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                           >  2:666338‑M2/1‑137 (MQ=255)
                                                                                           agaTGACGGTGAAAGCCCGCTAGAGCCCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAAGTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                           >  2:5242353‑M2/1‑137 (MQ=255)
                                                                                           agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                           >  1:5540173‑M2/1‑137 (MQ=255)
                                                                                           agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                           >  1:5540235‑M2/1‑137 (MQ=255)
                                                                                                         gCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCaa                                                                                                               <  2:922117/84‑1 (MQ=255)
                                                                                                           ccGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGGTAGGGCCCCAAAGCGGCCCCGACAAATAAAGCGGCGCAGCGGAAAGCCTGGGGCGGCTGTGCAGCCCAGTCGCGCCGCCACGCCGCGCCCCCCAAGAc                                           >  2:8410683/1‑150 (MQ=255)
                                                                                                              cTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGGGTGGCGATTGGGCCGCGAAGCGGGCCCGACCAATGAAGGGGGGGAGGTGTAAGGCTGGGGGTGGTGTTGAGGCCTGTGGGGTCGGAAGGGCGGTGCCCCAAAGaaacc                                        >  2:4031209/1‑146 (MQ=255)
                                                                                                                    ccccGCTCTCCGGGTAATGTGTGGGTGGGGCGTGGGGTCGCTATCGGCCCGCACCGCGGCCCCGACAATACACGCGGCGACGCTCAATGCCGGGGGCTGCTTTGCCGCCCTGGGGGGACGAACGGCCGCTCCCACAAAGACCGCGCAGAt                                  <  2:1236501/149‑1 (MQ=255)
                                                                                                                      ccGCTCTACGCGGAACGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGCCAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCATGGCCGCTCCCACAAAGACCGCGCAGATCa                                <  2:7477392/150‑1 (MQ=37)
                                                                                                                      ccGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCCCCAATCGGCCCCCAACAATAAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCAGTCGCGACGCAAGGCCCCTACCACAAAGACCGCGCAGAtaa                                >  2:2930525/1‑148 (MQ=25)
                                                                                                                          tctACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                         <  2:497841/105‑1 (MQ=255)
                                                                                                                                 ggAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCCCGCAGACCAGTGGATCAccc                     >  2:9018184/1‑150 (MQ=255)
                                                                                                                                           ggAGTGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcacc           <  2:4921382‑M1/150‑131 (MQ=255)
                                                                                                                                           ggAGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctgggcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggcttgcacc           <  2:2896398‑M1/150‑131 (MQ=255)
                                                                                                                                            gcgcgGCCGGGCGTCGCGAcagggcggcacagcagccccaggctttcagcttcgccgcttgatttgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccg          <  2:8724223‑M1/148‑132 (MQ=255)
                                                                                                                                             aGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggcccgatagcgacgcaaggccgcgcccacacattccgcgtagatcgggtctctagcgggcgttcaccgt         >  1:2137990‑M1/1‑18 (MQ=255)
                                                                                                                                             aGCGGCCTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgcgttgcggccctatagccacccaaggccgcgcccacacattacgcctagagcgggtctcaagcgggctttcaccca         >  1:2139676‑M1/1‑18 (MQ=255)
                                                                                                                                               cGGCCTTGCGTCGCGTTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCCAAGCTGAAATCCTGGGGCTGCTGTGAAGCCCTGTCGCGACGCAATGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCcac       >  1:8283449/1‑150 (MQ=37)
                                                                                                                                                ggCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCcaca      >  1:8537128/1‑150 (MQ=255)
                                                                                                                                                  ccTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgacttgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcggggctctagcgggctttcaccgtcatct    >  2:5242353‑M1/1‑13 (MQ=255)
                                                                                                                                                  ccTTGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatct    >  2:8678570‑M1/1‑13 (MQ=255)
                                                                                                                                                  ccTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGc    >  2:2991710/1‑150 (MQ=255)
                                                                                                                                                    ttGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcacggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcc  <  2:614634‑M1/150‑140 (MQ=255)
                                                                                                                                                    ttGCGTCGCGAcagggctgcacagcagccccaggctttcagcttcgccgcttgaattgtcggggccgctttgcggccctatcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctcc  >  1:6804730‑M1/1‑11 (MQ=255)
                                                                                                                                                    ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc  >  1:7532102/1‑150 (MQ=255)
                                                                                                                                                    ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc  >  2:6848268/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGTGACCCCGGGCACCTACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCC  >  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: