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	93.3%	116.3 / 0.5	45	V282V (GTA→GTG)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (20/22);  total (21/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 = 3.18e-01

AGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAAT  >  NC_000913/4638787‑4639265                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   aGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGGCGCccgccg                                                                                                                                                                                                                                       >  1:397812/1‑250 (MQ=255)       cccGGCGATCTTTCTGTGGTGCCTCCAGGGTGTAATAGCTAAATGCCCTTAACCAGCAACTGGCCGTATCCCAGTTTTTTAGGTGATTTTTTGCCCGCCAAACGTGAAGTTAGCGAGTCAGTTTTTCAGGCTCAATGTCAAAGCCGCTGGTTTGTTAATAATCTCGGTGAGCTTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTTGCGTCGCAGTGACGACGCCAGCCGATCAAt                                                                                                                                                                                                                                 <  2:524681/249‑1 (MQ=255)            cGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCCACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGGAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCCGCCCGCCCGTCAATACCact                                                                                                                                                                                                                           >  2:351637/1‑249 (MQ=255)                  ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTa                                                                                                                                                                                                                         <  1:20705/247‑1 (MQ=255)                    gtccgatctctgccgggcgTAATAGCGAAATGACCTTAACCAGCAACCGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAAcgga                                                                                                                                                                                                                    <  2:17474/237‑4 (MQ=255)                      gtggtgCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGAc                                                                                                                                                                                                                   <  1:420908/249‑1 (MQ=255)                             ctgccgggcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAAcgga                                                                                                                                                                                                                    >  1:17474/5‑238 (MQ=255)                               gccgggcgTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTa                                                                                                                                                                                                         <  1:342657/248‑1 (MQ=255)                                         acgtgttctcttccggtctcAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGGCGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCAtt                                                                                                                                                                                              <  1:443572/232‑1 (MQ=255)                                           aGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAAt                                                                                                                                                                                                                                 >  1:503138/1‑214 (MQ=255)                                               aaaTGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGAAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATAACGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGAGACGACGCCAGCCGATCAATACCCATTAACGGACCGGGCGAATAAGTACGCCATTctgctg                                                                                                                                                                                        >  1:622077/1‑251 (MQ=255)                                                        tagccAGAAACTGGCCTCATCCCATTTTTTTAGGTGATTTTCTACCATCCAAACGGGAAGGTATCTAGTCAGGGTGTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAAAAATCGCGGTGGGCGTTTGGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCTTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACt                                                                                                                                                                                <  2:649176/247‑1 (MQ=255)                                                            cAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTCCCACCCAACCGGGAGGTTAGCGGGTCAGGTTTTCAGGCGCATTGCCAAAGCAGCTGGTTTGATAATAATCTCGGTGAGCGTTTTGCTTCAGGCAAGGATACCGGCTGGGAAAACTTCCCCGCCTTTAGCGCCCCAGGGCCAACGCCACCCGATCAATACCAAATAAGGGAACGGGCGAATAAGTACGACAtataacaggaagagcgccggg                                                                                                                                                                           >  2:443572/1‑231 (MQ=255)                                                             aGCATCTGGCCGCAGCCCAGTTTTTGAGGTGATTTGTTGCGATCCAAACGGGAAGTGAGCGTGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGGTCGTTTTGCTTCTGGCAATGATACCGGCTGGGAAAACGTCCCCGGGTTTAGCGGCGCAGTGTGGACGTCAGCCGATCAACACCAATTAACGGTCCGGGCGCCTAAGGACGCCATTCTGCTGATTGCACTGACttt                                                                                                                                                                          <  1:415290/251‑1 (MQ=255)                                                                            cccAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTg                                                                                                                                                            <  1:477872/250‑1 (MQ=255)                                                                               aGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTTCTGATTGCCCTGACTTTTTTGGCGTTCTTTGtttt                                                                                                                                                        >  2:220162/1‑251 (MQ=255)                                                                                           gATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCa                                                                                                                                             <  1:219591/250‑1 (MQ=255)                                                                                                    ccAGCCAAACGGGACGTTAGCGAGTCAGGTTTTCCGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGGCGCCGCCAGCCGATCAATCCCCATTCACGGCCAGGGCGCCTAAGAACGCCATTCTGCTGATTGCACTTAATTTTATGGCGTTTTTTGttttttacaagctcacagcgacaac                                                                                                                                   >  2:614152/1‑231 (MQ=255)                                                                                                      aGCCAGACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                                                                                                                                 >  1:86822/1‑251 (MQ=255)                                                                                                      aGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGtt                                                                                                                                 >  1:203731/1‑251 (MQ=255)                                                                                                       gCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGttt                                                                                                                                <  1:351637/251‑1 (MQ=255)                                                                                                              gggAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGCTCAATACCAATTAACGGACCCGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGACAAGCTCACCGCGCAACGTTTACAccc                                                                                                                          >  2:471254/1‑250 (MQ=255)                                                                                                                       gCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCCatat                                                                                                                >  1:400170/1‑251 (MQ=255)                                                                                                                              aGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCtgg                                                                                                         <  1:220162/251‑1 (MQ=255)                                                                                                                                gTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCCATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGCCtt                                                                                                                                                                           >  2:108677/1‑183 (MQ=255)                                                                                                                                gTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACtt                                                                                                                                                                           <  1:108677/183‑1 (MQ=255)                                                                                                                                           caatGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTGCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCAtt                                                                                            <  1:212603/248‑1 (MQ=255)                                                                                                                                            aatGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCCGTGACGACGCCAGCCGATCAATACCAATTAAAGGACCGGGCGGCTAAGTACTCCATTCTGCTGCTTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACCACCACCGCGCAACGTTTTAAACCAATGCAATTTTTTATGGTGGGGCttttcttg                                                                                           >  2:93008/3‑246 (MQ=255)                                                                                                                                             atGGCAAAGAAGCTTGTTTGCTAATAATCGCGGTGCGCGTTTTGCTCTAGGCAATGATACCGGTTGGTAAAACTGCCCCGCTTTTAGCGTCGCAGTGACGCCGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGTTGGGGCTTTCATTgg                                                                                          <  2:510358/250‑1 (MQ=255)                                                                                                                                               ggCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGGACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTg                                                                                        <  1:640611/251‑1 (MQ=255)                                                                                                                                                   aaGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGCAGCCAGCCGATCAATACCAATTAACGGAACGGGCGACTAAATACGCCATTCTGCTGATTGCACTGCCTTTTTTGGCGTTTTTTGTTTTTGAAACGCTCAACGCGCAACGTTTACAACCAATGCAACTTTTTCTGGGGGGGCTTTTCTTGGtgttgt                                                                                    >  2:474330/1‑251 (MQ=255)                                                                                                                                                        gcTGGTTTGCTAATAATCTCGGTGAGCGTTTTGATTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACCCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCCTTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGGTTTTGAAACCCTCACCGCGCAACGTTTACACCCAATGCCATATTTGCTTGTGGGGCTTTCATTGGTGCTGTTTTAt                                                                               >  2:41578/1‑251 (MQ=255)                                                                                                                                                           ggTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTAtt                                                                              <  1:99315/249‑1 (MQ=255)                                                                                                                                                                  taataaTCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCACCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCAATTCTGCTGAATGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCCTTGGTGATGTTTTATTAGATCTTGc                                                                     >  1:374295/1‑251 (MQ=255)                                                                                                                                                                   aataatCTCGGTGATGGTTTTGCTTCAGGCAATGATACCGGCTGGTAAAACTTCCCCGCTTTTAGCGTCGCTGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATAGTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCt                                                                    <  1:327000/251‑1 (MQ=255)                                                                                                                                                                           cGGTTAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAAc                                                                                                                                    >  1:115832/1‑179 (MQ=255)                                                                                                                                                                           cGGTTAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTACCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAAc                                                                                                                                    <  2:115832/179‑1 (MQ=255)                                                                                                                                                                           cGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCttt                                                            <  1:471254/251‑1 (MQ=255)                                                                                                                                                                             gTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCt                                                          >  1:76486/1‑251 (MQ=255)                                                                                                                                                                                         cTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCCACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGt                                              >  1:641437/1‑251 (MQ=255)                                                                                                                                                                                                     aTACCGTCTGGGAAACCTTCCCCGCTTTTGGCTGCGCAGTGACGACTCCAGCCGATCAATACCAATTAACGGACCGGGCTACTAAGTACGCCATTCTGCTGATTGCACTGACTGTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTGACACCCAATGCAATATTTGCTGGTGGGGCTTTCAGTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCa                                  <  2:76486/251‑1 (MQ=255)                                                                                                                                                                                                         cGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGa                              <  1:614152/251‑1 (MQ=255)                                                                                                                                                                                                              gggAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATc                         <  1:430129/251‑1 (MQ=255)                                                                                                                                                                                                               ggAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCg                        >  1:638004/1‑251 (MQ=255)                                                                                                                                                                                                                                     gTCGCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATATCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGTTAAt  <  1:93008/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   AGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTTGTTTTTGAAACGCTCACCGCGCAACGTTTACACCCAATGCAATATTTGCTGGTGGGGCTTTCATTGGTGATGTTTTATTTGCTCTTGCTGGCGCTTTCTGAACATACCGGTTTTACCGTGGCATGGATAATCGCCAGTCTGATTGGGGCGATAAT  >  NC_000913/4638787‑4639265

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