Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,568,566	T→A	17.2%	intergenic (‑46/‑433)	mdtM ← / → yjiP	multidrug efflux system protein/pseudogene, transposase_31 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,568,566	0	T	A	17.2%	197.0 / 2.7	68	intergenic (‑46/‑433)	mdtM/yjiP	multidrug efflux system protein/pseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base T (7/48);  new base A (11/2);  total (18/50)
Fisher's exact test for biased strand distribution p-value = 1.27e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATACCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTAAAATTTGCGG  >  NC_000913/4568335‑4568794                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       cATACCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGACAGACTGACATCGGCATTAAAATCACGTTCCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATa                                                                                                                                                                                                                                  <  2:433676/238‑1 (MQ=255)   tACCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGAGGCAGAACAAACAACACCTTGTCCCGAGAACACCCACCAGGAAAACCAAAATTAATAca                                                                                                                                                                                                                                >  2:409528/1‑237 (MQ=255)   tACCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGCGGCATAACAAACAAATCATTGTCCGGAGAACCCCAACGGGGGAAACGAAATTCAATACCg                                                                                                                                                                                                                               >  1:12181/1‑239 (MQ=255)   tACCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGAGGCGTAACAAAAAAACCCTTTTCCGGGGAACCCCCACCGGGGAAAAGAAAATTAATaac                                                                                                                                                                                                                                >  2:173406/1‑236 (MQ=255)    aCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGt                                                                                                                                                                                                                              <  1:25397/239‑1 (MQ=255)         cAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTa                                                                                                                                                                                                                         <  2:506715/239‑1 (MQ=255)         cAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAATGGCCTAACAAACAACACCCTTTTCGGAGCACCCCCCACGGGAAAAAGAAATTCAATaaactcggga                                                                                                                                                                                                                         >  2:265606/1‑230 (MQ=255)         cAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACATACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTa                                                                                                                                                                                                                         <  1:177504/239‑1 (MQ=255)         cAGCAAGATAGAGACTGACAGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTa                                                                                                                                                                                                                         <  1:209116/239‑1 (MQ=255)           gCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGGATAACAAAAAACCCCTTTTCCGGAGAACCACAACAGGGAAAACGAAAATTAAAACCCTCGGGAAt                                                                                                                                                                                                                       >  2:240039/1‑239 (MQ=255)             aaGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTaa                                                                                                                                                                                                                        <  1:48949/236‑1 (MQ=255)                aTAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATGCCGTCGGTAATGACGa                                                                                                                                                                                                                  <  1:4614/239‑1 (MQ=255)                           cGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTga                                                                                                                                                                                                       <  1:48196/239‑1 (MQ=255)                                  agGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAg                                                                                                                                                                                                   <  2:209053/236‑1 (MQ=255)                                         aGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAg                                                                                                                                                                                                   <  1:404963/229‑1 (MQ=255)                                            cTGACATCGGCATTAAAATCACGTACCACATTAATTATCCCAGGCGCGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAGCTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc                                                                                                                                                                                        <  1:1415/237‑1 (MQ=255)                                               aCAACGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCaa                                                                                                                                                                                   <  1:174147/239‑1 (MQ=255)                                                cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCaa                                                                                                                                                                                   <  2:36923/238‑1 (MQ=255)                                                   cGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAAAAAACCCCTGTCCGGAGAACCCCCAAGGGGGAAACGAAAATCAATACCGCCGGTAATAACGAAACACAGGGGGGAAGGGCAACgaagcccgtctctt                                                                                                                                                                               >  2:404963/1‑227 (MQ=255)                                                    ggCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGGGGCATAACAAACAACTCCTTGTCCGGAGAACACCAAAGGGGGAAACGAAATTCAATACCGTCGGTAATGACGAAAAGCAGGTGAGAAGGGCAACGAGACTGCAAATGTa                                                                                                                                                                              >  2:77478/1‑239 (MQ=255)                                                     tttatataagagacagccACATTAAAGATCCCAAGCTGGAACAGATACGGCGACAGATACGCAGCAAAGTCAAACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTAc                                                                                                                                                                             <  1:468399/223‑1 (MQ=255)                                                                aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATg                                                                                                                                                                 <  1:173406/239‑1 (MQ=255)                                                                  gTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTa                                                                                                                                                               <  2:88454/239‑1 (MQ=255)                                                                  gTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGGCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTa                                                                                                                                                               <  1:9793/239‑1 (MQ=255)                                                                      cacaTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCa                                                                                                                                                           <  2:415261/239‑1 (MQ=255)                                                                      cacaTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCa                                                                                                                                                           <  1:429030/239‑1 (MQ=255)                                                                            aaagATCCCAGGCTGGATCAGAACCGTCGACAGATACGCAGCAAAGTCATACAAAAGCAACGCCAGCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCGTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑Ccaca                                                                                                                                                    <  1:521417/236‑1 (MQ=255)                                                                             aTGATCTCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCacac                                                                                                                                                   <  2:258504/239‑1 (MQ=255)                                                                             aTGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTa                                                                                                                                                               <  1:257569/228‑1 (MQ=255)                                                                             aTGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCacac                                                                                                                                                   <  2:89961/239‑1 (MQ=255)                                                                               gATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAACCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTg                                                                                                                                                 <  2:119019/239‑1 (MQ=255)                                                                                aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACTCCTTGTCCGGGGGACCCCAACAGGGGAAACGAAAATCAAAACCGTCGGTAATAAAGAAAAGAAGGGGGGAAGGGCAACGAGAATGCAAATGTACTGAATTT‑TACGTAATTCA‑acacacag                                                                                                                                                 >  2:37303/1‑236 (MQ=255)                                                                                 acccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGTGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGGCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGcc                                                                                                                                               <  2:440597/238‑1 (MQ=255)                                                                                  cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGCTGCATAACAAACAACTCCTTGTCCGGGGAACCCCAACGGGGAAAACGAAAATCAATACCGTCGGTAATGACGAAAGGCCGGGGGGGAGGGGACAGAGAACGCAAAAGTAATGAATTTACATG‑AATCAA‑CCACACCGcc                                                                                                                                               >  2:81299/1‑238 (MQ=255)                                                                                   ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGCGGGATAACAAACAAACCCTTTTCCGGAAAACCCCAACGGGGGAAACGAAAATCCAAAACCTCGGTCATGAGAGAAAGCAGGTGGGAAAGGGCACGAGAACGCAAAAGTTATGAAatttaacgtaacaatcacaactc                                                                                                                                                  >  2:410892/1‑216 (MQ=255)                                                                                       gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                                                                                                                                                                                     <  1:384221/181‑1 (MQ=255)                                                                                       gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAAAAACTCCTTGTCCGGGGAACCCCAACGGGGAAAACGAAAATCAATAACGTCGGTGATGACGCAAAGCAGGGGGGAAAGGCAACGAGACTGCAAAAGTACTTAAATT‑CATGTAATCAA‑Acaacctacccgtgga                                                                                                                                         >  2:184621/1‑226 (MQ=255)                                                                                       gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAACAAAAAAATCCTTGTCCGGGGAACCCCAACGGGGGAAACGAAAATCAAAACCGTCGGTAATGACGAAAAGCAGGGGAGAAGGGCAACGAGACTGCAAATGTAATGAATTT‑CATTGAATCAA‑Ccaacccacgctctta                                                                                                                                         >  2:257569/1‑226 (MQ=255)                                                                                          ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAacac                                                                                                                                      <  1:308549/239‑1 (MQ=255)                                                                                             tCAGATCCGGAGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAAcacaca                                                                                                                                   <  1:230921/239‑1 (MQ=255)                                                                                                 aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCCCAAAAAAAAACTCCCTGTCCCGGGAACCCCCACGGGGGAAAACAAAATCAAAACACGCGGGAAATACGAACAGCAGGGGAGAAAGGGAAAAAGACTGCAAAATTAATGAAATT‑TATGTAAACAACCCACACtcctttgtaacaacaaaca                                                                                                                                >  2:338551/1‑221 (MQ=255)                                                                                                  tCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACacaga                                                                                                                                 <  1:410892/236‑1 (MQ=255)                                                                                                      tCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGGCAAACAAACAAATCCTTGTCCGGGGAAACCCAACGGGGGAAAAGAAAATCAATACCGTCGGGAAAGAAGAAAAGCAGGGGGGAAGGGGAACGAGACTGCAAAAGTACTGAATTT‑CATGTAATCAA‑CCACAACGCCTGTGAACACACAGACAGaaaa                                                                                                                          >  2:62839/1‑239 (MQ=255)                                                                                                                      aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACt                                                                                                             <  2:505362/236‑1 (MQ=255)                                                                                                                        aGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCggg                                                                                                        <  1:93937/239‑1 (MQ=255)                                                                                                                           cagaCAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGGGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTc                                                                                                      <  1:240039/235‑1 (MQ=255)                                                                                                                           cATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTc                                                                                                      <  1:276428/238‑1 (MQ=255)                                                                                                                           cATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTc                                                                                                      <  1:175453/238‑1 (MQ=255)                                                                                                                           cATACAAAAACAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCa                                                                                                     <  1:338551/239‑1 (MQ=255)                                                                                                                             tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGa                                                                                                   <  1:15448/239‑1 (MQ=255)                                                                                                                             tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGa                                                                                                   <  1:455723/239‑1 (MQ=255)                                                                                                                             tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCGTGGCCGTAGAACCCCAACGGGGAAAACGAAAATCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATAT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGa                                                                                                   <  1:373057/239‑1 (MQ=255)                                                                                                                                     cAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATaa                                                                                           <  1:409528/239‑1 (MQ=255)                                                                                                                                         gCCATCGGGAAAAACAGCGTGGCGGCATGGCGTGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGAt                                                                                       <  1:215080/239‑1 (MQ=255)                                                                                                                                                  aaaaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGGCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATa                                                                               <  1:511237/238‑1 (MQ=255)                                                                                                                                                      aCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGtcat                                                                          <  2:264113/239‑1 (MQ=255)                                                                                                                                                           gTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCGTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGtcatca                                                                        <  1:37303/236‑1 (MQ=255)                                                                                                                                                                  cATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCTTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTac                                                               <  1:76801/238‑1 (MQ=255)                                                                                                                                                                   aTGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTacac                                                             <  1:348398/239‑1 (MQ=255)                                                                                                                                                                              aaaaaCGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTg                                                  <  2:429314/239‑1 (MQ=255)                                                                                                                                                                                    gTGGCATAACAAACAACTCCTTGTCCAGAGAACACCAACGGGGAAAACGAAAATCAATACCGTCGGTAATGACAAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGt                                            <  1:265606/239‑1 (MQ=255)                                                                                                                                                                                         aTAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCAt                                       <  1:62839/239‑1 (MQ=255)                                                                                                                                                                                                acaacTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATaa                                                                                           <  2:466099/180‑1 (MQ=255)                                                                                                                                                                                                acaacTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATaa                                                                                           >  1:466099/1‑180 (MQ=255)                                                                                                                                                                                                       cTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATaca                           >  1:377540/1‑237 (MQ=255)                                                                                                                                                                                                               gagaACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAgtgt                 >  1:246281/1‑239 (MQ=255)                                                                                                                                                                                                                    ccccAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTa            >  2:399121/1‑239 (MQ=255)                                                                                                                                                                                                                               aaaaCGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTAAAATTTGCgg  <  1:184621/238‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       CATACCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTT‑CATGTAATCAA‑CCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTAAAATTTGCGG  >  NC_000913/4568335‑4568794

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