Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,268	0	T	C	21.3%	88.5 / 12.6	49	intergenic (+474/+168)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (16/21);  new base C (7/3);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAA  >  NC_000913/4296004‑4296540                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                        tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                                                                                                                                                               >  1:835764/1‑275 (MQ=25)         cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGCTAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGAAAAAGCCTGTTTT‑AGCGTGATCTTGTAGGTCGGAAAAGGCGTGTACGTCGCAACCGACacta                                                                                                                                                                                                                                                                         >  2:657406/1‑270 (MQ=11)         cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGGGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATACGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATCAGGCGTTTACGCCGCATCCGACATCa                                                                                                                                                                                                                                                                         >  1:732187/1‑273 (MQ=255)             aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                                                                                    >  2:129053/1‑274 (MQ=21)                                     tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGA‑TGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                                                                          >  2:246676/1‑275 (MQ=17)                                          cgactcgttcATCAACGCCTGAGTGTACGCTTTCGTCTTCTTATCAGGCCTACGCCGGCATGCCCAATAGCCTGATTTAGCTGTATTTTGTGGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                     <  2:85373/266‑1 (MQ=11)                                            cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGA‑TGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                                                                                                                                   <  1:246676/275‑1 (MQ=14)                                             aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATACGACATTAACGCCTGATACGACGCTTAACGCGTCCTATCAGGCCTACGACAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGCAGGATAAGGCGGTTAAGCCGCATCCGACATAAATACCTGA‑TGAGACTCTTGACAGATATTATGAGGCCTa                                                                                                                                                                                                                                  >  2:773115/1‑275 (MQ=11)                                                   cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCg                                                                                                                                                                                                                                                                                     <  1:784795/219‑1 (MQ=14)                                                   cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCg                                                                                                                                                                                                                                                                                     >  2:784795/1‑219 (MQ=14)                                                                                                                                                                       cAACGCCTGAGGCGACGCTTACCACGTCTTATCGGCCCTACGCCAGACAGCGCAATGGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCACAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                                                                                    <  1:657406/275‑1 (MQ=11)                                                                                                                                                                                        cttgccgcgTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                   <  1:835861/270‑1 (MQ=255)                                                                                                                                                                                           gccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                                                 >  1:716941/3‑274 (MQ=255)                                                                                                                                                                                                gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGATCATTTTTCAAACTCTTca                                                                            >  1:349533/1‑274 (MQ=11)                                                                                                                                                                                                   ttATCAGGCCTACGCCCGACAGTGCAATAGCCTGAGTT‑AGCGTCATTTTGTAGGTCGAATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                        <  1:773115/275‑1 (MQ=14)                                                                                                                                                                                                    tATCAGGCCTACGCCAGTAAGCGCATGAGGCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TCCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                                       <  2:648074/275‑1 (MQ=14)                                                                                                                                                                                                    tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                        <  1:840681/274‑1 (MQ=255)                                                                                                                                                                                                    tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                        <  1:192399/274‑1 (MQ=255)                                                                                                                                                                                                         ggcctgcggcgtaCAGCCCAATGGTCTGCTTT‑GGCGAGATTTTCTAGGACGGACAAGGCCTTACGGGCGCATCCTAAATCAAACGCTGG‑TGCGACGCTTACCTCGTCTTATCAGGGCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCAGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                   <  2:425014/262‑1 (MQ=14)                                                                                                                                                                                                               cGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACCCCACCCAGCGCTATAGCCTGCTTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTT‑‑ACGC‑C‑GCGTCCAACATCACTGCCAGAGTCATGGTCTGACGCGTCTTAGCAGGCCTCTCTTAAACGTTGGTTAATATTTCCAGCTCTTCACACCCGGTacaccaca                                                            >  2:222411/1‑268 (MQ=11)                                                                                                                                                                                                                gCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                                           <  1:131812/275‑1 (MQ=255)                                                                                                                                                                                                                  gtgACAGCGCAATAGCCTGATTT‑GGCGTGATTTTGTAGATCGGATAAGGCGTTTATGCCGCATCCGACATCGACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGATCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                                         <  1:169818/273‑1 (MQ=255)                                                                                                                                                                                                                  cAGACAGCACAATAGCC‑GATTTGAGCTTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATACGACAGCAACGCCTGA‑TGCGACGCTTAACGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                                         <  1:642176/275‑1 (MQ=17)                                                                                                                                                                                                                        gcgcAATAGCCTGATTT‑AGCGTAATTTTGTATGTCGGATAAGGCGTTTATGTCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                                    <  1:124785/274‑1 (MQ=255)                                                                                                                                                                                                                         cgcAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                                                  >  1:678092/1‑275 (MQ=255)                                                                                                                                                                                                                         cgcAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGACGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTAACACCCGGTATCAAACCCTCCCAtcaa                                                  >  1:235206/1‑272 (MQ=255)                                                                                                                                                                                                                         cgcAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGAAATCAATGCCTGATGTGACGCTTTCCGCGTCTTATCAGGCCTATTTTAAGCGTTGGTTAATTTTTCAAACTCTTCACACCAGGTATCAAACCCTTCCAtca                                                   >  2:62036/1‑272 (MQ=14)                                                                                                                                                                                                                            gatcaccTGATTT‑AGCGTGATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATACGACATCAACCCCTGA‑TGCGACGCGTAACACGTCTTATCAGGCCTACGGCAGACAGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                               <  1:482222/270‑1 (MQ=17)                                                                                                                                                                                                                            aaatgCCTGATTT‑CGCGTGATTTTTGAGGTGGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCTACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                               <  2:296483/271‑1 (MQ=17)                                                                                                                                                                                                                            aaTATCCTGATTT‑ACCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                               <  1:504924/275‑1 (MQ=255)                                                                                                                                                                                                                            aaTAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                               <  1:129927/275‑1 (MQ=255)                                                                                                                                                                                                                            aaTAGCCTGATTT‑AGCGTGATTTTGTAAGTCGGATAAGGCGTTTATGCCGCATCCTACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                               <  1:275554/275‑1 (MQ=255)                                                                                                                                                                                                                                   tGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATTCAGTTCAGGtc                                         >  1:530396/1‑273 (MQ=255)                                                                                                                                                                                                                                   tGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTCTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAAAGCGTCTTATCAGGCCTATGCCAGACAGCGCAATAGCCTGATCCAGCGTGACTTTGTAGGTCGGATAAGGCGTTT‑‑ACTC‑C‑GAATCCGACATCAATGACTGATGCGACGCTTGCCGCGTCTTATCAGGTCTATCTTATCCGTTGGTTAATCTTTCAAAATCTTCACACCCGGTATCCAACCCTTACATACAGCTCAgttc                                         >  2:802513/1‑271 (MQ=17)                                                                                                                                                                                                                                    gATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGTCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTT‑‑A‑‑CGCTGCATTCGACCTCAATGACTGATGCGACGCTTGCCGTCTCTTATCAGGCCTATCTTAACCTTTTGTTCATTTTTCCAACTCTTTACACACTGTAACATATCCTTCCATAAAGcctacgtt                                         >  2:707176/1‑266 (MQ=14)                                                                                                                                                                                                                                    gATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTGGCCGCGTCGTATCAGCCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCGGGTATCAAACTCTTCCATACAGCTAAGGtt                                         >  2:223300/1‑273 (MQ=17)                                                                                                                                                                                                                                      ttt‑ACCGTGATTTTTAAGGTCGGATAAGGCGTTTATGCCGCATCCTACATCAACGCCTGT‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCG‑‑C‑C‑GCATCCAACATCAATGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCATTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                                     <  1:108020/275‑1 (MQ=17)                                                                                                                                                                                                                                       tt‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                                     <  1:129053/274‑1 (MQ=255)                                                                                                                                                                                                                                           gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                   >  1:162776/1‑257 (MQ=255)                                                                                                                                                                                                                                            cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GAATCCGAAATTAATGCCTGATGGGTAGCTTTACGCGTCTTATCAGGCCTAGCTTAACCGGTGGTTAATTTTTCAAACCCTTCACACCCGGTAAACAACCCTTCAATACAGCTCAGGGTAAACCAGtt                                >  2:311363/1‑274 (MQ=14)                                                                                                                                                                                                                                               gATTTTGTAGGTCGGTTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                               <  2:162776/257‑1 (MQ=255)                                                                                                                                                                                                                                                aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgc                            >  2:278194/1‑275 (MQ=255)                                                                                                                                                                                                                                                aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCATTCCATACAGCTCAGGTTCAACCAGTGCAg                             >  2:94195/1‑274 (MQ=255)                                                                                                                                                                                                                                                aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCGTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa                              >  1:779768/1‑273 (MQ=255)                                                                                                                                                                                                                                                      tAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCttt                      >  1:879790/1‑275 (MQ=255)                                                                                                                                                                                                                                                          gcGGTTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGAAGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGc                   <  1:417481/273‑1 (MQ=255)                                                                                                                                                                                                                                                                 atgCGTTTATGCCCACACCGCCACCAACGCCTGAGTGC‑TCGTGTACTGCGTCTTAGTAGCCTCAAGCCAGCCATAGCAATTGCCTGATTTAGCGTGATTTTTTAGGTCGAATATGTCGTTT‑‑ACGC‑C‑GCATCCTACATCAATGCCTGATTCGACGCTTGCCTCGTCTTATCAGGCCTATCTTAACCTTTTGTTAATTTTTCTAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttctt           <  2:757293/273‑1 (MQ=255)                                                                                                                                                                                                                                                                  ggCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttcttc          <  1:264921/275‑1 (MQ=255)                                                                                                                                                                                                                                                                        tATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATa    >  1:772130/1‑275 (MQ=255)                                                                                                                                                                                                                                                                          tGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaaa  <  1:94195/275‑1 (MQ=255)                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                        TGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGA‑TGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAA  >  NC_000913/4296004‑4296540

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