Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,240	G→A	L21L (TTG→TTA)	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,638,240	0	G	A	98.3%	195.1 / ‑3.9	59	L21L (TTG→TTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (26/32);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

AAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTG  >  NC_000913/4637994‑4638477                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              aaGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTc                                                                                                                                                                                                                                           >  1:353087/1‑251 (MQ=255)  agcagcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGACCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTcc                                                                                                                                                                                                                                          >  1:439384/1‑251 (MQ=255)   gccgtccgcTGGGGTTGGTGTTCGGCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAAtt                                                                                                                                                                                                                                            <  2:486656/240‑1 (MQ=255)           cTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAAtt                                                                                                                                                                                                                                            >  1:486656/1‑240 (MQ=255)             ggggTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTGAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCGCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTGGGTGCATTATTGCTGTTGTTAATTCCGATAATGc                                                                                                                                                                                                                                  <  2:66631/248‑1 (MQ=255)                     cGTTCGTCATTGAGGTCGCCCGTTTGTTTAACGTCGAAGTCACTTTGCGCTACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAg                                                                                                                                                                                                                       <  2:472808/251‑1 (MQ=255)                      gTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGTAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                                                                                                                                                                                                                      >  1:192060/1‑251 (MQ=255)                                   gTCGCCCGTTTGTTTAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg                                                                                                                                                                                                         >  1:254995/1‑251 (MQ=255)                                              gtttAACGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTGCGACTTCACCGTCACGTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTAc                                                                                                                                                                                               <  1:309472/250‑1 (MQ=255)                                                  aaCGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAg                                                                                                                                                                                          <  1:471578/251‑1 (MQ=255)                                                        gAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGAtgt                                                                                                                                                                                    <  1:113350/251‑1 (MQ=255)                                                                    cgcAACGTGCAGGAAGGTGGCGTGCTGGCCGCGCTTCGACTTCACCGGCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTTCCGTAGCGATGTGGAAGATGCGAt                                                                                                                                                                        <  1:289961/251‑1 (MQ=255)                                                                     gcAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGAtt                                                                                                                                                                       <  2:353087/251‑1 (MQ=255)                                                                          gtgaaGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGGGATTGTCGAACGTGCTGATTACCGTCACGATGTGGGAGATGCGGTTCGccc                                                                                                                                                                  >  2:488878/5‑250 (MQ=255)                                                                                         gagCTGGCCTACCTTCGACTTCCCCGTCCCTTCAAATAGCTTCACATTCTTCCCACATAGTCTTCGTCTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAACTCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTTCGATAATGCTGATTCGGCTGGTGATTGTCGAACGTGCTGGTTACCGAAGCGATGTGGAAGATGCGATTCGCCCAAGtcaccgcggacc                                                                                                                                                   >  2:52229/3‑240 (MQ=255)                                                                                         gTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACAGTAGCGGTGTGGAAGATGCGATTCGCCAAAGGACCAGCGGGcc                                                                                                                                                   >  2:529978/1‑251 (MQ=255)                                                                                              ggCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCaa                                                                                                                                               <  1:140039/250‑1 (MQ=255)                                                                                                    gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCg                                                                                                                                        >  1:459561/1‑251 (MQ=255)                                                                                                        ccaCTTTACCGTCACTTCACATAGCTTCAAATTCTTCCCACTTATTCTTCGTATCCTTCTGCCATTGCAAATGAGAAGACTATGTTGAAAGCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGGTGTTACTTCCGTTAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCATCGGGCCGCAAAAACTCGTTgg                                                                                                                                    <  2:329577/249‑1 (MQ=255)                                                                                                                    cacttcacttAGCTTCAAATTCTGCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGGAGACTATGTTGAAATCCCCCCGGTTCTGGAAATTGACTAGCCTGTTTGGTTCAGTATTGCTGTTGTTAATTCCGTGAATGCTGATTCGGCAGTTGATTGTCGTACGTGCTGATTACCGTAGCGATGGTGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATc                                                                                                                          <  1:537366/249‑1 (MQ=255)                                                                                                                     aCTTCTCATTGCTTCAAACTCTTCCCACATAGGCGTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCGGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCt                                                                                                                       <  2:420348/251‑1 (MQ=255)                                                                                                                      cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTAGTTGGGCCGCTCATCGCt                                                                                                                       >  2:58432/1‑250 (MQ=255)                                                                                                                         cacaTAGCTGCAAATCCTGCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCGCTGTTCTGGAAAATGAGTAGCCTGTTTGGTGCAGTATGGCTGTTGTTAATTCCGATAATGCTGAGTCGGCAGGTGATTGGCGAACGTGCTGATGACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACTAGCGGGCCGCAAAAACTCGGTGGGCCGCTCATCGCTATTc                                                                                                                   <  2:46048/251‑1 (MQ=255)                                                                                                                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCt                                                                                                                 <  1:182351/251‑1 (MQ=255)                                                                                                                            aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCtg                                                                                                                >  1:484401/1‑251 (MQ=255)                                                                                                                                  tcaaatatttcCCACATATTCTTCGTATCCTGCTGCCAGTGCAAAGGTGAAGTCAATGTTAAAATCCCCCCTGTTCTGGAAAATGACGATCCTGTTTGGTGCAGTATTTCTGTTGTTAATTCCTATAATGCTGATTCGGCAGGTGATTTGCGAACGTGCTGATTACCGTAGCTATTTGGAAGAGGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTGTTCCTGTGAcc                                                                                                           <  2:513258/242‑1 (MQ=255)                                                                                                                                          ttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTtat                                                                                                  >  1:353516/1‑251 (MQ=255)                                                                                                                                           tcctacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGTGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCttt                                                                                                    <  1:273590/244‑1 (MQ=255)                                                                                                                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTACTTCCGATAATGCTGCTTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACg                                                                                               >  1:523907/1‑251 (MQ=255)                                                                                                                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGCGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACg                                                                                               >  1:426128/1‑251 (MQ=255)                                                                                                                                                         cGTATCCTGCTGCCATTGCAAAGGAGAAAACTATGTTAAAATCCCCCCTGTTCTGGAAAATAACTAGCCTGTTTGATGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTCCCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCACAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGATCTTTATACGGTGCAGGAag                                                                                     >  1:404096/1‑249 (MQ=255)                                                                                                                                                              ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATa                                                                               >  1:246420/1‑250 (MQ=255)                                                                                                                                                                 gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAAc                                                                           >  1:99378/1‑251 (MQ=255)                                                                                                                                                                        ttGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCACTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTgg                                                                      >  1:110584/1‑249 (MQ=255)                                                                                                                                                                         tGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGc                                                                   >  1:94356/1‑251 (MQ=255)                                                                                                                                                                                   aagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGa                                                           >  1:408565/1‑249 (MQ=255)                                                                                                                                                                                     gattttgttatatctcccccTGTTCTGGCAAATGACTAGCCTGTTTGGTTTAGTATGGCTTTTGTTAATTCCGATAAGTCTGACTCGGCAGTTGAGTTCCGAACGGGCTTATTACCGTAGCGATTTGGAAGATTCGATTCTCCAAAGTGTCAGCGGGCTGCAAAAACTCGTTGTGCCGCGCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAATCGAAGtt                                                       <  2:448463/236‑1 (MQ=255)                                                                                                                                                                                     gaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGtt                                                       <  1:111585/251‑1 (MQ=255)                                                                                                                                                                                           gTTGAAATCCCCCCTTTTCTTGAAAATGACTATCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATc                                                  <  2:254995/250‑1 (MQ=255)                                                                                                                                                                                              gAAATCCCCCCTGTTCTGGAAAATGCCTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCAtt                                              >  1:399886/1‑251 (MQ=255)                                                                                                                                                                                                     ccccTGTTCTGGAAAATGACTAGCCTGGTTGGGGCAGTATGGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGtt                                       <  2:176518/251‑1 (MQ=255)                                                                                                                                                                                                         tGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGc                                                                                                                                                    <  2:40024/138‑1 (MQ=255)                                                                                                                                                                                                         tGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGc                                                                                                                                                    >  1:40024/1‑138 (MQ=255)                                                                                                                                                                                                          gTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCt                                   <  1:479194/250‑1 (MQ=255)                                                                                                                                                                                                            tCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGa                                 <  1:488878/250‑1 (MQ=255)                                                                                                                                                                                                            tCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAAATCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTg                                  >  1:554263/1‑249 (MQ=255)                                                                                                                                                                                                              tGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTc                              >  1:446744/1‑251 (MQ=255)                                                                                                                                                                                                               ggAAAATGACTAGCCTGTTTGGTTCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAg                                <  1:330382/248‑1 (MQ=255)                                                                                                                                                                                                               ggAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCt                                                                                                                 >  1:93178/1‑167 (MQ=255)                                                                                                                                                                                                               ggAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAg                                <  1:571309/248‑1 (MQ=255)                                                                                                                                                                                                               ggAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATCCGCCAAAGTACCAGCGGGCCGCAAAAACTCGCTGGGCCGCTCATCGCTATTCCt                                                                                                                 <  2:93178/167‑1 (MQ=255)                                                                                                                                                                                                                 aaaaTGACTAGCCGGTTTGGGGCAGTATTGCTGTTGGTAATTCCGATAATGCGGATTCGGCAGGTGATTGTCGAACGTGCCGAGGACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCGGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAtt                           <  2:399886/251‑1 (MQ=255)                                                                                                                                                                                                                    aTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgat                        <  1:169273/251‑1 (MQ=255)                                                                                                                                                                                                                     tGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatg                       >  1:274765/1‑251 (MQ=255)                                                                                                                                                                                                                         tAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATa                                                                               <  1:47276/191‑1 (MQ=255)                                                                                                                                                                                                                             cTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggttgatgg               >  1:239245/1‑251 (MQ=255)                                                                                                                                                                                                                               gTTTGGTGCAGGATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggttgatgg               <  2:239245/249‑1 (MQ=255)                                                                                                                                                                                                                                           aTTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCAGCGCGATTCCTGTGACCGAGCTTTATACg                                                                                               <  1:144366/157‑1 (MQ=255)                                                                                                                                                                                                                                           aTTGCTGTTGTTAATTACGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTg  <  1:111416/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              AAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTG  >  NC_000913/4637994‑4638477

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