Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,834 A→G intergenic (+117/‑28) bdcR → / → tabA putative transcriptional regulator BdcR/DUF386 domain‑containing toxin‑antitoxin biofilm protein TabA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,8340AG97.6% 109.4 / NA 42intergenic (+117/‑28)bdcR/tabAputative transcriptional regulator BdcR/DUF386 domain‑containing toxin‑antitoxin biofilm protein TabA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (18/23);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGA  >  NC_000913/4474589‑4475075
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
gtgACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttt                                                                                                                                                                                                                                              >  1:316734/1‑250 (MQ=255)
 tgACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAACCCCTGCCGCTAAGTTGTGATAAGACAATAAATCACGCGTTAAAGCTAGCGGGATTGATTTTCTTTTTTTATAACGATAGAAATTAAGCgttgt                                                                                                                                                                                                                                             >  2:133002/1‑251 (MQ=255)
    cTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAATCTTGAGGAGTGATGGGGGCGTTGGCCCGAAAGGCGGGAAAGGGCCTGCCGGTAAGTTGGGATAAGACAATAATTCACGCGTTAAGGGTAGCGGAATTGATTATCTTTTCGTATAAAGAAAGAAATGGAGCGGTGtttt                                                                                                                                                                                                                                          >  1:409237/1‑251 (MQ=255)
      gATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaa                                                                                                                                                                                                                                        <  1:343883/251‑1 (MQ=255)
      gATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCGATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGtttt                                                                                                                                                                                                                                          <  2:87670/249‑1 (MQ=255)
           tATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGACAGGCGGGAAAGCGCCTGCCGGTAAGTTGTGATAAGACAATAATTTACGCGTTAAAGCTAGCGGAATTGGTTATCTTTTCGGATAAAGAAAGAAATGGAGCGGTGTTttaattt                                                                                                                                                                                                                                     >  2:53293/1‑248 (MQ=255)
              gAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGGTAAGACAATAATTCCCGCGTTAAGGCTAGCGGAATTTATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAgg                                                                                                                                                                                                                                  >  1:84583/1‑249 (MQ=255)
                     gcgATGTCTGGGCTTTCTGCGAAGGCACGAGCGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATCGCGGGAGACAAGATCCAAACTATTCTTGTGGAATGATGCGGGCCCTGCCCCGGAAGGCGGGCAAGCGCCTGCCGCAAAGGTGTGATAAAACAACAATTCACGGGCTAAGGCTCGCGGGATTGATTCTCTTTTCGTATAAACATAGAAATGAAGCGGTGGTTTTATTAAgacggggaaa                                                                                                                                                                                                                          >  2:388002/3‑241 (MQ=255)
                             tGGGCTTTCTGAGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCCATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGGCCCGGAAGGGGGGAAAGCGCCTGCCGATAAGTTTTGATAAGACAATAAATCCCAGGTTAAAGCCAGCGGAATTGATTTTATTTTCGTATAACGCTAGAAATGAAGAGTTGTTTTAATTATGGAGTTGAAACgatgtgcc                                                                                                                                                                                                                 >  2:79898/1‑247 (MQ=255)
                                       gcaaaGGCACGAGAGGGGCACTCTATAGAACAACTCGTGGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGAGGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGTAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                                                                                                                                                                                                          <  2:264531/245‑1 (MQ=255)
                                           aGGCACGAGAGGGGCACTCTATAAAACAACTCTGTGCAACAGCTGCAATGTCGGGGGAAGCGAGAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGTAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATTAAGCGTTGTTTTAATTAAGGATTTGAAAAGATGATCATCGGAAATATTCAt                                                                                                                                                                                                   <  2:30106/251‑1 (MQ=255)
                                             gcaagagaGGGGCACTCTATAGAACAACTCTGTGCAAAAGCTGCAATGGGGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATaa                                                                                                                                                                                                 <  2:291700/247‑1 (MQ=255)
                                                     ggggCACTCTATAGACCAACTCTGGGCAAACGCTGCAATGGCGGGGGAAGCAATAAAAACTATTCTTGAGGAGTGATGCGGGTTTTGAACCGAAAGGCGGGAAAGCGCCTGCCGATATGTTGTGATAAGACAATAATACACGCGTTAAGGCTAGCGGAATTGATTAGCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGc                                                                                                                                                                                         <  2:305011/251‑1 (MQ=255)
                                                       ggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGTGGGGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGGTAAGTTGTGATAAGACAATAAATCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAAAGATAGAAATGAAGCGGTGTTTTTATTAAAGAGTGGGAAAGATGATCATCGGAAATATTCATAATCTTCAGCCg                                                                                                                                                                                       >  2:162995/1‑251 (MQ=255)
                                                              tataGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGGGATAAAAACTATTTTTGAGGGGTGATGCGGGGCTTGAACCGAAAGGCGGGGAAGCGCCCGCCGATAAGGTGGGATAAAACAATAAATTACGGGTTAAGGGTAGCGGGAATGTTTTACTTTTCGAATAACGAAAGAAATGAAGCGGTGTTTTTATTAAGGAGTGGAAAAGATGATCAACGGAAAAATTCATAATCTTCATCCGTGGCTAc                                                                                                                                                                                >  2:1520/1‑251 (MQ=255)
                                                                taGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTTCCGATAAGTTTGGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGGTTATCTTTTCGTATAAAGATAGAAATGAAGCGGTGTTTTAATTTAGGAGTGGAAAAGATGATCATCGGAAATATTTATAATTTTCAACCGTGGCTcccc                                                                                                                                                                              >  1:317603/1‑251 (MQ=255)
                                                                 aGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGACGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAcccc                                                                                                                                                                             >  1:143853/1‑251 (MQ=255)
                                                                              gCAACAGATGCAATGGCGGGGGAAGCGATAAAAACTAGTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCa                                                                                                                                                                <  2:316734/251‑1 (MQ=255)
                                                                                      cTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGGGGCTACCCCAGGAGTTACGCCAGGCGAtt                                                                                                                                                         >  1:168679/1‑251 (MQ=255)
                                                                                       tGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGGAAGGCGGGGAAGCGCCTGCCGCTAAGATTTGGTAAGACAATAAATCACGCGGTAAGGCTCGCGGGATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTTATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCCGCCGTGGGCACCCCCGGAGTTACGCCCGGCgcttt                                                                                                                                                        >  2:99628/1‑247 (MQ=255)
                                                                                          aaTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGc                                                                                                                                                     <  1:99628/251‑1 (MQ=255)
                                                                                                gggggAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATc                                                                                                                                                >  1:67871/1‑250 (MQ=255)
                                                                                                        cGATAAAAACTATGCTTGGGGATGGAGGGGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCacac                                                                                                                                       <  2:9620/251‑1 (MQ=255)
                                                                                                                                   gCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATACGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAAATGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGACCAAGGGCa                                                                                                              >  1:254949/1‑249 (MQ=255)
                                                                                                                                         cTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATAAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACgct                                                                                                      >  1:243455/1‑248 (MQ=255)
                                                                                                                                                 gAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAgg                                                                                              <  1:162995/251‑1 (MQ=255)
                                                                                                                                                            aaGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGt                                                                                   >  1:183790/1‑251 (MQ=255)
                                                                                                                                                                    gCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTGTCGTATAACGATAGAAAGGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTttatct                                                                            <  2:137014/250‑1 (MQ=255)
                                                                                                                                                                     ccGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGTACGATATCGAAGGCAATCGACTGTTttatctta                                                                          <  1:388002/251‑1 (MQ=255)
                                                                                                                                                                     ccGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTttatctta                                                                          <  1:133002/251‑1 (MQ=255)
                                                                                                                                                                      cGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGTAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTttatcttat                                                                         <  2:315032/251‑1 (MQ=255)
                                                                                                                                                                                  aTAAGACAATAATTCACGCGTGAAGGCTAGCGGAATTGATTATCTTGTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAATGAGTGGAAAAGATGATCATCGGAAATATTCATAATCGTCAGCCGTGGCTCCCCCAGGAGTTACGCCAGTCGATTGTGCATAGCAAAGCACACGTTACGGCAGAACCGCCAAAGGGCAAGCACGATATCGAATGCAATCGACTGTTttatcttatcgctgtagatat                                                             <  2:399453/251‑12 (MQ=255)
                                                                                                                                                                                                        aaGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGc                                       >  1:409399/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTa                            <  2:143853/251‑1 (MQ=255)
                                                                                                                                                                                                                       gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCAt                        <  1:53293/251‑1 (MQ=255)
                                                                                                                                                                                                                               tttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCt                 <  2:67871/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACAtt      <  2:84583/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGGCATTc     <  1:324296/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGtt                                                                                                                                                                      <  1:367343/88‑1 (MQ=255)
                                                                                                                                                                                                                                            tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGtt                                                                                                                                                                      >  2:367343/1‑88 (MQ=255)
                                                                                                                                                                                                                                             aGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAg   <  1:145381/250‑1 (MQ=255)
                                                                                                                                                                                                                                             aGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTAGCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTGGACATTCAGa  <  2:169037/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
GTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCATATCAAAGCACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGA  >  NC_000913/4474589‑4475075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: