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,6320AG97.9% 148.1 / ‑4.4 48I225I (ATT→ATC) arcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (25/22);  total (26/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 = 5.75e-01

ACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATG  >  NC_000913/4639396‑4639872
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
actagtgttgtgggggtggtGCTGGCCGGCAGGATGTTTGTGCAGCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCtatgat                                                                                                                                                                                                                                       <  2:171693/237‑1 (MQ=255)
      gtgCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTtccgg                                                                                                                                                                                                                                 <  1:42172/251‑1 (MQ=255)
             tggtggCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTtccggcgtatcc                                                                                                                                                                                                                          >  1:15708/1‑251 (MQ=255)
                                          actcGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcg                                                                                                                                                                                             >  1:280951/4‑251 (MQ=255)
                                            tcGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggc                                                                                                                                                                                           >  1:613519/2‑251 (MQ=255)
                                             cGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggcg                                                                                                                                                                                          <  1:262603/251‑1 (MQ=255)
                                                 atatCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATc                                                                                                                                                                                      <  2:641644/251‑1 (MQ=255)
                                                     cGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTc                                                                                                                                                                                   <  2:167574/250‑1 (MQ=255)
                                                        cTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGAGGCGATGATTTCCGGCGTATCAGGCGTAGCTTCGAATTGTTTACGAATACGGCGGACCGTccacc                                                                                                                                                                               >  2:499732/1‑247 (MQ=255)
                                                           gTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCt                                                                                                                                                                             >  1:244225/1‑250 (MQ=255)
                                                                cGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAACTACGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCCGATCACCACCGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAGTGTTTACGAATACGGCGGATCgcaccggctcccgt                                                                                                                                                                        >  2:241149/1‑237 (MQ=255)
                                                                   tttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACg                                                                                                                                                                     <  1:207851/250‑1 (MQ=255)
                                                                          gCCGAAAATGAAAGCCAGGAAAGAAGTTACAATGGACGATTAATTACGTATCTGGAAATAAGGGTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGACTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGAGTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAAGACGGCGGATCGTCACGTCTACAGTACGGTC‑GTgtg                                                                                                                                                            <  2:218694/251‑3 (MQ=255)
                                                                            cGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGt                                                                                                                                                          <  1:209645/251‑1 (MQ=255)
                                                                              aaaaTGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCAACGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTCCAGTC‑GTGAGGttt                                                                                                                                                        >  1:655240/1‑251 (MQ=255)
                                                                                aaTGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGCTCGTCACGTCTACCGTACGGTC‑GTGTGGTTTc                                                                                                                                                       >  1:680361/1‑250 (MQ=255)
                                                                                  tGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGGTTCGAAATGTTTACGAATACGGCGGATCGTCACCTCTACAGTACGGGT‑GTGGGGTTTCatg                                                                                                                                                    >  2:355967/1‑249 (MQ=255)
                                                                                  tGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GGGTGGTTTTAGc                                                                                                                                                    >  2:231789/1‑251 (MQ=255)
                                                                                    aaaGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGCAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGTTGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACCATACGGTcagtgggttttaaccc                                                                                                                                                   >  2:158408/1‑236 (MQ=255)
                                                                                     aaGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCa                                                                                                                                                 <  1:356523/251‑1 (MQ=255)
                                                                                     aaGCCAGTAAAGAAGGTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCa                                                                                                                                                 <  2:613519/251‑1 (MQ=255)
                                                                                                           ggACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTCAAAAGCGCCGTTTTTTTTGACGGTGGAAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCAAGTCAACAGTACGGTC‑GTGTGGATTCAGCTAACGGCCGG‑TCATTttctttggcc                                                                                                                          >  1:430142/1‑246 (MQ=255)
                                                                                                              cgatcaattACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGtt                                                                                                                       <  2:411397/244‑1 (MQ=255)
                                                                                                              cgatcaattACGTAGCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGtt                                                                                                                       <  1:524876/244‑1 (MQ=255)
                                                                                                                 tcaattACGTATCTGGAAATAAGGTTGAAACATAAAAACGGCGCTAAAACGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGAGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGAGTAGATTCGAAATGTTTACGAATACGGCGGCCAGTCACGTCTTCAGTCCGGTCTGTG‑GGTTTCAGCTCACGGCCGA‑TATTTTTCTTcagaggttca                                                                                                                     >  2:314456/5‑243 (MQ=255)
                                                                                                                     ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACg                                                                                                                >  1:418390/1‑251 (MQ=255)
                                                                                                                        cGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGTTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACgaga                                                                                                             <  1:175628/251‑4 (MQ=255)
                                                                                                                                gAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAAttt                                                                                                     <  1:309043/251‑1 (MQ=255)
                                                                                                                                       ggTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑gg                                                                                             <  1:644765/251‑1 (MQ=255)
                                                                                                                                           gAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑GCATTTTCTTCAGCAGTTCAGCCCGCGGCTGAATTTTTACT‑GGGttt                                                                                         >  2:184132/1‑251 (MQ=255)
                                                                                                                                           aaaaaaTAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGttt                                                                                         <  2:144326/250‑1 (MQ=255)
                                                                                                                                                aTAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGGGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTcac                                                                                     >  1:554439/1‑250 (MQ=255)
                                                                                                                                                   aaaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGGCTACAGGACGGGC‑GTGGGGTTTCATCCCACGGCCGGATC‑TTTTCTTCCGCCGTTCAGCCCGGGACTGAATTTTGGCT‑TGGTTTTcaccga                                                                                  >  2:134791/1‑247 (MQ=255)
                                                                                                                                                              aaaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATgg                                                                      <  1:507914/251‑1 (MQ=255)
                                                                                                                                                               aaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTgtgggttc                                                                                                                                                       >  2:240354/1‑164 (MQ=255)
                                                                                                                                                               aaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTc                                                                                                                                                       <  1:240354/171‑1 (MQ=255)
                                                                                                                                                                 aaGCGCCGTTTTTTTTGCCGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTAGGGTTTT‑ACAGAAGTGAAGCATGgcgc                                                                   >  1:537919/1‑251 (MQ=255)
                                                                                                                                                                    cgccggtgttttttACGGTGTTACAGCCGTTTAATTTGCCAGGTCACCACAGCATCGTTATCCCTCACCTTGGATGGTGGCGATGATTTTCGGCGTATCCGGCGTAGATTCAAAATTTTTCCTAATCCGGCGGATCGTCCCGTCTTCAGTACTGTC‑CTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GCGTTTTCACTGAAGTGAAGCATGGCGCGGa                                                                <  2:176875/243‑1 (MQ=255)
                                                                                                                                                                                    ggtggtAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGGTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGGGAAGCATGGCGCGGAACTAGCTGggcggc                                                  >  2:524625/1‑249 (MQ=255)
                                                                                                                                                                                            aGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTc                                                                                                                                                  <  1:364380/147‑1 (MQ=255)
                                                                                                                                                                                               cGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTATATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTCTCACAGAAGTGAAGCATg                                                                       <  2:97576/217‑1 (MQ=255)
                                                                                                                                                                                               cGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATg                                                                       >  1:97576/1‑217 (MQ=255)
                                                                                                                                                                                                gATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTc                                    >  1:331189/1‑251 (MQ=255)
                                                                                                                                                                                                   tAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGGCCCGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTTAGCAGTTCAGGACGAGACTGAATTTTTCCT‑GGGTTTTTACCGGAGGTAAGCATGTCGCGGGACACACTGCGCGGcggtttttacgcctcgcc                                 >  2:418005/1‑234 (MQ=255)
                                                                                                                                                                                                            aGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGcc                        >  1:641416/1‑251 (MQ=255)
                                                                                                                                                                                                                    acAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCaacg                <  1:484887/250‑1 (MQ=255)
                                                                                                                                                                                                                            gATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGACACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCCAGCGACTGAATTTTGCCT‑GGGTTTTCCCAGAAATGAAGCATGGCGCGGGACTCGCTGCGCGGCCGCTTGTACTGCTCGCCCTCAGGGCCGGTCAACGCACGGCTg        >  1:375605/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATg  >  1:315858/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
ACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATG  >  NC_000913/4639396‑4639872

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