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