Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,543,571 A→T 6.0% intergenic (+64/+57) PP_4871 → / ← PP_4872 lysine decarboxylase family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,543,5710AT 6.0% 688.0 / 17.9 245intergenic (+64/+57)PP_4871/PP_4872lysine decarboxylase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (121/101);  new base T (7/10);  total (133/112)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

GCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTACGCG  >  NC_002947/5543423‑5543719
                                                                                                                                                    |                                                                                                                                                    
gcagcaGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAt                                                                                                                                                     <  2:7713311/150‑1 (MQ=255)
gcagcaGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAt                                                                                                                                                     <  1:5016356/150‑1 (MQ=255)
  agcagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgg                                                                                                                                                   <  1:6806923/150‑1 (MQ=255)
    cagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                                 <  2:3531735/150‑1 (MQ=255)
    cagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                                 <  2:1148153/150‑1 (MQ=255)
    cagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                                 <  1:543208/150‑1 (MQ=255)
      gcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGa                                                                                                                                               <  1:4400358/150‑1 (MQ=255)
             gAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCt                                                                                                                                        <  1:123907/150‑1 (MQ=255)
             gAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCt                                                                                                                                        >  2:6003535/1‑150 (MQ=255)
             gAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCt                                                                                                                                        <  1:3924191/150‑1 (MQ=255)
              cgctgctGGAAGGGATGGACAGCTGTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                                                                       <  1:5784220/149‑1 (MQ=255)
              aGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCtt                                                                                                                                       <  1:7586848/150‑1 (MQ=38)
                ctgctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTtgt                                                                                                                                     <  2:3015724/150‑1 (MQ=38)
                ctgctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTtgt                                                                                                                                     <  1:236671/150‑1 (MQ=255)
                     ggAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcgc                                                                                                                                <  1:8065337/150‑1 (MQ=255)
                     ggAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcgc                                                                                                                                <  1:2484140/150‑1 (MQ=255)
                        aGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACACGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACACAGGCAATGGGGGAGCGGCCTTGTGTCGCGAt                                                                                                                             <  1:7226512/150‑1 (MQ=255)
                           gATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATggg                                                                                                                          <  2:5940257/150‑1 (MQ=255)
                               gACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGc                                                                                                                      >  2:3700778/1‑150 (MQ=255)
                                aCAGGTTTGTGGCGCCGGTCGCGCCTAGGTGGGGCGACACGCAGCCTGAGTGAAGGTGAGAATTTTGGGACCGCTTCGGGGTCCTTTCGCGACACAAGGCGGCTCCCAAAAAGGAATTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCa                                                                                                                     <  1:5381876/150‑1 (MQ=255)
                                aCAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCa                                                                                                                     >  2:5972285/1‑150 (MQ=255)
                                 cAGCTTGGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCaa                                                                                                                    <  2:3439322/150‑1 (MQ=255)
                                  aGCTTTGTGGCGCGGGGCGCGCCTAAGTGGGTCGACACGCAGCCTGACTGATGCGGTGACTTGTGGGACCGCTTCGGGGACCTTTCGCGACACAAGGCCGCTCCCACAATGGCATTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCaaa                                                                                                                   <  1:4785616/150‑1 (MQ=37)
                                    cTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCCAGGGGGGAGCGGCCGTGTGGCGCGATGGGCTGCAAagc                                                                                                                 >  1:1505396/1‑150 (MQ=255)
                                    cTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTGGTGTCGCGATGGGCTGCAAagc                                                                                                                 >  2:2782710/1‑150 (MQ=255)
                                         tgGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGccc                                                                                                            >  2:7455667/1‑150 (MQ=255)
                                           gcgcCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGGGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCa                                                                                                          >  2:3352965/1‑150 (MQ=255)
                                           gcgcCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCa                                                                                                          >  1:6341633/1‑150 (MQ=255)
                                           gcgcCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCAAAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCa                                                                                                          >  1:3630762/1‑150 (MQ=255)
                                              ccGGTCGCGCCTACGGGGGCCTACAAGCAGCCTGACTGACGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTGCGCGACACAAGGCCGCTCCCAAAACGGCAAGGGGGGAGCGGCCTGGTGCCGCGAGGGGCTGCTTTGCAGcccatcgc                                                                                                       <  1:8804701/150‑6 (MQ=255)
                                                      gcCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAATGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCa                                                                                               <  2:79422/150‑1 (MQ=255)
                                                      gcCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCCATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCa                                                                                               >  1:8185642/1‑150 (MQ=255)
                                                       ccTATGGGGGTCGACAAGCAGCCTGACTGATGCTTAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCTACACACGGCCGCCCCCCCAAAGGAAAGGGGGGCGGGGCCTTGGGCCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAg                                                                                              <  1:3793827/150‑1 (MQ=255)
                                                                     cAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCACTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGACTGa                                                                                >  2:4980169/1‑150 (MQ=255)
                                                                     cAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCCGTTAACGCGGAATGa                                                                                >  2:3225588/1‑150 (MQ=255)
                                                                     cAAGCAGCCTGACTGAAGCTGAGACTGTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTAGATCAGTTCACGCGGAATgc                                                                                >  2:6824573/1‑149 (MQ=255)
                                                                      aaGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCCATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCCAAAGTCGATCAGTTAACGCGGAATGAc                                                                               >  1:3013194/1‑150 (MQ=255)
                                                                      aaGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGAc                                                                               >  1:6326044/1‑150 (MQ=255)
                                                                      aaGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAAAGCGGAATGAc                                                                               >  1:7456471/1‑150 (MQ=255)
                                                                        gcagcCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAg                                                                             >  2:5182422/1‑150 (MQ=255)
                                                                             ctgactgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACACAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCCCTAGTCGATCCGTTCACGCGGAATGACAGGCTGa                                                                        >  2:851562/1‑150 (MQ=39)
                                                                                    aaGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCCATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGACTGACAGGCTGACGCGGct                                                                 >  1:7973702/1‑150 (MQ=255)
                                                                                        tGAGACTTTTGGGCCCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttctt                                                             <  1:8264688/150‑1 (MQ=255)
                                                                                        tGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttctt                                                             >  1:5558830/1‑150 (MQ=255)
                                                                                         gagaCTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTGAACGCGGAATGACAGGCTGACGCGGCTTCttg                                                            >  2:2334548/1‑150 (MQ=255)
                                                                                         gagaCTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATCGGCAGCAAAGCAGCCCCAATAGTCGAGCAGTGAACGCGGAATGACAGGCAGACGCGGCTTCttg                                                            >  2:5625335/1‑150 (MQ=39)
                                                                                         gagaCGTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCAAGAGTCGAGCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttg                                                            >  2:5625195/1‑150 (MQ=255)
                                                                                             cTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGc                                                        <  2:7524023/150‑1 (MQ=255)
                                                                                              ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCa                                                                                               >  1:2068845/1‑110 (MQ=255)
                                                                                              ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCa                                                                                               <  2:2068845/110‑1 (MQ=255)
                                                                                              ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTTcc                                                       >  1:4328127/1‑150 (MQ=255)
                                                                                              ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGcc                                                       >  2:5516304/1‑150 (MQ=255)
                                                                                              ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAAAGCGGAATGACAGGCTGAAGCGGCTTATTGTTGcc                                                       >  1:7845686/1‑150 (MQ=255)
                                                                                              ttttGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCACCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAAACCGGAATGAAAGGCTGACGCGGCTTCTTGTTGcc                                                       >  1:4824899/1‑150 (MQ=255)
                                                                                                aaaggACCGCGAAGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCGTCTTGTTGCCCt                                                     >  1:6036270/4‑150 (MQ=39)
                                                                                                 tGGGACCGCTTCGCGGGCCTTTCGCGACCCACGGCCGCTCCCACAACGGCAATGGGGGAGCGGCCTGGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCtt                                                    <  1:8344022/150‑1 (MQ=39)
                                                                                                  gggACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCATTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTg                                                   <  1:6592775/150‑1 (MQ=255)
                                                                                                   ggACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGc                                                  <  2:6116505/150‑1 (MQ=255)
                                                                                                          ttCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCtt                                           >  2:7820548/1‑150 (MQ=255)
                                                                                                            cgcgGCCCATCCGCGCCCCACGGCCGCCCCCACAACGGCACTGGGGGAGCGGCCTTGTGACGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGc                                         <  1:6122374/150‑1 (MQ=38)
                                                                                                              cgGTCCTTTCGCGACACATGGCCGCTCCCAAAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTTACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTc                                       <  2:117344/150‑1 (MQ=255)
                                                                                                                gTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGc                                     <  2:1826204/150‑1 (MQ=255)
                                                                                                                  ccTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCct                                   <  2:5337680/150‑1 (MQ=255)
                                                                                                                     ttCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCGTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCtt                                >  2:6589888/1‑150 (MQ=255)
                                                                                                                      tCGCGACACCAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCttt                               <  2:4828628/150‑1 (MQ=255)
                                                                                                                      tCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGACAGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCttt                               >  2:5902129/1‑150 (MQ=25)
                                                                                                                      tCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCttt                               <  2:1901175/150‑1 (MQ=255)
                                                                                                                        gcgACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgc                             >  1:5888335/1‑150 (MQ=255)
                                                                                                                        gcgACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgc                             >  2:3931258/1‑150 (MQ=255)
                                                                                                                        gcgACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgc                             >  2:843524/1‑150 (MQ=255)
                                                                                                                        gcgACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCGTCTTGTTGCCCGTGCCACCCTTGCTCGCCTCTTTgc                             >  1:4529992/1‑150 (MQ=255)
                                                                                                                          gACACAATGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCGTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgcgc                           >  2:7967473/1‑150 (MQ=255)
                                                                                                                           acacAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGAGCCGTTAACGCGGAATGACAGGCTGACGCGGCTTCGTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCt                          >  2:5343844/1‑150 (MQ=255)
                                                                                                                            cacaAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTc                         <  2:1964148/150‑1 (MQ=255)
                                                                                                                               aaGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCtt                      >  2:3033880/1‑150 (MQ=255)
                                                                                                                               aaGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCtt                      >  2:7018645/1‑150 (MQ=255)
                                                                                                                                aGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTg                     >  1:7865989/1‑150 (MQ=38)
                                                                                                                                aGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTg                     >  1:1715858/1‑150 (MQ=255)
                                                                                                                                aGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTg                     >  1:4203970/1‑150 (MQ=255)
                                                                                                                                 ggCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTgg                    >  1:2250291/1‑150 (MQ=38)
                                                                                                                                    cGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAg                 >  2:4082249/1‑150 (MQ=255)
                                                                                                                                    cGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAg                 >  2:8419964/1‑150 (MQ=255)
                                                                                                                                    cGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTGGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAg                 >  2:250159/1‑150 (MQ=255)
                                                                                                                                     gCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGc                >  1:4803469/1‑150 (MQ=255)
                                                                                                                                       tCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCct              >  2:404270/1‑150 (MQ=255)
                                                                                                                                       tCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCct              >  1:6911669/1‑150 (MQ=255)
                                                                                                                                         ccACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgc            >  1:1408097/1‑150 (MQ=255)
                                                                                                                                         ccACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgc            >  2:4394712/1‑150 (MQ=255)
                                                                                                                                          cacaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttg                                                            >  2:4615248/1‑101 (MQ=255)
                                                                                                                                          cacaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCcggcg           >  1:7054891/1‑149 (MQ=255)
                                                                                                                                           acaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgctg          <  2:6657651/150‑1 (MQ=255)
                                                                                                                                           acaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgctg          >  2:4226359/1‑150 (MQ=255)
                                                                                                                                           acaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgctg          >  1:3894017/1‑150 (MQ=255)
                                                                                                                                           acaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgctg          <  2:861715/150‑1 (MQ=255)
                                                                                                                                               aGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTa      <  2:734967/150‑1 (MQ=255)
                                                                                                                                               aGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTa      >  2:8083656/1‑150 (MQ=255)
                                                                                                                                                ggCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttc                                                               >  1:7363840/1‑92 (MQ=255)
                                                                                                                                                ggCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTAc     >  2:6235489/1‑150 (MQ=255)
                                                                                                                                                   aaTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTTGCAGCCTGCTGGTTAcgcg  >  1:3211997/1‑150 (MQ=255)
                                                                                                                                                   aaTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTAcgcg  >  1:2822356/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGTTACGCG  >  NC_002947/5543423‑5543719

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