Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,639,632 A→G I225I (ATT→ATC 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,6320AG88.9% 106.4 / 5.2 45I225I (ATT→ATCarcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (23/17);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCC  >  NC_000913/4639392‑4639875
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
gggAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTTGCTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAATTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCga                                                                                                                                                                                                                                            <  1:286148/251‑1 (MQ=255)
      cAGTGTGGTGGGGGGGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGAttt                                                                                                                                                                                                                                      <  2:208651/251‑1 (MQ=255)
                                  atggtgTTTGTGCCTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGTTCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAACCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGTTGGTGGCGATGATTTCCGGCGTATCCGGGGTAGATTCGAAATg                                                                                                                                                                                                          >  2:290076/5‑251 (MQ=255)
                                  atgatgTTTGGGAAACATAATAACGAATGGGATGCGTTTTCACTGCCGAAAAGGAAAGCCAGTAAAGAAAGTAAAATGGACGATCAATTACGTAACTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATg                                                                                                                                                                                                          <  2:480687/251‑1 (MQ=255)
                                     atgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTGCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGttt                                                                                                                                                                                                       <  2:246786/251‑1 (MQ=255)
                                                       atCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGt                                                                                                                                                                                     >  1:374085/1‑251 (MQ=255)
                                                         cGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGGGGAAAAGCCGATTATTCTTCCAGTTCACCCCAGAGCCGATAACCTTCCCCGTGAATGGGGGCGTTGATTTCCGCCGTACCCGCCGTAGATTCGAATGGTTTAGGAATCCGCCGGATCgccc                                                                                                                                                                                   >  1:396217/1‑248 (MQ=255)
                                                                tATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTAc                                                                                                                                                                            >  1:6052/1‑251 (MQ=255)
                                                                   gCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAg                                                                                                                                                                          <  2:158824/250‑1 (MQ=255)
                                                                         tCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTc                                                                                                                                                                   >  1:115596/1‑251 (MQ=255)
                                                                             tGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGCCGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAAGAATACGGCGGCGCGTCACGGCTACAGTAAGGTCgggg                                                                                                                                                               >  2:418837/1‑248 (MQ=255)
                                                                                 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:114320/251‑1 (MQ=255)
                                                                                                                                               gAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGttt                                                                                             >  1:482474/1‑251 (MQ=255)
                                                                                                                                                        aaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAg                                                                                    <  1:37466/251‑1 (MQ=255)
                                                                                                                                                            ggCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCTTGGGTTTTCACAGAAGTGaa                                                                                <  1:63999/251‑1 (MQ=255)
                                                                                                                                                               gcTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGc                                                                              >  1:264188/1‑250 (MQ=255)
                                                                                                                                                                 tAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACCGAAGTGAAGCATg                                                                          >  1:217239/1‑251 (MQ=255)
                                                                                                                                                                     aaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGgcgc                                                                      >  1:271737/1‑251 (MQ=255)
                                                                                                                                                                                 ttttGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTCCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTg                                                           <  1:190141/250‑1 (MQ=255)
                                                                                                                                                                                 ttttGAAGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTACCGTCTACAGTAAGGTCGTGTGGGTTTAGCTCACGGACGGTCCTTTTCTTCAGCCGATCAGCACCAGGCTGAACTTTGCCTGGGGTTTTACCAAAGTGAAGCATggaggggaaatccctgc                                                          >  2:473346/1‑236 (MQ=255)
                                                                                                                                                                                     gACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCgg                                                      >  1:119919/1‑251 (MQ=255)
                                                                                                                                                                                           ggtAAAGCCGATTAATCTTCCAGACCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCt                                                 <  2:406261/250‑1 (MQ=255)
                                                                                                                                                                                              aaaGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGt                                              >  1:326595/1‑250 (MQ=255)
                                                                                                                                                                                                ggCCGATTAATCTGCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACt                                           <  2:119919/250‑1 (MQ=255)
                                                                                                                                                                                                           tttCCAGATCCCCCCAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATc                                 <  2:271737/249‑1 (MQ=255)
                                                                                                                                                                                                                     accacAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCa                      <  1:484733/247‑1 (MQ=255)
                                                                                                                                                                                                                        acAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCCCAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCaacg                   >  1:495230/2‑251 (MQ=255)
                                                                                                                                                                                                                         cAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCTCAGAAGTGAAGCATGGCGCGGAACTCGCTTCGCGGCAGCTTGTACTTCTCGCCATCAGGGCCGATCaaca                   >  1:69549/1‑249 (MQ=255)
                                                                                                                                                                                                                            aaGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGGGAAGCATGGCGCGGAACTCTCTGCGCGGCAGCTTGTTCTGCTCGCCATCAGGGCCGCTCaaagaacg               >  2:493189/1‑251 (MQ=255)
                                                                                                                                                                                                                               cGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGACTACGGCGGATCGTCACGTCTACAGTACGGGCGGGTGGGTTCAGCTCACGGGCGGTCATTTTCTTCAGCAGTTCCGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATTGGCCGGAACTCGCTGCGCGGCAACTTGTACTTCCCGCCATCTGGGCCGATCaccgaaaggat            >  2:297535/1‑246 (MQ=255)
                                                                                                                                                                                                                                  aaaCCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGtt         <  2:396217/250‑1 (MQ=255)
                                                                                                                                                                                                                                       ttCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCg                                                                                                                                             <  2:146880/115‑1 (MQ=255)
                                                                                                                                                                                                                                       ttCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCg                                                                                                                                             >  1:146880/1‑115 (MQ=255)
                                                                                                                                                                                                                                         cACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAACGTCACGTCTCAAGTACGGTCGTGTGGTTTCAACTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGGTTTCACAGAAGGGAAGCCTGGAGCGGAACTCGCTGCGCGGCAGCTTGTTCTGCTCGCCCTACGGGCAGATCACAGAACGGCTGTTGgtgttc  >  2:423915/1‑249 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
GGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTCC  >  NC_000913/4639392‑4639875

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