Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,719 C→T P196P (CCG→CCA)  arcA ← DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,7190CT89.4% 103.7 / 4.0 49P196P (CCG→CCA) arcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (28/14);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

GAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGT  >  NC_000913/4639473‑4639965
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
gAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACCGCTACAGTACGGGCGTGTGGtt                                                                                                                                                                                                                                                       >  2:484733/1‑251 (MQ=255)
     tGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGTTAAAAAACGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATTTTCCAGCTCACCACACAAGCGTTAACCTTCACCGTGGATGGTGGCGCTGATTTCCGGCGTATCCGGAGTAGTTTCGAATTGTTTACGAATACGGCGGATCGTCTCGTCCAAAGTATGGACGTGTGGTTTCAGc                                                                                                                                                                                                                                                  >  1:499163/1‑251 (MQ=255)
                    aagTTACAATGGACGATCAATGACGTATCTGGAAATAAGGGTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTGTTGTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGCAGCGAGATCCTTCACCGTGGATGGTGGCGATGAGTGCCGGCGTATTCGGCGGAGATTCGAAATGTTTACGAATACGGCGGATCGACACGTCTTCAGTACGGTCGTGCGGTTTCGGCGTACTGCCGGTCAtt                                                                                                                                                                                                                                   <  2:57454/251‑1 (MQ=255)
                         aCAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTttctt                                                                                                                                                                                                                              <  1:185456/251‑1 (MQ=255)
                                 cgatcaattACGTAACTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGtt                                                                                                                                                                                                                      <  1:245464/244‑1 (MQ=255)
                                      aattACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAACGTCAAGTATACAGAACGGTCGTGTGGTTTAAGCTCACGGCCGGTCATTTTCTTCAGCAGTTAAGc                                                                                                                                                                                                                  >  1:153002/3‑250 (MQ=255)
                                         tACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCCTTTTCTTCCGCAGTTccgcccgg                                                                                                                                                                                                              >  2:37466/1‑244 (MQ=255)
                                          aCGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACg                                                                                                                                                                                                               >  1:314562/1‑249 (MQ=255)
                                                  ggAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCTCGTCTACAGTACGGTCGTGTGGTTTCAGCTCACTGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAAt                                                                                                                                                                                                      <  1:497951/250‑1 (MQ=255)
                                                        aaGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGAGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTAAGGTCGTGTGGTTTCAGCTCACGGCCGGGCATTTTCTTCAGCAGTTcctccagagactgaattttggct                                                                                                                                                                                               >  1:196767/1‑229 (MQ=255)
                                                         agggtGAAAAATAAATACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCGTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATTTTTACGAATACGGCGGATCGTCACGTCTACAGTACGTTCGTGTGGTTTCAGCTCTCGGCCGGTCAGTTGCTGCAGCAGTTCAGCACGAGACTGAATTTTGCCTg                                                                                                                                                                                              <  1:187269/247‑1 (MQ=255)
                                                              gAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGttt                                                                                                                                                                                         >  1:482474/1‑251 (MQ=255)
                                                              gAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGttt                                                                                                                                                                                         <  1:114320/251‑1 (MQ=255)
                                                                       aaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAg                                                                                                                                                                                <  1:37466/251‑1 (MQ=255)
                                                                           ggCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCTTGGGTTTTCACAGAAGTGaa                                                                                                                                                                            <  1:63999/251‑1 (MQ=255)
                                                                              gcTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGc                                                                                                                                                                          >  1:264188/1‑250 (MQ=255)
                                                                                tAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACCGAAGTGAAGCAT‑g                                                                                                                                                                      >  1:217239/1‑251 (MQ=255)
                                                                                    aaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑Ggcgc                                                                                                                                                                  >  1:271737/1‑251 (MQ=255)
                                                                                                ttttGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTCCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTg                                                                                                                                                       <  1:190141/250‑1 (MQ=255)
                                                                                                ttttGAAGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTACCGTCTACAGTAAGGTCGTGTGGGTTTAGCTCACGGACGGTCCTTTTCTTCAGCCGATCAGCACCAGGCTGAACTTTGCCTGGGGTTTTACCAAAGTGAAGCAT‑ggaggggaaatccctgc                                                                                                                                                      >  2:473346/1‑236 (MQ=255)
                                                                                                    gACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCgg                                                                                                                                                  >  1:119919/1‑251 (MQ=255)
                                                                                                          ggtAAAGCCGATTAATCTTCCAGACCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCt                                                                                                                                             <  2:406261/250‑1 (MQ=255)
                                                                                                             aaaGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGt                                                                                                                                          >  1:326595/1‑250 (MQ=255)
                                                                                                               ggCCGATTAATCTGCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACt                                                                                                                                       <  2:119919/250‑1 (MQ=255)
                                                                                                                          tttCCAGATCCCCCCAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATc                                                                                                                             <  2:271737/249‑1 (MQ=255)
                                                                                                                                    accacAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCa                                                                                                                 <  1:484733/247‑1 (MQ=255)
                                                                                                                                       acAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCCCAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCaacg                                                                                                              >  1:495230/2‑251 (MQ=255)
                                                                                                                                        cAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCTCAGAAGTGAAGCAT‑GGCGCGGAACTCGCTTCGCGGCAGCTTGTACTTCTCGCCATCA‑GGGCCGATCaaca                                                                                                              >  1:69549/1‑249 (MQ=255)
                                                                                                                                           aaGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGGGAAGCAT‑GGCGCGGAACTCTCTGCGCGGCAGCTTGTTCTGCTCGCCATCA‑GGGCCGCTCaaagaacg                                                                                                          >  2:493189/1‑251 (MQ=255)
                                                                                                                                              cGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGACTACGGCGGATCGTCACGTCTACAGTACGGGCGGGTGGGTTCAGCTCACGGGCGGTCATTTTCTTCAGCAGTTCCGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATTGGC‑CGGAACTCGCTGCGCGGCAACTTGTACTTCCCGCCATCT‑GGGCCGATCaccgaaaggat                                                                                                       >  2:297535/1‑246 (MQ=255)
                                                                                                                                                 aaaCCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGtt                                                                                                    <  2:396217/250‑1 (MQ=255)
                                                                                                                                                      ttCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCg                                                                                                                                                                                                                                         >  1:146880/1‑115 (MQ=255)
                                                                                                                                                      ttCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCg                                                                                                                                                                                                                                         <  2:146880/115‑1 (MQ=255)
                                                                                                                                                        cACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAACGTCACGTCTCAAGTACGGTCGTGTGGTTTCAACTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGGTTTCACAGAAGGGAAGCCT‑GGAGCGGAACTCGCTGCGCGGCAGCTTGTTCTGCTCGCCCTAC‑GGGCAGATCACAGAACGGCTGTTGgtgttc                                                                                             >  2:423915/1‑249 (MQ=255)
                                                                                                                                                                                ccggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑t                                                                      >  1:293955/1‑249 (MQ=255)
                                                                                                                                                                                 cggcgtatccggcgtaGATTCGCAATGTTTACGAATACGGCGGCTCGTCACGTCTACAGTACGGGCGTGTGGTTTCAGCTCACGGCCGGTCCTTTTCTTCAGAAGTTCAGCACAGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGGACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCAGTTGATGTCCACGTCCCCACCATACAACTTG‑TAGc                                                                   >  2:309927/1‑251 (MQ=255)
                                                                                                                                                                                                          tGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCCGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTAACAGAAGTGAAGCAT‑GGCGCGGCACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGACGTCCAGTTCCCAACCATTGTACTTG‑TAGCTTTCAACGCTACGAAGttcttccct                                          >  1:281554/1‑248 (MQ=255)
                                                                                                                                                                                                              tACGAATACGGCGGATCGTCCCGTCTACAGAACGGTCGGGTGGTTTCAGATCACGGCCGGTCATTTTCTTCAGCAGTTCCGAAAGAGACTGAGTTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCA‑CATGGGCCGATCAACGAACGGCGGATGATGTCCCGTTCCCCACCATTGAAACTG‑GAGCTGTCAACGCTACGACGTTCTTCGcgggac                                      >  2:341375/1‑246 (MQ=255)
                                                                                                                                                                                                              tACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGAc                                       <  2:159240/250‑1 (MQ=255)
                                                                                                                                                                                                                cGAATACGGCGGATCGTCACGTATACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGCATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGCACGGCTGTTGATGTCCAGTTCCCAACCATTAAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGt                                    >  1:460936/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTAcc                                 >  1:460011/1‑251 (MQ=255)
                                                                                                                                                                                                                           gATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTc                          >  1:290659/1‑250 (MQ=255)
                                                                                                                                                                                                                               gTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGt                     >  1:305426/1‑251 (MQ=255)
                                                                                                                                                                                                                                tCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGCTTCATGGt                     >  1:187402/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aCAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGAc              >  1:182848/1‑249 (MQ=255)
                                                                                                                                                                                                                                        aCAGTACGGACGTGAGGATTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGATTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTAGCCATAA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTACACACAATTGAACTTGATA‑CTTTAAAAGCTACGGCGTTCTTAGCTGACAGTACCCAGATTCCTGGGACgagtc              >  2:331660/1‑245 (MQ=255)
                                                                                                                                                                                                                                         aaGTACGGGCGTGGGGGTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAgc           <  2:182848/250‑2 (MQ=255)
                                                                                                                                                                                                                                               ggTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCATCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCTATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGTCAGCAGGTTg     <  1:371925/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGACCGGCTGTTGATTTCCAGTTCCCAACCGTTGAACTTG‑TTGCTTTCAACGCTACGACGTTCTTCGCTGACCGTACCCATATTCATGGTACGGGACAGCCGGTTGCgg  >  2:115121/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
GAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCAT‑GGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCA‑GGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTG‑TAGCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGT  >  NC_000913/4639473‑4639965

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