Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,268	0	T	C	26.0%	64.0 / 26.3	53	intergenic (+474/+168)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (21/16);  new base C (7/6);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATC  >  NC_000913/4295983‑4296565                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                              aCGCCAGACAGCGCAATAGCCGTGTTCAGAGTGATTTTTTAGGTCGGATAAGGCGTTTACG‑CCGCATCCTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAtn                                                                                                                                                                                                                                                                                                <  2:259262/301‑2 (MQ=255)   cccaacgcacgacatggCCTGATTTAGCGTGTTTTTGTAGGCCGGTTAAGGCGTTTACG‑CCTCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaa                                                                                                                                                                                                                                                                                              <  2:871100/285‑1 (MQ=255)         aaacgcAATATCCTGATTTAGCGTGATTTTGTAGTTCGGATAATGTATTTGCG‑CTGCATACAACACCAACGCCTGCTGCGACGCTTACCGTCTCTTATCAGGCCTACCGCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                                                                                                                                                                        <  1:509177/298‑1 (MQ=11)            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATCCctgctg                                                                                                                                                                                                                                                                                     >  2:643990/1‑301 (MQ=19)                            gtggTGATATTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                                                                                                                     <  2:888469/298‑1 (MQ=12)                              gatagactatatgcatcgctgagggcgttACGCCCGCATCCGCCAGCCACGCCATATGCGACGCTTTCCGGGGCTTGTCAGACATGTGCCAGACAACGCAATAGCCTAGTTTCGCGTGATGTTGGAGGTCGTATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                                                                                                                     <  1:350198/273‑1 (MQ=11)                              cGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACGCCAGACAGCGCCATAGCCCGATTTCTCGTGATTTTGTAGGTCGCATAAGGCGTTTACGCCGCATCCGGCATCAATGCCTGAGTCTACGCtcgccccgtctt                                                                                                                                                                                                                                                                   >  2:851525/1‑290 (MQ=11)                               gTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGACTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGACTGCTGCGACGATTGCCGCGTCttt                                                                                                                                                                                                                                                                  >  2:193501/1‑300 (MQ=255)                                    tttGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGGGTCTTATCAgg                                                                                                                                                                                                                                                             >  2:69283/1‑301 (MQ=17)                                          ggTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaat                                                                                                                                                                                                                                                                                             >  2:376922/1‑262 (MQ=25)                                           gccGGATAAGGCGTTTACG‑CCGAATCCGACATCAACGCCTGAGGGGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCACAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                                                                                                                                                                      <  2:227822/299‑4 (MQ=11)                                           gTCGGATGCGGCATAAACG‑CCTTATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGCGATGTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGACTGATGCGACGCTTACAGCGTCttcttcggca                                                                                                                                                                                                                                                           >  2:521823/1‑288 (MQ=11)                                            tCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGTGTTTACGCTCCCGTCGACAGCAAGGCCGGATGCGACCCTTCCCTCGTCTTTGCAGGCCtttctt                                                                                                                                                                                                                                                     >  1:229392/1‑296 (MQ=11)                                            tCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTCCTCCACCAATGCCTGATGCCACgtgtctcctgtcttatcgaaagagtccc                                                                                                                                                                                                                                                     >  2:71276/1‑274 (MQ=11)                                                                              aggcctgataagacgcggcaagcgTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATg                                                                                                                                                                                                                                                                                     <  2:214267/215‑1 (MQ=11)                                                                                                                                                                   taagtcgttctggcggcttctgatttCAACCCCTGTTGCGCCGCTTGCCGCGTCTTTTCGGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                              <  1:789924/275‑1 (MQ=11)                                                                                                                                                                              cgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATATGGCGTTT‑‑ACGCCGCATCCGACATCCATGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCGATCTTAACCGTTGGTTAAtgt                                                                                                                  >  1:446151/2‑299 (MQ=255)                                                                                                                                                                              cgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGAGGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCGTTGGTTAAttt                                                                                                                  >  2:242191/2‑301 (MQ=255)                                                                                                                                                                               gCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑AGGCCTCATCCGACCTCAATGCCTGATCCGACGCTTGCCACTTCGTATCAGGCCTATCTTAACCGGTGGTTGAtttt                                                                                                                 >  2:441330/1‑301 (MQ=11)                                                                                                                                                                               gCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCCCATCCGACCTCCCTGCCTGCTCCGACCTTTGCCGCCCCTTCTCAGGCCTATCGTAACCGTTGGTTAAtttt                                                                                                                 >  1:647238/1‑301 (MQ=11)                                                                                                                                                                                       cGACGTCAACGCCTGATGCGACGCTTGCCCGGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                                                                                                         <  2:446274/301‑1 (MQ=255)                                                                                                                                                                                                cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCTACGGTTGCCGCTGCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                                                                >  1:339902/1‑301 (MQ=255)                                                                                                                                                                                                         cGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCAGCCGACATCAATGCCTGATGCCACGCTTGCCGCGTCTTATCCGGCCTATCTTAACTGTTGGTTAATTTTTCAAACTCTTCACACCAGGTATCaa                                                                                        >  1:571203/1‑300 (MQ=14)                                                                                                                                                                                                         cGACGCTTACCGCTTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                       <  2:398219/301‑1 (MQ=255)                                                                                                                                                                                                              cttgccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGACGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCGTTGGTTAATTTTTCAATCTCTGCACACGCGTTATAAAACCCtt                                                                                  >  1:717189/6‑301 (MQ=14)                                                                                                                                                                                                                ggccgcgTCTTATCAGGCCTACGCCAGACAGCGTAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                                                                  <  2:374755/296‑1 (MQ=255)                                                                                                                                                                                                                  ccgcgTCTTATCAGGCATACGCCGAACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAt                                                                              <  2:821551/300‑1 (MQ=255)                                                                                                                                                                                                                     cggttgcgcggacctaagtcgaaccgcgcGATAGGCTGGTTAAACGTGGTTTTGTAGGCAGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGTTGCGACGCTTAACGCGTCTTATCAGGCCACCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTCCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                                                                           <  1:441209/276‑1 (MQ=17)                                                                                                                                                                                                                           aacAGGCCGACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                                                                     <  2:717388/299‑1 (MQ=255)                                                                                                                                                                                                                                  caacaccAGACAACGCAAAAGCCTGATTTAGCGTGAGATTGGAGTGCGGATAAGGCGTGTAAGCCACATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCAGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                                                              <  2:376054/296‑1 (MQ=17)                                                                                                                                                                                                                                    aCGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACCTCAATGCCTGCTGCGACGCTTGCCGCTTCCTATCAGGCCTATCTTAACCGTTGGTTCATTTTTCACACCCTACACCCCCGGTTACAAACCCTTCCATAAAGCTCAGGTCCAAcc                                                             >  1:767295/1‑300 (MQ=14)                                                                                                                                                                                                                                        gcagatggcctatgttgtGATTTAGGGGGGTTTGTTGGGTGGGCTATGGCGTTTCTGCCTCAGTCTACATCAACGCCTGATGCTACGCTTAACGCGTCTTATCAGGCATACGCTAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc                                                        <  1:604285/284‑1 (MQ=18)                                                                                                                                                                                                                                            cAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAGGCGACGCTTAACGCGACTTATCATTCCTACGCCAGACAGCGCAATAGCCTGATTCAGCGGGATATTGTAGGTCTGATAAGGCGGTT‑‑ACGCCGCATCAGAAATCAATGCCTGATGCGACGCGTGCCTCGTCTTATCAGGCCTATCTTAACCTATGTTTAATTTTGTAAACTCTTCACACACGATGTCAAACACTTCAATACCACTCAGGTTCacgcacttcac                                                      >  2:763701/1‑289 (MQ=14)                                                                                                                                                                                                                                            cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCGTCGGTTAATTTTTCAAACTCTTCACCACCGTGTTCAAACCCTTCCAAACAACTCCAGTTCAATCAGttaagcc                                                    >  1:662678/1‑295 (MQ=14)                                                                                                                                                                                                                                             aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATATAGCGTGACTTTATAGGTAGGATAAGGCGTTT‑‑ACGCCGCATCCTACATCAATGCCTGATGCTACTCTTGCCGCTTCTTATAAGGCCTAACTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAACTAAGGTTCCAcccttgctactc                                                   >  2:126048/1‑291 (MQ=17)                                                                                                                                                                                                                                                      gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCCTACAGCTCAGGTTAAACCAGAGCAGCGCCTCTTGCtt                                          >  1:135648/1‑301 (MQ=255)                                                                                                                                                                                                                                                       ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTTTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCGTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACGCGGTATCAAACCCTTCCATACAGCTCAGGTTCCACCAGTGCAGCGCCTTTTGCttt                                         >  2:685004/1‑300 (MQ=255)                                                                                                                                                                                                                                                       ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCttt                                                                                                      >  1:526560/1‑239 (MQ=255)                                                                                                                                                                                                                                                       ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                      <  2:526699/240‑1 (MQ=255)                                                                                                                                                                                                                                                        cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCACTGCACTGCCTTTTTCTTGt                                        >  1:463402/1‑301 (MQ=255)                                                                                                                                                                                                                                                           tttggccatGATTTTGTAGGTCGGATAAGGCGTTTAGGCGTAATCCGCCATCAACGCCTGATGCGACACTTAACGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttc                                      <  1:69274/292‑1 (MQ=18)                                                                                                                                                                                                                                                                gCGTGATTTTGGAGGTCGGATAAGGCGTTTAGTCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                <  2:662858/301‑1 (MQ=255)                                                                                                                                                                                                                                                                   tGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAACCTTCCATACCGCTCAGGTTCAACCAGTGAGCCGCCCTTTGCTTTTTCTTCTCGac                              >  1:396538/1‑299 (MQ=18)                                                                                                                                                                                                                                                                            gggtgtgATAAGGCGTTTATGCCGCATACGACATCAAGTCCTGATGCGACGCGTAACGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttct                    <  2:396641/295‑1 (MQ=255)                                                                                                                                                                                                                                                                             ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCTCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTACGGTTCAACCAGTGCAGCGCCTTTTGCTTTTTCTTCTCGATAAAACcttctt                   >  2:91402/1‑301 (MQ=255)                                                                                                                                                                                                                                                                              gTCGGATAATGGGTTTATGCCGCATCCGACACCAACGCCTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCttcttc                  <  2:695214/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                 ggttAAGGCGTTTATGCCGCATCTCACACCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAcc               <  2:229443/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTAGCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCAAACAGCTCATGTTCAACCAGTTCATCTCCTTTTGCTTGTTTTCCTCGATAAAAACTTCTTCagcc              >  1:656848/1‑298 (MQ=18)                                                                                                                                                                                                                                                                                        gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACGCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGAAGCGCCTTTTGCTTGTTCTTCTCTATAAACCCTTCTTCACCGTTTaat        >  1:888041/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                        gCGTTGATGTCGGATGCGGCATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACACTTCCCTACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCACGAGAAAACCTTCTTCACCGTTTac         >  1:287762/1‑299 (MQ=21)                                                                                                                                                                                                                                                                                         cGTTTATTCCGCACCCGACATCAACGCCTGATGCGACGCTTAACTCGTCTTATCAGGCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGa       <  2:888281/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                         cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGa       <  2:339986/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                          gtgtATGCCGCACACGACATCAACGCCTGAGGCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa      <  2:571357/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                               tGCGGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCCCCGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATACAACCTTCTTCACCGTTCAAGAAcatc  >  2:677602/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                              ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATC  >  NC_000913/4295983‑4296565

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