Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 418,428 Δ1 bp 50.0% coding (540/1296 nt) phoR → sensory histidine kinase in two‑component regulatory system with PhoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913418,4280G.50.0% 72.7 / 82.0 42coding (540/1296 nt)phoRsensory histidine kinase in two‑component regulatory system with PhoB
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base . (10/11);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAAAACGCCGCCGTGAACTGGGCAATCTGATTAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCT  >  NC_000913/418145‑418709
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                         
aTAAAAAACGCCGCCGTGAACTGGGCAATCTGATTAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTta                                                                                                                                                                                                                                                                                      <  1:61320/289‑1 (MQ=255)
                       ggCAATCTGATTAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGAtgg                                                                                                                                                                                                                                                               >  1:270614/1‑289 (MQ=255)
                            tCTGATTAAACGCTTTCGTAGCGGCGCGTAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTGGCGCGCAACCAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTCCGTTACCCGGAGTTTACGCAATATCTGCAAATGTGTGATTTTTCTCGCCCGCCCAATCTGGTTCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACACACAAAACTTTGCTGATGGTGgc                                                                                                                                                                                                                                                           >  2:173766/1‑289 (MQ=255)
                                  tAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATg                                                                                                                                                                                                                                                   <  1:20505/289‑1 (MQ=255)
                                          ttCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCa                                                                                                                                                                                                                                            >  1:305813/1‑289 (MQ=255)
                                                                                   aCCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCa                                                                                                                                                                                                   >  2:235915/1‑289 (MQ=255)
                                                                                               ggCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATg                                                                                                                                                                                       >  1:99069/1‑289 (MQ=255)
                                                                                               ggCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAGACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGCGCCATGa                                                                                                                                                                                      >  1:318333/1‑289 (MQ=255)
                                                                                                                tAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCa                                                                                                                                                                                                                                                <  2:278710/214‑1 (MQ=255)
                                                                                                                tAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCa                                                                                                                                                                                                                                                >  1:278710/1‑214 (MQ=255)
                                                                                                                 aaCGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGa                                                                                                                                                                     <  2:270614/289‑1 (MQ=255)
                                                                                                                 aaCGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGAc                                                                                                                                                                    <  1:110587/289‑1 (MQ=255)
                                                                                                                   cGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGAcc                                                                                                                                                                   <  1:235915/289‑1 (MQ=255)
                                                                                                                       cTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGgc                                                                                                                                                                                                                       <  2:67510/233‑1 (MQ=255)
                                                                                                                       cTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGgc                                                                                                                                                                                                                       >  1:67510/1‑233 (MQ=255)
                                                                                                                        tG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTta                                                                                                                                                                                                                                                                                      >  2:213592/1‑169 (MQ=255)
                                                                                                                        tG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTta                                                                                                                                                                                                                                                                                      <  1:213592/169‑1 (MQ=255)
                                                                                                                        tG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTtata                                                                                                                                                                                                                                                                                    <  1:110777/170‑1 (MQ=255)
                                                                                                                        tG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTtata                                                                                                                                                                                                                                                                                    >  2:110777/1‑170 (MQ=255)
                                                                                                                                aacaaATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAggg                                                                                                                                                        <  2:58311/288‑1 (MQ=255)
                                                                                                                                 acaaATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGtt                                                                                                                                                      >  2:120496/1‑289 (MQ=255)
                                                                                                                                       tCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTg                                                                                                                                                 <  2:99069/289‑1 (MQ=255)
                                                                                                                                                  gcTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGatgaat                                                                                                                                     <  1:120496/289‑1 (MQ=255)
                                                                                                                                                                aaCGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAACATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGa                                                                                                                       >  1:303360/1‑289 (MQ=255)
                                                                                                                                                                     gCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGACGCTGGAAGgc                                                                                                                   >  1:39776/1‑288 (MQ=255)
                                                                                                                                                                          aCATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGt                                                                                                              >  1:101059/1‑289 (MQ=255)
                                                                                                                                                                                                  cgaggTTTACCCACTATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAACCTGGTGCTCAACCCCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAGACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATg                                                                                     <  2:318333/285‑1 (MQ=255)
                                                                                                                                                                                                     aGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAg                                                                                                   <  1:69998/273‑1 (MQ=255)
                                                                                                                                                                                                     aGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAg                                                                                                   >  2:69998/1‑273 (MQ=255)
                                                                                                                                                                                                               atatCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGAcc                                                                         <  2:305813/289‑1 (MQ=255)
                                                                                                                                                                                                               atatCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGa                                                                           <  2:304756/286‑1 (MQ=255)
                                                                                                                                                                                                               atatCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGa                                                                           >  1:304756/1‑286 (MQ=255)
                                                                                                                                                                                                                          cgcgTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGaaaaa                                                                                                    <  2:231654/250‑1 (MQ=255)
                                                                                                                                                                                                                          cgcgTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGaaaaa                                                                                                    >  1:231654/1‑250 (MQ=255)
                                                                                                                                                                                                                            cgTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGTAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCgg                                                                   >  1:219023/1‑282 (MQ=255)
                                                                                                                                                                                                                            cgTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCgg                                                                   <  2:219023/282‑1 (MQ=255)
                                                                                                                                                                                                                                      ctcGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAg                                                 <  2:303360/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  aTCTGGTGCTCAGCAACGGGCGGCAACTGGAAATTCGCGTCATCCTTATACCCACAAACAGTGGCTGGGGGTGGCGCGGGATGGCACGGAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACg                                     <  2:296585/289‑1 (MQ=255)
                                                                                                                                                                                                                                                         gCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCAAGAGGATGGAAGGACTGGTGAAGCCATTGCTGACGCTGTCGa                              >  2:108864/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                   ttCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGaaaaa                                                                                                    <  1:236614/194‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   ttCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGaaaaa                                                                                                    >  2:236614/1‑194 (MQ=255)
                                                                                                                                                                                                                                                                                     cgcgTCAGCCTGAATCCCCCAAGCAGTGGCGGAGGGTGGAGCGTGATGTGACGAGAATGCATAAACGGGAAGGGACGCGGCGTAACGTATTTGCCAACGTGAGCCATGAGTTACGTACGCCAGTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAAGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCt  <  1:173766/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                         
ATAAAAAACGCCGCCGTGAACTGGGCAATCTGATTAAACGCTTTCGTAGCGGCGCGGAGTCGCTGCCCGACGCGGTGGTGCTGACCACGGAAGAGGGCGGTATTTTCTGGTGTAACGGTCTG‑‑GCGCAACAAATTCTTGGTTTGCGCTGGCCGGAAGATAACGGGCAGAACATCCTTAACCTACTGCGTTACCCGGAGTTTACGCAATATCTGAAAACGCGTGATTTTTCTCGCCCGCTCAATCTGGTGCTCAACACCGGGCGGCATCTGGAAATTCGCGTCATGCCTTATACCCACAAACAGTTGCTGATGGTGGCGCGTGATGTCACGCAAATGCATCAACTGGAAGGGGCGCGGCGTAACTTTTTTGCCAACGTGAGCCATGAGTTACGTACGCCATTGACCGTGTTACAGGGTTACCTGGAGATGATGAATGAGCAGCCGCTGGAAGGCGCGGTACGCGAAAAAGCGTTGCACACCATGCGCGAGCAGACCCAGCGGATGGAAGGACTGGTGAAGCAATTGCTGACGCTGTCGAAAATAGAAGCCGCACCGACGCATTTGCT  >  NC_000913/418145‑418709

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