Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,660,610	0	A	G	36.5%	40.9 / 51.4	53	E66G (GAG→GGG)	mdtF	multidrug efflux pump RND permease MdtF
Reads supporting (aligned to +/- strand):  ref base A (12/21);  new base G (8/11);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAAACGCCAGCACCGAATCGAAACAATAACGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTATTTCTGATAACGGCAGCCTC  >  NC_000913/3660361‑3660860                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           aaaaCGCCAGCACCGAATCGAAACAATAACGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAg                                                                                                                                                                                                                                                           <  1:341959/251‑1 (MQ=255)       ccAGCACCGAATCGAAACAATAACGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAat                                                                                                                                                                                                                                                     <  1:196233/251‑1 (MQ=255)        cAGCACCGAATCGAAACAATAACGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAat                                                                                                                                                                                                                                                     >  1:31731/1‑250 (MQ=255)            aCCGAATCGAAACAATAACGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGaa                                                                                                                                                                                                                                                <  1:246565/251‑1 (MQ=255)              cGAATCGAAACAATACCGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATg                                                                                                                                                                                                                                              <  1:312748/251‑1 (MQ=255)                 aTCGAAACAATAACGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGc                                                                                                                                                                                                                                           >  2:65848/1‑251 (MQ=255)                         aaTAACGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACGCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGc                                                                                                                                                                                                                                   <  2:138802/251‑1 (MQ=255)                          aTAACGTTGCAGGCTTAAGGGGACTTTCATGGCTACCTATTATACTGATCGCCCGGCTTTTGCCTGGGTACTTCCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGCTTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGcc                                                                                                                                                                                                                                  >  1:217640/1‑251 (MQ=255)                            aaCGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTg                                                                                                                                                                                                                                >  1:186072/1‑251 (MQ=255)                                         tAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTc                                                                                                                                                                                                                     <  2:253405/249‑1 (MQ=255)                                          aaGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTc                                                                                                                                                                                                                  >  2:287981/1‑251 (MQ=255)                                          aaGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGATTTCTGCCTGGGAACTTGCCATTATTATGATGAATGCCGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATACGCAGATTGCGACACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGATAGAAGCCTCGGTCACTCAGGTGATTGAGCAAAATATGCATGGGCTTGATGGCCTGATGTCCATGTCTTc                                                                                                                                                                                                                  >  1:184271/1‑251 (MQ=255)                                                aCTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGCCCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAAcc                                                                                                                                                                                                              >  2:133613/1‑249 (MQ=255)                                                     cATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGTTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCTCAAACTGTAGAAGACTCTGTCACTCAGGTGATTGAGCAAAATATGAATGTGCTGGATGGCCTGATGTACATGTCTTCAACCAGTGAt                                                                                                                                                                                                        <  1:27864/250‑1 (MQ=255)                                                                   ttATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCACTGcc                                                                                                                                                                                         >  1:134177/1‑251 (MQ=255)                                                                                  tttttGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGCCGATCATGAACTTACCGGTTGCACAGTATCCGCAGATAGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGCCACTCAGGTGATTGGGCCAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGAcc                                                                                                                                                                          >  1:242950/1‑251 (MQ=255)                                                                                      tGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTc                                                                                                                                                                       <  2:210179/250‑1 (MQ=255)                                                                                       gCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCg                                                                                                                                                                      <  2:129687/250‑1 (MQ=255)                                                                                                  tGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGAc                                                                                                                                                            <  2:362624/249‑1 (MQ=255)                                                                                                  tGCCAATAATATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTGCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGAc                                                                                                                                                            <  2:362634/249‑1 (MQ=255)                                                                                                                       aGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTc                                                                                                                                     <  2:79596/251‑1 (MQ=255)                                                                                                                       aGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTc                                                                                                                                     <  2:99093/251‑1 (MQ=255)                                                                                                                       aGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTc                                                                                                                                     <  1:351263/251‑1 (MQ=255)                                                                                                                         gtggtCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCaa                                                                                                                                   <  1:65848/251‑1 (MQ=255)                                                                                                                         gtggtCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACGCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACAGCTCCTGATATCGCACAGGTTCaa                                                                                                                                   <  2:33198/251‑1 (MQ=255)                                                                                                                            gtCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTg                                                                                                                                <  1:266869/251‑1 (MQ=255)                                                                                                                            gtCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACGCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTg                                                                                                                                <  2:205104/251‑1 (MQ=255)                                                                                                                                      cATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATaaa                                                                                                                       <  1:95663/250‑1 (MQ=255)                                                                                                                                       aTGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACt                                                                                                                     <  1:168984/251‑1 (MQ=255)                                                                                                                                          aaCTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCa                                                                                                                  <  2:31731/251‑1 (MQ=255)                                                                                                                                                ccGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGc                                                                                                            <  2:66837/251‑1 (MQ=255)                                                                                                                                                           gTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCGGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCAt                                                                                                 <  2:212349/251‑1 (MQ=255)                                                                                                                                                                cGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCt                                                                                            <  1:353321/251‑1 (MQ=255)                                                                                                                                                                          gcCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTgc                                                                                  <  2:277087/251‑1 (MQ=255)                                                                                                                                                                             aCCGACCATTACCGTCAGCGCTAACTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCCCTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTgcagc                                                                               >  2:59017/1‑251 (MQ=255)                                                                                                                                                                                 aCCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTgcagcagc                                                                            >  1:163801/1‑250 (MQ=255)                                                                                                                                                                                   cATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACCGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAgg                                                                         >  2:53667/1‑251 (MQ=255)                                                                                                                                                                                    aTTACCGTCAGCGATACCTATCCAGGTGCCGCTGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAggg                                                                        >  1:296967/1‑251 (MQ=255)                                                                                                                                                                                                       agccAGGTGCCGATGCGCATACGGTAGAAGACGCGGTCACTCAGGTGATTGGGCAAAATGTGAATGGGCTTGATGCCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTTTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAATCTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGTGGATTATCGTCGATAAGt                                                       <  2:62721/247‑1 (MQ=255)                                                                                                                                                                                                       aTCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTTCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGTATTAGCGTCGATAAGTCg                                                     <  1:134919/251‑1 (MQ=255)                                                                                                                                                                                                       aTCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTc                                                      <  1:287981/250‑1 (MQ=255)                                                                                                                                                                                                         ccAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGa                                                    <  1:52460/250‑1 (MQ=255)                                                                                                                                                                                                           aGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTAACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGGTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCa                                                 <  2:262374/251‑1 (MQ=255)                                                                                                                                                                                                                       gcAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTg                                      >  1:165475/1‑250 (MQ=255)                                                                                                                                                                                                                            cGGTAGAAGACTCGGTTAATCAGGTGATTGTGCAAAAGATGAATGGGATTGATGGCCTGATGTACATGTCTTAAACCAGTGATGCGGCGGGCAATGCGTCTAGCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTGCAAGTGTAAAATAAACAGCAGCTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGCCGAGCAGGAATATCCTGATGGtn                                <  2:275228/251‑2 (MQ=255)                                                                                                                                                                                                                            cGGTAGAAGACTCGGGCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGCTGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATTTCCTGATATCGCACAGGTATACGTGCAAAATAAACAGCAACTCGGCATGCCTTCATTAACTGAATCAGTGCAGCAGCAGGGGCTTAACGTCGAGAAGGCGGGCAGAAATATTCTGATGGt                                 >  2:273127/1‑250 (MQ=255)                                                                                                                                                                                                                             ggTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCGATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAg                               <  2:186072/251‑1 (MQ=255)                                                                                                                                                                                                                                  aagaCTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAgcggcg                          >  1:241106/1‑251 (MQ=255)                                                                                                                                                                                                                                   agaCTCGGTCACTCAGGTGATTGGGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTTCTGATATCGCACCGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGGAGCGGCGt                         >  2:256168/1‑251 (MQ=255)                                                                                                                                                                                                                                        cGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGATACTGGGACATCTCCTGATATCGCACAGGTTCAAGttca                                                                                                                              <  1:282126/145‑4 (MQ=255)                                                                                                                                                                                                                                        cGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGcc                                                                                                                              >  2:282126/1‑144 (MQ=255)                                                                                                                                                                                                                                           tCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTATTTCt                 >  1:238989/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTATTTCTGATAACgg         >  1:266902/1‑251 (MQ=255)                                                                                                                                                                                                                                                          aGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTATTTCTGATAACGGCAGCCTc  <  2:143949/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           AAAACGCCAGCACCGAATCGAAACAATAACGTTGCAGGCTTAAGGGGACTTTCATGGCTAACTATTTTATTGATCGCCCGGTTTTTGCCTGGGTACTTGCCATTATTATGATGCTTGCAGGTGGTCTGGCGATCATGAACTTACCGGTTGCGCAGTATCCGCAGATTGCGCCACCGACCATTACCGTCAGCGCTACCTATCCAGGTGCCGATGCGCAAACGGTAGAAGACTCGGTCACTCAGGTGATTGAGCAAAATATGAATGGGCTTGATGGCCTGATGTACATGTCTTCAACCAGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAATAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTATTTCTGATAACGGCAGCCTC  >  NC_000913/3660361‑3660860

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