Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,114,389	C→T	A368T (GCA→ACA)	mdtG ←	efflux pump MdtG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,114,389	0	C	T	98.1%	151.2 / ‑4.6	52	A368T (GCA→ACA)	mdtG	efflux pump MdtG
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/20);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

GCGTCAATGGCTTCTTCAAGCGTTGCGTACATGGTCATCACGATTCTCCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGATCAACAGCAGTACAGAAAAGATCAGCGCTGTAATCAGGATCTTTTCGG  >  NC_000913/1114144‑1114631                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     gCGTCAATGGCTTCTTCAAGCGTTGCGTACATGGTCATCACGATTCTCCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCa                                                                                                                                                                                                                                                 <  1:194792/251‑1 (MQ=255)      aaTGGCTTCTTCAAGCGTTGCGTACATGGTCATCACGATTCTCCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAAt                                                                                                                                                                                                                                            >  1:206476/1‑251 (MQ=255)                  aaGCGTTGCGTACATGGTAATCACGATTCTCCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGCAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCAAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACg                                                                                                                                                                                                                                >  1:61179/1‑251 (MQ=255)                          cGTACATGGTCATCACGATTCTCCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTCGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAat                                                                                                                                                                                                                        >  1:171282/1‑251 (MQ=255)                                       caCGATTCTCCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATc                                                                                                                                                                                                                     >  1:8355/1‑241 (MQ=255)                                       caCGATTCTCCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATc                                                                                                                                                                                                                     <  2:8355/241‑1 (MQ=255)                                            ttCTACCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTGCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGGTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGt                                                                                                                                                                                                       <  2:205025/250‑1 (MQ=255)                                             tctcCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGtt                                                                                                                                                                                                      >  1:147812/1‑250 (MQ=255)                                             tctcCCGTTGAGTTGAGGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTTCAAGTATGAAAGGCGAAAAATCAGGTCGATATCTGTGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGGCAGCGGTGACGAGAAATACCGCTCTGAAACCTTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGtt                                                                                                                                                                                                      <  1:480669/250‑1 (MQ=255)                                                 ccTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCGCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTTAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGc                                                                                                                                                                                                 <  2:429814/251‑1 (MQ=255)                                                        gTTGATGAGGGTTCAGGGAAAAGAATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGCAACCGTAGTTAGCTGAAATCGCTGTTCCAATCAATGGTCAGGTACCGTTGCCAATATCACGAACCGCTTGGtaaaagcgga                                                                                                                                                                                              >  2:419125/1‑238 (MQ=255)                                                            atgaGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGATTCCAATATCACGAAACGATTGGTtatcactgaatataagc                                                                                                                                                                                      >  2:354860/1‑237 (MQ=255)                                                              gaGGTTTCAGGGAAAAGGATAGCTGATTCTCAGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGcc                                                                                                                                                                                     <  1:271034/250‑1 (MQ=255)                                                                               gATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGAATGCGTTGAATAACACTACGCCAGCGGTGCCGCGAACTACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGCAACGATTGGTTATAACTGAAGATACGCCCGGCGATCTGGTTGCTCg                                                                                                                                                                   >  2:260799/1‑251 (MQ=255)                                                                                   gCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGtt                                                                                                                                                               >  1:413660/1‑251 (MQ=255)                                                                                                         tATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTAcg                                                                                                                                         >  2:466295/1‑251 (MQ=255)                                                                                                             aaaGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTAcggcgg                                                                                                                                     >  1:90926/1‑251 (MQ=255)                                                                                                               aGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTCCGGCgggg                                                                                                                                   >  2:446527/1‑251 (MQ=255)                                                                                                                  cGAAACATCAGGTCGATACCTGGGGCATTCGACGACGACGGAGACTGTTCCATGAATAGACTGCGGTGACTAACACTACGCCAGCGGCGACGAGAAATACCGCTATGAAACCGCAGATCGCTGAAATAGCTGTTCCCATCAATGGTCCGGCAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGGTACGCACGGAGTTCTGGATGCTAGAGTTGGAAACCACCAGTGTCTGGACGGCGGGGatg                                                                                                                                >  2:259639/1‑249 (MQ=255)                                                                                                                      aaaccAGTTCGATACCTGGGGTAGTCGACGACGACGTAGACTGCTCCATGAATAGACTGCGTTGAAGAACACTCCGCTAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTACCGGTGCCAATATCACGAACCGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGTTCGAGTGGTAAACCAACAGTGTCTGTACGCCGGGGAGTAGTg                                                                                                                            <  2:12354/247‑1 (MQ=255)                                                                                                                          cAGTTCGATACCTGGGGTATTCCTCAACGACGTAGCCTGTTCCATGAATCGACTGCGTTGAATAACTCTACGCCAGCGGTGACGAGAAATTCCGCTCTGAAACCGTTGTTCGCGGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCTCGAAACGATTGGTTATAGCTGAAGATACGCCCGGTGATCTGTTTGCTCGAGTTGTAAACCACCAGTGTCTGTACGGCGGGGAGTAGTGCAc                                                                                                                         <  2:387908/250‑1 (MQ=255)                                                                                                                             ttCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGAGACGAGAAATACCGCTCTGAAACCGAAGTTCGCTGAACTCGCTGTTCCCATCAGTGGTACGGTAACGTTGCCAATATAACGAAACGATTGGTTTTAGCTGAAGATACGCCCGGCGAGCTGGTTGCTCAGGTTGTAAAACAACAGTGTCTGTACGGCGGGGaggggtgccacgt                                                                                                                      >  2:426280/1‑239 (MQ=255)                                                                                                                                   aCCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGTGGTGAGTAGTGCACCATCggcgg‑c                                                                                                              >  1:118491/1‑251 (MQ=255)                                                                                                                                        gggTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGa                                                                                                         <  1:168376/251‑1 (MQ=255)                                                                                                                                                   gacgacgTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGGGTAGGGCACCATCGGCGG‑CACCGAGCAAAAAACGt                                                                                              >  1:255745/1‑251 (MQ=255)                                                                                                                                                       acgTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAACCGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTACACCAACAGTGTCTGTACGGCGGGGGGGAGTGCACCCTCGGCGG‑GCGCGGGCAAAAAACGCaaa                                                                                           >  2:10573/1‑250 (MQ=255)                                                                                                                                                       acgTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTaaaa                                                                                          <  1:254538/251‑1 (MQ=255)                                                                                                                                                                        aaTAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑Gc                                                                        <  2:29556/251‑1 (MQ=255)                                                                                                                                                                                 cGTTGAATAACACTACGCCAGCGGTGACGAGAAATCCCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCCTCGGCCG‑CCCCGAGCAAAAAACGTAAAATCCCAAATTTCAATG‑GGGTccggccg                                                               >  2:201876/1‑245 (MQ=255)                                                                                                                                                                                    tGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGATGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGGAGTGGCCCATCGGCGGGCACCG‑GCAAAAAACGGAAAAAGACAAGTTGCAATG‑GCGGCTGAAAGTa                                                             >  2:339438/1‑250 (MQ=255)                                                                                                                                                                                         aaCACTACGCCAGCGGTGTCGAGAAATACCGCTCTGAAACCGTAGTTTGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACTGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACa                                                        <  2:275293/250‑1 (MQ=255)                                                                                                                                                                                         aaCACGACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCTTAGTTCGCTGAAATCGCTGTTCCCATCAATGGGCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCGCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTATTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACAt                                                       <  2:235773/251‑1 (MQ=255)                                                                                                                                                                                          acccTACGCCAGCGGTGACGTGAAATCCCTCTCTGAATCCGTTGTTCGCTTAAATCGCTGTTTCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACAtt                                                      <  2:91706/248‑1 (MQ=255)                                                                                                                                                                                                 gCCAGCGGTGACGAGAAATACCGCCCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGATCa                                               <  2:147812/251‑1 (MQ=255)                                                                                                                                                                                                  ccAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGAAAACTCCCACGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGATCa                                               >  1:17022/1‑250 (MQ=255)                                                                                                                                                                                                       ggTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGGGTAGTGCACCATCGGCGG‑CACCGCGCAAACAACGTAAAATCCCAAGATGCAATG‑GCGTCTGAACGTTAGAAAATGGGATCaaacgcc                                         >  2:117762/1‑246 (MQ=255)                                                                                                                                                                                                       ggTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCCAGCAAAAAACATAAAATCCCAAGATGCAAAG‑GCGGCCGAAACTAAGAAATTggggcacaaagc                                          >  2:28152/1‑241 (MQ=255)                                                                                                                                                                                                             gagaAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGGGTAGTGCACCCTCGGCGG‑CACCGAGCACAAAAAGTAAAATCCCAAGATGCAATG‑GCGTCTGAACGTAAGACATTGGGATCAACAACAGTaa                                     >  2:445411/1‑248 (MQ=255)                                                                                                                                                                                                              agaAATACCGGTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCCATGGTCCGGAAACGTTGACAATATCACGAAACGATTGGTTATAGCTGAAGGTACGACCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGGAGGGCAACCTCGGCGG‑CACAGAGCAAAAACCGTACAATCACATGATGCAAAGTGC‑ACTGAACGAAAAACATTGGGAGCAACAGCAGTAc                                     >  2:48335/1‑248 (MQ=255)                                                                                                                                                                                                              agaAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGAAATTGGGATCAACAGCAGTAc                                     >  1:355531/1‑248 (MQ=255)                                                                                                                                                                                                                     ccGCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGGTCAAACGCAGTACAGaaa                                >  1:213497/1‑246 (MQ=255)                                                                                                                                                                                                                       gCTCTGAAACCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGATCAACAGCAGTACAGAAAAGATc                           >  1:446050/1‑249 (MQ=255)                                                                                                                                                                                                                            gAAACCGTTGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGAACTGATTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCAG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAATACACTGGGATCAACAGCAGTACAGAAAAGATCAgcgcgg                    <  2:253252/251‑3 (MQ=255)                                                                                                                                                                                                                             aaaCCGTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCACAAATTGCAATG‑GCGTCTGAAAGTAAGGCATTGGGATCGAAAGCCGTACAGAAAAGATCAGCGCTGt                   >  1:189210/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATc                                                                                                                     >  1:59716/1‑150 (MQ=255)                                                                                                                                                                                                                                  gTAGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATc                                                                                                                     <  2:59716/150‑1 (MQ=255)                                                                                                                                                                                                                                    aGTTCGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGATCAACAGCAGTACAGAAAAGATCAGCGCTGTAATCAgg            <  1:445411/251‑1 (MQ=255)                                                                                                                                                                                                                                    aGTTCGCTGAAATCGCTGTTCCCATCAAGGGTCCGGTAACGTTGCCAATATCACGATACGATTGGTTATAGCTGAAGATACGGCCGGCGATCTGGTTGTTCGAGTTGTAAACCATCAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAATCGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGGCATTGGGAGCAACAGCAGTACAGAAAAGATCATCGCAGTAATCAg             <  1:124342/250‑1 (MQ=255)                                                                                                                                                                                                                                        cGCTGAAATCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGGAAGACATTGGGGTCAACAGCAATACAGAAAAAGTCAGCGCTGTAATCAGGATCt        >  1:404067/1‑251 (MQ=255)                                                                                                                                                                                                                                            gAAATCGCTGTTACCATCAATTGGCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGATCAACAGCAGTACAGAAAAGATCAGCGCTGTAATCAGGATCttt      <  2:187131/249‑1 (MQ=255)                                                                                                                                                                                                                                             aaaTCGCTGTTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAAAGTAAAATCCCAAGTTGCAATG‑GCGTATGAACGGAAGACATTGGGATAAACAACAGTACAGAAAAAATCAGCGCTGTAATCAGGATAtttt     >  2:429078/1‑249 (MQ=255)                                                                                                                                                                                                                                              aaTCGCTGTTTCCATCAATGGGCCGGTGACGATGCCAATATCACGAAACGAGTGGTTATAGCTGAAGATACGCCCGGCGGTCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCAGCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGATCAACAGCAGTACAGAAAAGATCAGCGCTGTAATCAGGAGCTTTTCgg  <  2:337920/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     GCGTCAATGGCTTCTTCAAGCGTTGCGTACATGGTCATCACGATTCTCCCTTTGAGTTGATGAGGTTTCAGGGAAAAGGATAGCTGATTCTCCGCTTTTGCAAGTATGAAAGGCGAAAAATCAGTTCGATACCTGGGGTATTCGACGACGACGTAGACTGTTCCATGAATAGACTGCGTTGAATAACACTACGCCAGCGGTGACGAGAAATACCGCTCTGAAACCGTAGTTCGCTGAAATCGCTGCTCCCATCAATGGTCCGGTAACGTTGCCAATATCACGAAACGATTGGTTATAGCTGAAGATACGCCCGGCGATCTGGTTGCTCGAGTTGTAAACCAACAGTGTCTGTACGGCGGGGAGTAGTGCACCATCGGCGG‑CACCGAGCAAAAAACGTAAAATCCCAAGTTGCAATG‑GCGTCTGAACGTAAGACATTGGGATCAACAGCAGTACAGAAAAGATCAGCGCTGTAATCAGGATCTTTTCGG  >  NC_000913/1114144‑1114631

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