Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,286	0	C	T	24.1%	68.8 / 28.7	61	intergenic (+492/+150)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base T (8/6);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

GACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCA  >  NC_000913/4295989‑4296576                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                         aaCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                                                                                                                       <  2:378987/300‑1 (MQ=21)  aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                                                                                                                                                                      >  2:872890/1‑301 (MQ=21)   cAGCGCAATACCCGGATTTCGCATGATTTTGTAAGCCGAATAACGCGTTTACGCCGGATCCGACAGCAACGCCTGATTCGACGCTTACCGCGTATTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGATATCAACGCCTGTTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                                                                                                                                                                                     <  2:138994/301‑1 (MQ=11)    aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGAATAAGGCGTTTACGCCGCAGCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                                                                                                                                                                                    <  2:359955/301‑1 (MQ=21)        cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGCTGCg                                                                                                                                                                                                                                                                                                >  2:903244/1‑301 (MQ=21)             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTGAGCGTTCTTTTGTAGGCCGGATAAGGCGTTTACGCCGAATCCGACATCCATGCCTGATTCGACGcc                                                                                                                                                                                                                                                                                           >  1:708628/1‑300 (MQ=255)             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAAAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAGGccggtagattgga                                                                                                                                                                                                                                                                                            >  2:822498/1‑289 (MQ=255)               cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACCTCAATGCCTGATGCGACACttt                                                                                                                                                                                                                                                                                         >  2:624722/1‑300 (MQ=19)                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTacc                                                                                                                                                                                                                                                                                       >  2:394848/1‑299 (MQ=21)                     tAGCGTGTGTTTGTAGGTCTGATAAGGCGTTTACGCCGCTTCCGACATCAACGCCTGATGCGCCGCTTGCCCCGTCTTATCAGGCCTACGCCAGACAGCGAAAAAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                                                                                                                                                                                                                                                                    <  1:872663/300‑1 (MQ=11)                       cccacgcgctgcccatccgacctgccGTTTACGCCGCCTCCGACATCACCGCCTGATGGGGGGCTTGCCGCTTCTTATCAGGCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTATGCCGCATCAGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                                                                                                                                                                                 <  1:478233/276‑1 (MQ=11)                                   atgTCGGATGCGGCGTAAACGCCTTATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCACACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGAATCCGACCTCACCGCCTGATGCGACGCTTAACGTGTCTTATCAGGCCTACTCCCTACCGCGTCATAGCCTGATTGAGCGTTACTTGGTAGGTCGGTCAAGGGATTTGCACAGCATCCGTCACCAATGCGTGATGCGAGGCTTGCCACGTCTTTTACGGCCtgt                                                                                                                                                                                                                                                                     >  2:98837/3‑299 (MQ=2)                                       cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAATGCCGGATGCGACGCTTGCCGCGTCTTCTCAggtatctcttt                                                                                                                                                                                                                                                                 >  1:419151/1‑292 (MQ=255)                                             aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGCCCTACGCCAGACAGCGCAATAGGCTGATTATGCGTTATTTTGTTGTCCGGTACAGGCGTTCACGCCGCAGCCGACACCAACGCTTGACGCGACTCTTGCCGGGTATTATCCGGCCtgcctttatccct                                                                                                                                                                                                                                                           >  1:350890/1‑289 (MQ=11)                                             aGGCGTTGACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                                                                                                                                                                                   <  2:4329/293‑4 (MQ=9)                                              ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCAGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtcgatcggaagggcccacgtctg                                                                                                                                                                                                                                                           >  2:478366/1‑278 (MQ=11)                                                                                                                                                                  gTTTACGCCGCTCCCGACATAAACCCCGTATTCGACGCTTGCCCCGTCTTTTCAGGCATACGCCCTAAAGCGAAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                                                                                  <  2:528486/301‑1 (MQ=11)                                                                                                                                                                          cGCATCCTACATCAACGCCTGAGGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                                                                          <  2:45680/301‑1 (MQ=255)                                                                                                                                                                                 caactcaacgcctgagtgagctctggccgcgTCTTATCATGCCTACGCCAGACAGCTCAATATCCTGAGTTTGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                                                                                   <  1:808502/274‑1 (MQ=11)                                                                                                                                                                                   gatcgaggcggggggagagggtgggcgtgTCTGATCGGGCCTAGGCGAGACAGCGCAAGAGGCTGAGTTGGGGTGGTTTTGTAGGTCGGAGAAGGCGTTTAGGCCGCATCCGACAGCAAGGTGTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                                 <  1:394735/273‑1 (MQ=14)                                                                                                                                                                                      cAACGCCTGAGGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                               <  2:419267/300‑1 (MQ=255)                                                                                                                                                                                      cAACCCCTGATGCTACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                               <  2:708823/300‑1 (MQ=255)                                                                                                                                                                                       aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCCTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCCGATGCGACGCTTGCCGCGCCCTATCAGGCCTACCTTAACCGTTGGTTAATTTTTAAAACACTATacac                                                                                                             >  1:825268/1‑301 (MQ=11)                                                                                                                                                                                        gcgcgtgatgtggatcagtccgcgTCTGATCACGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGAATCCTCCATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                            <  1:903007/281‑1 (MQ=14)                                                                                                                                                                                            cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                                                                                        >  2:922628/1‑301 (MQ=255)                                                                                                                                                                                             tgataagacgcggcaagcgTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCCGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATCAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTGGGTTAATTTTTCAAACTCTTCACACCCTGtt                                                                                                       >  1:447857/17‑300 (MQ=11)                                                                                                                                                                                              gATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCCATGCCTGATGCGACGTTTCTCGCGTCTTATCAGGCCTCTATTAACCGTTGGTTAATCTTTCACACTCTTTACAccctctat                                                                                                      >  1:1039764/1‑296 (MQ=11)                                                                                                                                                                                                tGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTGCGGACAAGGCGTTT‑‑ACGC‑C‑GCATCCGACCTCAATGCCTGATGCGACGCTTGCCGCGTTTTATGAGGCCTCACTTAACAGTGGGGTAATTTTTCaacatccttacagctggcaccc                                                                                                    >  1:292930/1‑281 (MQ=11)                                                                                                                                                                                                   gACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑CCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACACTTCACACCGGTTCTCAAAc                                                                                                 >  1:217782/1‑301 (MQ=11)                                                                                                                                                                                                          gccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCTCGTATTATCAGGCCTACGCCAGACAGCGCAAAATCCTGATTTATCGTGATTTCGTAGTTCGGATAAGGCGTTT‑‑A‑‑CGCAGCATCCGACACCCATTCCAGAGGAGACGCTCTCCGCGTCGTATCAGGCTTATCTTAATCGTTGGTTCATTTCTCAAACTCCTCCCCACCGGTATCaacctctccc                                                                                           >  1:564953/3‑292 (MQ=11)                                                                                                                                                                                                               gggttttcaagccgagtccAAACAGCGTATTTGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGCCGTTTATGCCGCATGCGACATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                                                                                     <  1:922388/284‑1 (MQ=17)                                                                                                                                                                                                                  ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGTGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGGTAATTTTTCAAGCTCTTCACACCCGGTACCAAACCCTTCCAAACAGCTc                                                                                  >  2:311230/1‑301 (MQ=255)                                                                                                                                                                                                                    aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCCAGCTc                                                                                  >  2:855289/1‑299 (MQ=255)                                                                                                                                                                                                                        ggCCTACGCCAGACAGCGCAATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTATGCAGCATACGACATCAACGCCTGATGCGAGGCTTAACGGGTATTATCAGGCCTACGGCCGCCAGCGCAACAGCCTGAGTTCTGGTGAGTTTGTATGTCGGAGAAAGAGTGT‑‑ACGC‑C‑GCATCCGACATCAATGACTGATGCTTGATTTCGCGAGTCTTCCCAGGCCTGCGTTAGCCGTCTGTTCAGTTTTCACTTTCTTCACACCGCGCCTCAAGCACTTCCATCCAgtccagcttt                                                                            >  1:576434/1‑292 (MQ=2)                                                                                                                                                                                                                        ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCGTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATTAAACCCTTCCATACACCtatgcttc                                                                            >  1:456371/1‑294 (MQ=255)                                                                                                                                                                                                                             aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa                                                                       <  2:825484/301‑1 (MQ=255)                                                                                                                                                                                                                                     cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGACTGATGCGAGGTTTGCCGCGTCTTATCAGGCCGATCTTCACCTTTGGTTAAATTTTCAAACTTTTCTATACCGGTCTCAAACCCTTCCAAAAAGCTCAGGTTAAACCacgccagcg                                                               >  1:118127/1‑293 (MQ=14)                                                                                                                                                                                                                                       gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTCTCCGTCCTCTCTTACCCGTTGTTTAATTTTTCAAACTCTTCACACCCGGTATCACACCCTTCCATACAGCTCAGGTTCAACCCGttccgcccc                                                             >  2:142790/1‑293 (MQ=17)                                                                                                                                                                                                                                          cagtcgtctggtttttcgggATTTTGTAGGGCGGATAAGGCGTTTGCGGCGAATCAGAAATAAAGGGCTGATTGGACGCTTAAAGCGTCTGATCGGGGCTACGCCAGAAAGCGCAATAGCCTGATGTAGCGGGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCttt                                                          <  1:142758/282‑1 (MQ=18)                                                                                                                                                                                                                                          cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCATTCACAACAGCTCAGGTTCATcccactcaccccattt                                                          >  1:660103/1‑287 (MQ=17)                                                                                                                                                                                                                                                ccttattcatcttgatgtgttAGGTCGGATAAGTCGTTTATGCCTCATCCGACATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTg                                                    <  2:830549/281‑1 (MQ=18)                                                                                                                                                                                                                                                   gATTTAGCGTTGTTTTGTAGGTCGGATAAGGCGTGTAAGCCGGAACCGACATAAACGCATGATGCGACGCTTAACGCGTCTTATGAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGt                                                   <  2:587088/299‑1 (MQ=18)                                                                                                                                                                                                                                                    gtttGGCGTGATTTTTGCGGCAGGATAAGGCGTTTACGCCGCCGCCGACATCCACGCCTGCTGCGACCCTTAACGCGTCTTAGCAGGCCTACGCCAGAAAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttc                                                 <  2:350993/299‑1 (MQ=18)                                                                                                                                                                                                                                                        gtgttatttgtttcggTCTGATAAGTCGTGTATGCCGCATTCTAGATGCACGCGTGATGCGACGCTTACCGCGTTTAATCAGGCCTACGCCTGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                            <  1:39453/287‑1 (MQ=18)                                                                                                                                                                                                                                                        aGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTCTCACCAGTGCAGCGCCTTTTCCTTGTTTTtctc                                            >  1:264377/1‑301 (MQ=255)                                                                                                                                                                                                                                                            gaaTTTTGTGAGTCGGAGAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATa                                        <  2:118156/299‑1 (MQ=255)                                                                                                                                                                                                                                                             gtgtgtgtggcatgaaTAAGGCGTATATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCGTATCAGGCCTACGCCCGACAGTGCACTAGCCTGATTTAGCGTGATTTTGTCGGTCTGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaa                                       <  2:981014/286‑1 (MQ=21)                                                                                                                                                                                                                                                                 gtGTATGTCGGAGAAGGCGTTTTGTCCCCATCCGACGTCAACGCCGTATGCGCCGCTAACCGCGTCTTATCAGGCCTACGCCAGACAGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAc                                    <  2:293007/299‑1 (MQ=21)                                                                                                                                                                                                                                                                    cccctcccccacacgaccctcctccaaTCTGACATCACCGCCTGATGCGACGTGTAACCCGTCTTCTCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTTTAGGTCAGATAAGGCGTTTCCG‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                <  1:300649/275‑1 (MQ=21)                                                                                                                                                                                                                                                                     aGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATTCGACGCTGGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACATCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAATGCAGCGCCTTTTGCTTGTCTTTCCTGATAAAAActttt                               >  1:968678/1‑299 (MQ=18)                                                                                                                                                                                                                                                                      ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTCTCAAACCCTTCAATACAGCTCAGGTTCAACCCGTTAAGCGCCCTTTCCTAGTTCTTCTTGTTAAAACcttatt                              >  1:978536/1‑301 (MQ=18)                                                                                                                                                                                                                                                                         cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCCCAGATTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCCCGATAAACCcttcttacc                           >  1:1111984/1‑298 (MQ=255)                                                                                                                                                                                                                                                                              aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGttt                      <  2:621981/301‑1 (MQ=255)                                                                                                                                                                                                                                                                               agacGTTTATGACGCATCCGACATCAACGCCTGGTTCGACGCGTAACCCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACGCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTa                     <  2:157077/298‑1 (MQ=255)                                                                                                                                                                                                                                                                               aGGCGTTGATGCAGCATCCGACAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTa                     <  2:456496/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                  cGTTGATGCCGCCACCGACATTAACGCCTGTTGCTAAGCTTAACCCGTTTTCTCAGGCCTACGACAGACAGCGCAATGGCCTGATGTAGCGGGATTTTGGAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGa                  <  1:33515/301‑1 (MQ=25)                                                                                                                                                                                                                                                                                   gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCATTCCAGCGCCTTTTGCTTGTTCTTCTCGTTAAAACCTTCTTCCCCTTTTaaaa                  >  1:510468/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                    tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTACCACCCGGTATCACACCCTTCCATACACCTCACGTTCAACCAGTGCGGCGCCTTTTGCTTTTTTTtctctaatcggaagggggtttggtcgggac                  >  1:984244/1‑273 (MQ=17)                                                                                                                                                                                                                                                                                    tgtaagcccgaTCCGCCATCAACGCCGAAGGCGACGCTTAACGAGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                            <  2:984498/263‑1 (MQ=21)                                                                                                                                                                                                                                                                                        tGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTCACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACATCTCAGGTTCAACCACTGCAGCGCCTTTTTCTTTTTCTTCCAGATAACACCTTCTTCACCGATTAAGAAcatc             >  1:1010807/1‑300 (MQ=21)                                                                                                                                                                                                                                                                                                  gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCCATGACTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCCGTTTCAACCAGGGCAGCGCCTTTTGCTTGTTCTTCTCGATAATACcttcagatcggaagagcacacctttgacctccac  >  2:300737/1‑271 (MQ=18)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                         GACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCA  >  NC_000913/4295989‑4296576

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑