Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,268	0	T	C	20.0%	107.2 / 21.4	65	intergenic (+474/+168)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (34/18);  new base C (7/6);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAA  >  NC_000913/4295983‑4296561                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                             cctccAGACAGCGCCATGTCCTGTTTTAGCGTGATTTGGTAGGTAGGATAAGGCGTTTACGCCGCCTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                                                                                                                                                                                                                 <  2:787803/296‑2 (MQ=11)  cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACCCCGCATCCGActgc                                                                                                                                                                                                                                                                                               >  2:122550/1‑297 (MQ=21)   gcgagatggcgaaaatgCCTGATTTAGGTTGATTTTGTAGGTCGGAGAAGGCGTTTACGCCGCATCCGACATCAGCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaa                                                                                                                                                                                                                                                                                             <  2:366185/285‑1 (MQ=255)         cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCACCCCCTGATGCGCCGCTTACCGCGTCTACTCAGGCCTACCCCACACCGCACAAAATCCTGCTGTATCGTGACTTTCTAGGTCGGGTGCGGCGTTTACGCCTCATccataacaaatgcctt                                                                                                                                                                                                                                                                                       >  1:543900/1‑287 (MQ=11)                  gatatcatatggCGTGATTTTGTGGGTCGAATGAGTCGTGTACCCCGCATCCGACTACCAGTCCTAATGCTACGCTTGCCCCGTCTTATCAGGCCAACCCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTATGCCGCATCCTACATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGc                                                                                                                                                                                                                                                                              <  2:1005934/290‑1 (MQ=11)                  aGCCTGATTTGGCGTGAGTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCCACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                                                                                                                                                                                                                                               <  2:93831/300‑1 (MQ=255)                      tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGGTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGACTGATGGAACGCtttc                                                                                                                                                                                                                                                                          >  1:15766/1‑299 (MQ=255)                               gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCttt                                                                                                                                                                                                                                                                 >  2:307684/1‑300 (MQ=12)                                   ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTTCGCCGCATCCTACATCAATGCCTGAGCCGACGCTTGCCGAGTCTTATCat                                                                                                                                                                                                                                                             >  2:640917/1‑300 (MQ=255)                                                               cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGATTAACGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACCTCAATGCCTGATTGGACGCTTCCCGCCTCCTATTAGGcacaacttaacctattgttcatttttt                                                                                                                                                                                                                                 >  1:796615/1‑275 (MQ=11)                                                                      gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                                                                                                                                                                <  2:139575/247‑1 (MQ=17)                                                                      gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                                                                                                                                                                >  1:139547/1‑247 (MQ=17)                                                                                                                                                                    aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGTCCTATCTTAACCGtt                                                                                                                       >  2:587037/1‑301 (MQ=9)                                                                                                                                                                       ggTTTACGCCGAATCTGACAGCAGCGCGTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGGGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGt                                                                                                                    <  2:66890/300‑1 (MQ=11)                                                                                                                                                                                  cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTCACCGTTGGTTAATTTTTCa                                                                                                          >  2:264603/1‑300 (MQ=255)                                                                                                                                                                                    tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACCTCAATGCCTGATGCGACGCTTGCCGGCGC‑TTAACAGGCCTATCTTAACCGTTTGTTATTTTTTCaaa                                                                                                        >  1:852581/1‑300 (MQ=11)                                                                                                                                                                                         cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTCCG‑‑C‑‑CGCAGCCGACATCAATGCCTGATGCGCCGCTGGCCGCGTC‑TTACCAGGCCTATCTTAACCGTGCGTTCATTTTTGAAAcctttc                                                                                                  >  1:260104/1‑296 (MQ=11)                                                                                                                                                                                            gtaggCCTGATAAGACGCGGCAAGCGTCGCATCAGGCGTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTTGCGTGATTTTGTAGGTCGGATCAGGAGTTT‑‑ACGC‑‑CGCATCCGAAATCAATGCCTGATGCGCCGCTTGCCGCGTC‑TTATCAGGACTGTATTAACCGTTGGTTAATTTTGTAAAatctcaacc                                                                                               >  2:313214/5‑296 (MQ=11)                                                                                                                                                                                            cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                >  2:821941/1‑300 (MQ=255)                                                                                                                                                                                            cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGGGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACCTACATGCCTGATGCGACGTTTTCCCCGTC‑TTAACGGGCTTTATTTAACCGTTGGTTAATTTTTCAAACCCtttccc                                                                                               >  1:84723/1‑297 (MQ=11)                                                                                                                                                                                             aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCTCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTAacac                                                                                              >  2:630668/1‑301 (MQ=255)                                                                                                                                                                                                   tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGCCGCTGGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGtt                                                                                        >  2:651056/1‑300 (MQ=255)                                                                                                                                                                                                    gATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGCCACCAATGCCCGATGCGACGCTTGCCGCGTC‑TTATCAGGCCCATCTTAACCGTTGGTTAATTTTTCAACCTTTTCACACCCCGttt                                                                                       >  1:1160924/1‑299 (MQ=11)                                                                                                                                                                                                     aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                                                                      >  2:182394/1‑301 (MQ=255)                                                                                                                                                                                                           cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                                                                <  2:1123789/301‑1 (MQ=255)                                                                                                                                                                                                             cttgccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTT‑‑ACGC‑‑CGCCTCCGACATCAATGCCTGATGCTACGCTTGCCGCCTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACCCTTCACACCCGGTATCAAACCCtt                                                                              >  1:924945/6‑301 (MQ=255)                                                                                                                                                                                                                gccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CACATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCCGGCCTATCTTACCCGTTGGTTACTTTTTCAAACTCCTCAAACCCGGTATCCAACCCTTccc                                                                           >  2:1100026/3‑300 (MQ=14)                                                                                                                                                                                                                       ttttgacggCCGACGCCGGAAAGAGAATTAGCCAGATTTGGCGTGAGTTTGTAGGTCGTTTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                     <  1:307610/293‑1 (MQ=17)                                                                                                                                                                                                                        ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAGGCGTCTTATAGGGCTTACGCAACAAGGCGCAAAAGCCGGATTAAGGGTGTTTTGGTAGGGCGAAAAGGGGGTTA‑‑AGCC‑‑CGCAT‑CGAAAACAAAGGCTTAAGGGAAGGTTGTCCGGTCTTTATCAGG‑CTATCTTTAACGTTGGTTAATTTTTTAAAATATTCAAAACAGGTAACAAAACCTTCCAAACaagcac                                                                  >  2:791513/1‑296 (MQ=2)                                                                                                                                                                                                                        ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                                                                   >  2:802684/1‑301 (MQ=255)                                                                                                                                                                                                                         tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGCCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTTATACCCGCTATTAAACTCTTCCCTCCAGccc                                                                   >  2:127766/1‑298 (MQ=14)                                                                                                                                                                                                                          aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGACGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACTCTTAACGCGTCTTATCACGCCTACGCCATACACCGCAATAACATGCTTTCGCGTCATAGTGTAGACCGTAACATGTGTTT‑‑A‑‑CGTCCCATCGGCCATCAATGCCGGAAGCGACGCGTGCCGCGTC‑GTATCAGCCCCTTTCTAACCGACGCGTAATCTCTCAACCCCTTCACACCCgattgccagtccatccatcaaccccc                                                                  >  1:912238/1‑275 (MQ=2)                                                                                                                                                                                                                            cagtcatACGCCAGCAAGGGGTATTGCTTGTTTTTGCGTGATTGTGTAGGTCGGATAAGGCGGTCACGCAGCATCTGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                                                                <  1:747271/294‑1 (MQ=17)                                                                                                                                                                                                                            cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACCTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                                                               >  2:145900/1‑301 (MQ=255)                                                                                                                                                                                                                              ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTCACCGTTGGGTAAttt                                                                                                              >  2:566814/1‑252 (MQ=12)                                                                                                                                                                                                                                   ccacaagaaggcgaataaGCCTGATTTGGCTTGTTTTTTAAGGTCGGTTAAGTCGTTTATGCCCCAACCGCCACCAACGCCTGATGCGACGATTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAcc                                                         <  1:499008/283‑1 (MQ=17)                                                                                                                                                                                                                                    cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTCTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTACCAGGCCTATCTTAACCGTTCGTAAATTTTTCAAACTCTTCACTCGCGATATCAAACCCTTCCACACAGCTCAGtttacaccc                                                        >  1:1170816/1‑294 (MQ=14)                                                                                                                                                                                                                                    cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAACa                                                        >  2:499136/1‑300 (MQ=255)                                                                                                                                                                                                                                    cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGGTGGTTAATTTTTCAAACTCTTCACACCCGGTACCAAACTCTTCCATTCACCCCAGGTTCAACCac                                                       >  1:356671/1‑300 (MQ=255)                                                                                                                                                                                                                                     gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCCGt                                                      >  2:232925/1‑301 (MQ=255)                                                                                                                                                                                                                                             gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCCCAGGTTCAACCAATTCAgcg                                                >  2:1096935/1‑299 (MQ=255)                                                                                                                                                                                                                                             gcacagtAGCCGGATTTATCGTGTTTTTGTGGGTCGAGTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCGGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg                                                <  1:1096659/293‑1 (MQ=18)                                                                                                                                                                                                                                             cctcAATACGCTGTTGTAGCGTGATTTGTTAGGTCGGGTAAGGCGTTTATCCCGCTACATCAATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc                                              <  1:205074/298‑1 (MQ=18)                                                                                                                                                                                                                                                    aGCCTAATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGc                                        <  2:356762/300‑1 (MQ=255)                                                                                                                                                                                                                                                       cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACAGCTCAGGTTCAACCAGTGCCGCGCCTTTTGCTTGt                                    >  2:210544/1‑301 (MQ=255)                                                                                                                                                                                                                                                           ttggtttggattgtggagcgcGGATAAGGCTTTTATTCCGCCTCCCACATCAGCGCCTGATGCGACGCGTCACGCGTCTTATCGTGCCTACGCCAGACAGCGCAATAGGCTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttct                                 <  1:650882/280‑1 (MQ=18)                                                                                                                                                                                                                                                               gagtgtTTTTGTAGTGCGTATAAGGCTTTTATGCCGCATCCTCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCCACGCCAGACACCGCAATAGCCTGGTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCt                                       <  1:347504/288‑1 (MQ=18)                                                                                                                                                                                                                                                               gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCACACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                     >  1:459930/1‑260 (MQ=255)                                                                                                                                                                                                                                                               gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCGTTCCATACAGCTCAGTTTCCACCAGTACCGCGCCTTTTGCtagatacgccg                             >  2:347595/1‑290 (MQ=17)                                                                                                                                                                                                                                                               gCGTGATTTTGTAGGTCGAATAAGGCGTTTATGCCACATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                            <  2:1161214/301‑1 (MQ=255)                                                                                                                                                                                                                                                                cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCCGTGCAGCGCCTTTTGTCTTTTCTTTTCg                            >  1:1154049/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 gTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCCGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCCACCAGTGCCGCGCCTTTTGCTTGTTCTTCTCGAt                          >  1:88389/1‑301 (MQ=255)                                                                                                                                                                                                                                                                  tGATTTTGTAGGTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGAAATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACACTTCCATACAGCTCAGACTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTAGAt                          >  2:168850/1‑300 (MQ=18)                                                                                                                                                                                                                                                                  ggATTTTGTAGGTGGGATAAGGGGTTTATGCCGCATCCGCAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATa                         <  1:232870/300‑1 (MQ=255)                                                                                                                                                                                                                                                                   gATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCCCACCTGGTATCAAACCTTTCCATACAGCTCAGGTTCAACCACGGAAGCGCACTTTGCTTTttctttttgtta                         >  1:941053/1‑293 (MQ=17)                                                                                                                                                                                                                                                                       ttGTAGGTCGGATAAGTCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAc                     <  2:941290/300‑1 (MQ=255)                                                                                                                                                                                                                                                                       ttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAc                     <  2:88408/300‑1 (MQ=255)                                                                                                                                                                                                                                                                         gTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACct                   <  2:1154339/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          tAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                     <  2:152652/249‑1 (MQ=255)                                                                                                                                                                                                                                                                           aGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACACGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttct                >  2:561094/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                     ggacGTTTGTGCCGCATCCGACCTTAAAGCCTGATGCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAAGTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTa      <  2:260168/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                        catttATGCCGCATCCTGCATCAACGCCTGATGCGCCACTTACAGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGa   <  1:182357/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                        cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTATTCACCGTTTa      >  1:1162678/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                         gTTTATGCCGCCTCCGACATCACCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa  <  2:852805/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                         gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCCAACCCTTCCTTACAGCTAAGGCTCAACCCGTTCCGCGCCTTTTGCTTTGTTTTCTCGATAAAACCTCCTTCCCCGTTTAAGa   >  1:862484/1‑300 (MQ=18)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                             ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAA  >  NC_000913/4295983‑4296561

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