Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,268	0	T	C	23.5%	39.8 / 13.0	35	intergenic (+474/+168)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (11/15);  new base C (2/6);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

GCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAA  >  NC_000913/4295979‑4296561                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                            gCCTACGCCAGACAGCGCAATAGCCGGAGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                                                                                                                                                                                   <  1:291246/300‑1 (MQ=25)  ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                                                                                                                                                                                                 <  1:396170/301‑1 (MQ=25)        ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaat                                                                                                                                                                                                                                                                                           >  1:453379/1‑300 (MQ=21)                                 gctgtaTTGTGTAGTGCGGATAAGGCGCTCACGCCCNCTCTGACATTACCGCCTGCTGCTACGCTTGCCTGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCGGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                                                                                                                   <  2:73353/295‑1 (MQ=11)                                    tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                                                                                                                                                                <  1:140283/300‑1 (MQ=17)                                             gtggtgaaaaagacgtgacccccgtaacctaCATCAACGCCTGATGCGCAGCTGGCTGCGTATTATCGGGCCTACGCCTGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCTTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                                                                                                                                                       <  2:453369/270‑1 (MQ=11)                                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAACCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCGCCGCATCCGACATCAATGCCTGATGCGACGCTTTCCCCGTCTTATCAGGCCTatcttaatcgtttgt                                                                                                                                                                                                                                         >  2:267949/1‑287 (MQ=11)                                                                                                                                                                          gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCTCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                                                                      <  1:431763/300‑1 (MQ=9)                                                                                                                                                                           cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCAACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                                                                     <  1:103126/300‑1 (MQ=255)                                                                                                                                                                           cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTCGGt                                                                                                                    >  2:57447/1‑301 (MQ=9)                                                                                                                                                                             ttgtagccgcgacgaaCATAAACTCCTGATGCGACGGTTGCCGCGTCTTACAAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                                                                  <  2:99659/286‑1 (MQ=11)                                                                                                                                                                               aacgcagCATCACACATCAAGGCCCGAAGCGACGCTGGCGTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGTTTTGGCGTGATTTTGTAGGTCGAATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATTCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                <  1:347658/295‑1 (MQ=11)                                                                                                                                                                                ccgccgccACCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                                                               <  1:57448/299‑1 (MQ=255)                                                                                                                                                                                           gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCCTAACGCGTTTTATCAGGCCTACGCCAGACCTCGCAATAGCCCGATTTAGCGTGATTTTGTCGGTCGGATAAGCCGCTTACGC‑‑CGC‑ATCCGACATCAATTCCTGATGCGATGCTTG‑CCGCGTCTCTACAGGCCTTTTTTAACCGTGTGTTAATTTCTAAAAct                                                                                                      >  1:243504/1‑299 (MQ=11)                                                                                                                                                                                              tctacaccccgaagagttgctagtcgagtgttATCAGGCATAGACAAGAAGGAGCAATAGCCTGATTTGGCGTGATTTTGTATGTCGGATATGGAGTTTATGCCGCATACGACATCAACGCCTGATGCGGCGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                                                                 <  2:373351/271‑1 (MQ=14)                                                                                                                                                                                                     ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATACCCTGATTTCGCTTGATTTTGGAGGTCGGAGAAGGCGTTTA‑‑CG‑CCCAATCTGCCATCAATGCCTGCTGCGAC‑CTTGACCGCTTTCTATCCGGCCTATGTTACACGTTGGTTTATTTTTCatatctttcttcatctt                                                                                          >  2:267357/1‑285 (MQ=1)                                                                                                                                                                                                        gaggCGACGCTTGCCGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCATATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                                        <  1:191091/297‑1 (MQ=255)                                                                                                                                                                                                               cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATACGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACTATTACGGCGGCTTATCAAGACTATGCCCGACAGCCCAATAGCCTGATTCAGGGTGACTTTGTTGGACGGATCACACGTTTA‑‑CGACGC‑CTCCTACATGAATACCCGATGCCACCCTTT‑CCGCATCATAACATGCATATCCTCCACCTTGGATAATTGTTCaccctcatacccgcaggtttcacaca                                                                                 >  2:490657/1‑275 (MQ=2)                                                                                                                                                                                                               cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTC‑ACGCGTCTTCTCAGGCCTATATTAACCGTTGGTTAATTTTTAAAACTTTTAACCCCCAGTTTCAAAcc                                                                                 >  2:490661/1‑300 (MQ=11)                                                                                                                                                                                                                    tcagacgtgtgctcttccgatctcAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCACACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                                                           <  1:394301/278‑1 (MQ=17)                                                                                                                                                                                                                      cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGTCAGACAGCGCAATAGCCTGCTTTATCGTGATTTTGTAGGTCGGATAAGTCGCTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACCCTGT‑CCGCGTCTTATCAGGCCTATTTTAACCGTTGGTTAATTNTTCAAACCCTTCCCACCCGGtcacaacacccacacca                                                                         >  2:70553/1‑285 (MQ=14)                                                                                                                                                                                                                          tcttagacggCCTTCGCCACACGGCGCAATAGCCTGATTTGGCCGTATTTTGTTGGTCGGAATAGGCGTTTATGCCGCAGCCCAACTCAACGCCTTAAGCGACGCTTTACGCGTCTTATCAGGCCTAGGCCAGACAGAGGAATAGCTCGATTTAGCGTGAGTTTGTAGGTCGGTTAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                     <  1:488256/293‑1 (MQ=17)                                                                                                                                                                                                                                  ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACACTTCCATACAGCTCAGGtt                                                              >  1:489064/1‑300 (MQ=255)                                                                                                                                                                                                                                     cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaat                                                          <  1:98527/301‑2 (MQ=255)                                                                                                                                                                                                                                           cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCACACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTcaaagatcggaagaactgctgttagg                                                    >  2:394292/1‑276 (MQ=14)                                                                                                                                                                                                                                               aaggcgaaTAGCCTGATTAGGCGTGGATTTGAAGGGCGGAAAAGGGTAGTAGTCCGAATCAGACAGGAACGCCGGATGCGACGCTTAGCGCTTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg                                                <  2:21113/295‑1 (MQ=18)                                                                                                                                                                                                                                                   gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt                                            >  1:181074/1‑301 (MQ=255)                                                                                                                                                                                                                                                    ctttccctccacgccggtcttccgatctggTCGGAGAAGGCTTTCATGCCGCTTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCttt                                           <  2:145560/273‑1 (MQ=18)                                                                                                                                                                                                                                                          acTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTATGCCGCATCCGACATCAACGCCTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCt                                       <  1:267361/298‑1 (MQ=255)                                                                                                                                                                                                                                                            tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACATCTCCGGTTCAACCAGTGCAGCGCCTTTTGCTTGtt                                   >  2:24533/1‑301 (MQ=255)                                                                                                                                                                                                                                                                    cgggATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                            <  1:267953/297‑1 (MQ=255)                                                                                                                                                                                                                                                                                ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTAGCtcggaagagcacacgtctgaactc                >  1:145562/1‑277 (MQ=18)                                                                                                                                                                                                                                                                                           gCGTTTATGCAGCAGCCGACATAAACGCCGGATGCGACGCTTAACGCGTCTTATCAAGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa     <  2:181072/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                            cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGa   <  1:24533/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                             ttttATGCCGTATCGGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCATTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa  <  1:2813/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                            GCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAA  >  NC_000913/4295979‑4296561

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