Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,020,956 A→G 32.8% intergenic (‑37/‑182) sulA ← / → sxy SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,020,9560AG32.8% 78.2 / 58.7 61intergenic (‑37/‑182)sulA/sxySOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base A (23/18);  new base G (11/9);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GTAGCCCAGATGCCTGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGA  >  NC_000913/1020670‑1021237
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                         
gTAGCCCAGATGCCTGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtac                                                                                                                                                                                                                                                                                          <  1:781833/288‑1 (MQ=255)
    cccAGATGCCTGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagtac                                                                                                                                                                                                                                                                                     >  2:90684/1‑289 (MQ=255)
       aGATGCCTGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACAt                                                                                                                                                                                                                                                                                   >  2:541500/1‑288 (MQ=255)
         aTGCCTGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCa                                                                                                                                                                                                                                                                                >  2:290720/1‑289 (MQ=255)
           gCCTGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGt                                                                                                                                                                                                                                                                              <  1:466663/289‑1 (MQ=255)
             cTGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGa                                                                                                                                                                                                                                                                            <  2:720875/289‑1 (MQ=255)
                       tCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGAt                                                                                                                                                                                                                                                                  <  1:90684/289‑1 (MQ=255)
                        cccGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATc                                                                                                                                                                                                                                                                 <  2:165375/289‑1 (MQ=255)
                                            ggTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAaga                                                                                                                                                                                                                                              >  1:860850/1‑288 (MQ=255)
                                                ttAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTTAGTTACTGTATGGATGTACAGTACATCCAGTGACCACAACGATCAACCCTAATTTCGGCAAGAGCCt                                                                                                                                                                                                                                          >  1:63966/1‑288 (MQ=255)
                                                     cAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGc                                                                                                                                                                                                                                       >  2:724971/1‑286 (MQ=255)
                                                     cAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGc                                                                                                                                                                                                                                       <  1:724971/286‑1 (MQ=255)
                                                                                   tgtgggaaCAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCGGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAt                                                                                                                                                                                                      <  1:4927/283‑1 (MQ=255)
                                                                                           cAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTg                                                                                                                                                                                              <  2:348781/289‑1 (MQ=255)
                                                                                            aaTGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc                                                                                                                                                                                             <  2:439638/289‑1 (MQ=255)
                                                                                               ggCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGcc                                                                                                                                                                                            >  2:752416/1‑287 (MQ=255)
                                                                                                   aCAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa                                                                                                                                                                                      >  2:789361/1‑289 (MQ=255)
                                                                                                               gCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTg                                                                                                                                                                            >  2:529567/1‑287 (MQ=255)
                                                                                                               gCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTg                                                                                                                                                                            >  1:199283/1‑287 (MQ=255)
                                                                                                                      catGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTc                                                                                                                                                                   <  1:261812/289‑1 (MQ=255)
                                                                                                                       atGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATTAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCg                                                                                                                                                                  <  2:199283/289‑1 (MQ=255)
                                                                                                                       atGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCg                                                                                                                                                                  <  2:429111/289‑1 (MQ=255)
                                                                                                                               aTCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGcc                                                                                                                                                                          <  2:138657/273‑1 (MQ=255)
                                                                                                                               aTCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGcc                                                                                                                                                                          >  1:138657/1‑273 (MQ=255)
                                                                                                                               aTCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAAc                                                                                                                                                          >  1:641752/1‑289 (MQ=255)
                                                                                                                               aTCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAAc                                                                                                                                                          >  2:387085/1‑289 (MQ=255)
                                                                                                                                       gATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGt                                                                                                                                                  <  2:268066/289‑1 (MQ=255)
                                                                                                                                       gATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCg                                                                                                                                                   <  1:541500/288‑1 (MQ=255)
                                                                                                                                                  tCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATc                                                                                                                                       <  1:290720/289‑1 (MQ=255)
                                                                                                                                                   cACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCa                                                                                                                                      >  2:400290/1‑289 (MQ=255)
                                                                                                                                                   cACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCa                                                                                                                                      <  2:388191/289‑1 (MQ=255)
                                                                                                                                                   cACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCa                                                                                                                                      <  1:752554/289‑1 (MQ=255)
                                                                                                                                                        aTAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagta                                                                                                                                                                                                                                                                                      <  1:762977/140‑1 (MQ=255)
                                                                                                                                                        aTAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagta                                                                                                                                                                                                                                                                                      >  2:762977/1‑140 (MQ=255)
                                                                                                                                                          aaGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTTGTAGTTAACGGATCCGTTTATGTGAAATATTCtttt                                                                                                                               >  1:194540/1‑289 (MQ=255)
                                                                                                                                                                                gTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAAt                                                                                                         <  2:194540/289‑1 (MQ=255)
                                                                                                                                                                                                          tgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCt                                                                                                                                  <  2:91218/238‑1 (MQ=255)
                                                                                                                                                                                                          tgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCt                                                                                                                                  >  1:91218/1‑238 (MQ=255)
                                                                                                                                                                                                                                   tgtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTg                                                      >  2:83822/1‑289 (MQ=255)
                                                                                                                                                                                                                                     tgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc                                                                                                                                                                                             <  2:603022/152‑1 (MQ=255)
                                                                                                                                                                                                                                     tgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc                                                                                                                                                                                             >  1:603022/1‑152 (MQ=255)
                                                                                                                                                                                                                                        aTAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAAc                                                 >  2:70630/1‑289 (MQ=255)
                                                                                                                                                                                                                                         tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACg                                                <  2:641752/289‑1 (MQ=255)
                                                                                                                                                                                                                                            ccTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTg                                             >  1:668250/1‑289 (MQ=255)
                                                                                                                                                                                                                                               gAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCtata                                                                       >  2:380333/1‑260 (MQ=255)
                                                                                                                                                                                                                                               gAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCtata                                                                       <  1:380333/260‑1 (MQ=255)
                                                                                                                                                                                                                                               gAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGc                                          >  2:423837/1‑289 (MQ=255)
                                                                                                                                                                                                                                               gAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGc                                          >  1:582246/1‑289 (MQ=255)
                                                                                                                                                                                                                                                aaGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGc                                          >  2:603135/1‑288 (MQ=255)
                                                                                                                                                                                                                                                 aGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGca                                         <  1:400290/288‑1 (MQ=255)
                                                                                                                                                                                                                                                 aGTATACAACATCAATCCGACCCCAACGAGTTACTGTATGGGTGTACAGTACATCCAGTGACAACAAAGCTAAACCCTATTTTCCGAAAGAACCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCATTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATCTTATGATTTTAAAAGGAATTTTATGAAAAGACTCTCCTATAAACGTATCTATAAATCACAAGAATACCTGGCAACGTTGGGca                                         <  1:266212/288‑1 (MQ=255)
                                                                                                                                                                                                                                                      aCATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTaa                                                                                                                                                           <  2:803911/169‑1 (MQ=255)
                                                                                                                                                                                                                                                      aCATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTaa                                                                                                                                                           >  1:803911/1‑169 (MQ=255)
                                                                                                                                                                                                                                                               atcCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGa                           >  1:774706/1‑288 (MQ=255)
                                                                                                                                                                                                                                                               atcCAGCCCCTGTGAGTTACTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGa                           >  1:238920/1‑288 (MQ=255)
                                                                                                                                                                                                                                                               atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTaa                                                                                                                                                           >  1:157791/1‑160 (MQ=255)
                                                                                                                                                                                                                                                               atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTaa                                                                                                                                                           <  2:157791/160‑1 (MQ=255)
                                                                                                                                                                                                                                                               atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGa                           >  2:793402/1‑288 (MQ=255)
                                                                                                                                                                                                                                                               atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGa                           >  2:286397/1‑288 (MQ=255)
                                                                                                                                                                                                                                                               atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAAAGTTGGGCACAATTGAATACCGa                           >  1:546561/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                         tgtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTgg                <  1:603135/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCGTTAAGAATAAGTTGCCTATGTTCGTAGTTAACGGATCTGTTAATGTGAATCATTCTTTTATGTTCTGATTTTGAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATGTATAAATCACAAGAATACTTGGCAACGTTGGGCACAATTGAATACGGATCATTGTTTGGc               >  2:574742/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                          gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGc               >  2:60196/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                 aCTGTATGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTCCa        >  2:767520/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                       tGGATGTGCAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTg   >  2:12056/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                             tacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGa  <  2:743340/283‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             tacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATACGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGa  >  1:743340/1‑283 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                         
GTAGCCCAGATGCCTGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGA  >  NC_000913/1020670‑1021237

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