Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,023	A→G	V282V (GTA→GTG)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,023	0	A	G	97.6%	115.8 / ‑5.7	41	V282V (GTA→GTG)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (16/24);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

ATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATG  >  NC_000913/4638782‑4639266                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    aTTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTGCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGGTTAGCGTCGCAGTGACGACGCCa                                                                                                                                                                                                                                            <  1:378572/251‑1 (MQ=255)     aaGCGGTCCCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGcc                                                                                                                                                                                                                                         <  2:283157/250‑1 (MQ=255)      aGCGGTCCCGGCGATCTTGCTGTGGTGCCGGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTGTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGa                                                                                                                                                                                                                                       <  2:445126/251‑1 (MQ=255)          gttccGGCGATCTTTCTGGGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGGCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATTCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCaa                                                                                                                                                                                                                                   <  1:418483/248‑1 (MQ=255)                       ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTTGCGTCGCCGTGACGCCGGCAGCCGGTCAAttcacattaa                                                                                                                                                                                                                         >  2:165486/1‑239 (MQ=255)                       ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTaa                                                                                                                                                                                                                         <  1:165486/248‑1 (MQ=255)                            gggggCCTGCCGGGCGTAAAAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGATTCAGGTTTTCAGGCGCAATGGAAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGGGACGACGCCAGCCGAGCAATACCAATTAAcggtccgg                                                                                                                                                                                                                 <  2:367211/247‑8 (MQ=255)                                        gggggaaaaGCGAAATGACCCTAACCAGCAACTGGGCGAATCCCAGGTTTTTAGGTGAATTTTTGCTAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAACGTCAAAGCAGCTGGTTTGCTGATAACCTCGGTGAGCGTTTTGTTTCGGGCAATGATACGGGCTGGGAAAACGTCCCCGCGTTTAGCGGCGCAGTGACGAGGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTAc                                                                                                                                                                                                     <  2:45351/243‑1 (MQ=255)                                                                ccAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGa                                                                                                                                                                               >  1:25410/1‑249 (MQ=255)                                                                ccAGCAACTGGCCGCAGCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGa                                                                                                                                                                               <  2:25410/249‑1 (MQ=255)                                                                            cGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGt                                                                                                                                                                                                       >  1:355533/1‑213 (MQ=255)                                                                            cGCATCCCAGTTTTTTAGGTGAGTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTGCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATGCCAATGAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTc                                                                                                                                                                 <  1:223283/251‑1 (MQ=255)                                                                             gCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTAAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGCCTTTTATGGCGTTCt                                                                                                                                                                >  1:476684/1‑251 (MQ=255)                                                                                 cccAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGt                                                                                                                                                            >  1:495646/1‑251 (MQ=255)                                                                                            aGGTGATTTTTTGCTAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAAGGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTATCGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACg                                                                                                                                                  <  1:465693/250‑1 (MQ=255)                                                                                             ggTGATTTTTTGCCAGCCACACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTGGCGTAGCCGTGGCGACGCCAGCCGATCCATACCAAATCACGGAACGGGCGAATTAATACCACACTTTGCTGATTTCACTGGCTTTTCTGGCGTTCTTTTTTTTTGAAAcacc                                                                                                                                                >  2:61603/1‑248 (MQ=255)                                                                                              gtgagttgttgccAGCCAACCGGGAAGTTAGCGAGTCAGGTTTGCCGGCGCAAGGGCAAAGCAGCTGGTTGGTTAATAATTGCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTGGAAACGCTc                                                                                                                                               <  2:149178/240‑1 (MQ=255)                                                                                                       tgccAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGGCTAAGTACGCCCTTCTGCTGATTTCCATGGCTTTTATGGCGTTCTTTTTTTTTGAAACGCTCACCTCGCa                                                                                                                                       >  2:229749/3‑250 (MQ=255)                                                                                                           aGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                                                                                                                                  <  1:211303/251‑1 (MQ=255)                                                                                                                     gAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAAt                                                                                                                        >  1:40980/1‑251 (MQ=255)                                                                                                                          tAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAat                                                                                                                   >  1:331507/1‑251 (MQ=255)                                                                                                                                         tCAGGCGTAATGGCAAAGCAGCTGGTTTGCTAATAATTTCGGTGGGCGTTTTGCTACAGGCAAGGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTgggg                                                                                                     <  2:274168/250‑1 (MQ=255)                                                                                                                                              cgcAATGGCAAAGCAGCTGGTTTGCTACTAATCTCGGTGAGCGTTTTGCTTCAGGCACTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCCTTCTGCTGATTGCACTGCCTTTTATGGCGTTCTTTGTTTTTTAAACGCTCACCGCGCAACGTTTAAACCCAATGCAATATTTGCTGGTGGGGCTTTc                                                                                                >  1:163111/1‑250 (MQ=255)                                                                                                                                                      cAAAGCAGCTGGTTTTCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCTCCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGAt                                                                                       <  1:242792/251‑1 (MQ=255)                                                                                                                                                      cAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGAt                                                                                       >  1:38320/1‑251 (MQ=255)                                                                                                                                                       aaaGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATg                                                                                      <  1:324496/251‑1 (MQ=255)                                                                                                                                                            agcTGGCTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAACCTTCCCCGCGTTTAGCGGCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTa                                                                                 <  1:266253/251‑1 (MQ=255)                                                                                                                                                                ggTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTAttt                                                                              <  1:229749/250‑1 (MQ=255)                                                                                                                                                                        aataatCGCGGTGGGCGTTTTTCTTCAGGCAAAGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGGCGCAGTGACGACCCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGc                                                                      <  2:390238/250‑1 (MQ=255)                                                                                                                                                                          taatCTCGGTGAGCGTTGTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCCGCCAGCCGATCAATACCAATTAACGGACCGCGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAAAGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTTATTGGTACTGTTTTATTTGCTCTTGATgg                                                                   >  2:101359/1‑251 (MQ=255)                                                                                                                                                                                 ggTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCCAACCGCGCAACGTTTACACCCAATGCCATATTTGTTGGTGGGGCTTTCATTGGGGGTGTTTTATTTGCTCTTGCTGTCGCTTTc                                                            >  2:314011/1‑251 (MQ=255)                                                                                                                                                                                           ttGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATat                                                   <  2:210857/250‑2 (MQ=255)                                                                                                                                                                                                tcaagCAATGATACCGGCTGGGAAAACTTCCCCCCGGTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCAtt                                                                                                                                                                                               >  2:393248/5‑105 (MQ=255)                                                                                                                                                                                                tCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCAtt                                                                                                                                                                                               <  1:393248/105‑1 (MQ=255)                                                                                                                                                                                                      aaTGTTACCGGATGGGAAACCTCCCCCGCGTTTAGCGGCGCAGTGACGTCCCCAGCCGATCAATACCTATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGt                                       <  2:382535/251‑1 (MQ=255)                                                                                                                                                                                                            aCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGGGGGGCTTTCATTGGTGGTGTTTTATTTGCTCTTGATGGAGCTTTCTGAACATATCGGTTTTACCGTGGCATg                                 >  1:16395/1‑251 (MQ=255)                                                                                                                                                                                                              cGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATgg                                >  1:416810/1‑250 (MQ=255)                                                                                                                                                                                                                     gAAAACTTCCCCGCGTTTAGCGGCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGc                        <  1:83155/251‑1 (MQ=255)                                                                                                                                                                                                                       aaaCTTCCCAGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCa                      <  2:33660/251‑1 (MQ=255)                                                                                                                                                                                                                                     ttAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCTCCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGTTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGc         <  1:128913/250‑1 (MQ=255)                                                                                                                                                                                                                                            cGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAATg  <  1:251357/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    ATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAATG  >  NC_000913/4638782‑4639266

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