Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,608	G→A	C233C (TGC→TGT)	arcA ←	DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,608	0	G	A	97.9%	143.8 / ‑4.3	48	C233C (TGC→TGT)	arcA	DNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/24);  total (24/24)
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

TCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCC  >  NC_000913/4639364‑4639850                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ccaattctGCCGAGAGCGCGGTGTTGTGGGGAACCAGTGTGCTGGTGGTGGCGCTTGCCGGCATGATGTTTGTGAATCGTAATATCGACTGGTAAGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGCCGATCAATTTCGTAGCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGGTTTTGTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCTCCACAGAAg                                                                                                                                                                                                                                                 <  2:428324/246‑1 (MQ=255)         gCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAAcc                                                                                                                                                                                                                                         <  1:41578/251‑1 (MQ=255)                gcgcgcTTTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCg                                                                                                                                                                                                                                  <  1:474040/251‑1 (MQ=255)                                 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)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTC‑GTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCC  >  NC_000913/4639364‑4639850

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: ‑