Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,191	0	A	C	28.6%	35.3 / 18.1	42	intergenic (+397/+245)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base C (9/3);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

CGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA  >  NC_000913/4295899‑4296480                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    ggagattgcttAGGTGGGATAAGGCGTTAAGGCCGCCTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCtt                                                                                                                                                                                                                                                                                             <  2:38278/291‑1 (MQ=255)  gggagtttGTGGGTTGGATAAGGCGTTAAGTGCGCATCCTACATCAACGCCTGATGCAACGCTTAACGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAAGCCGTTTACG‑CCGCATCCGACATCATCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTa                                                                                                                                                                                                                                                                                            <  2:858256/296‑1 (MQ=11)          gTAGGTCGGATCAGGCGTTCATTCCTGAGCGGTCAGCAACGCCTGATGCGACGATTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGACTGATTTAGTGTGATTTTTTAGGTCGGATAAGGCTTTTACG‑CCGCATCCGACATCAACGCCTTATGCTACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGCTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                                     <  2:42528/300‑1 (MQ=11)          gTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGTCAGGATAAGGCGTTTATGCCGCATCCGACATGAACGCCTGATGCCACGCTTAACGCGCCTTATCAGGCCt                                                                                                                                                                                                                                                                                    >  1:31912/1‑301 (MQ=255)             ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCGGGCCTa                                                                                                                                                                                                                                                                                   >  2:884883/1‑299 (MQ=25)             ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACCCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGGCTGAATTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGGCATCAACGCCGGATGCACAGCTTACCGCCTCTTATACGGCTTAcc                                                                                                                                                                                                                                                                                 >  2:92840/1‑300 (MQ=11)                       ggCGTTGATGTCGGATGCGGCATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCATGAAGCGACGCTTAACGAGTCGTATAAGGCCAACGCCAGAAAgcg                                                                                                                                                                                                                                                                       >  2:350282/1‑301 (MQ=11)                              aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACCTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCTAGCGAAATAGACTGAATTAGCGTGTTTTTGTCGGTAGGATAAGGCGGTTTTGCCGCATCGGAAATCAACGCCGCAGTAGCCGCTTACCGCGTCTTTTTAGCCTTGCGCCAGACCgcgcctacgc                                                                                                                                                                                                                                                                >  1:888229/1‑295 (MQ=11)                                        cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTGAGc                                                                                                                                                                                                                                                      >  2:405622/1‑301 (MQ=21)                                         atcATCAACTCCTGATGCTACGCTTAACGCTGCTTATCAGGCCTAAGCCATACAGCGAAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                     <  2:749873/299‑1 (MQ=255)                                         aacacCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                     <  2:31916/300‑1 (MQ=21)                                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACACCGCAATAGCCTGATTTAGCGTCACTTTGTCGGTCGGATAAGGCGTTTATGCCCCATCCGACATCAACGCCGGCTGCGTCTCTTAACGCCTCTTATCAGGCCTACCCACAACCGCCCAATCGCCCGCTTTTGCGTGAttt                                                                                                                                                                                                                                               >  1:830103/1‑301 (MQ=11)                                                   ccTTATGCGACGCTTAACGCGTCTTATCAGGCTTACGCCAGACAGCGCAATGGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                                                                                                                                                                                           <  2:781955/301‑1 (MQ=21)                                                    cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CAGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                                                                                            >  2:299149/1‑299 (MQ=25)                                                      gataggACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                     <  2:845316/283‑1 (MQ=255)                                                      gataagACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAACAGGCCTACCCAAGACAACCCAACAGCCTGATTTAGCgaacactgtagtgc                                                                                                                                                                                                                                        >  1:845088/6‑288 (MQ=11)                                                          cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCTGATTCGCCGCTTAACGCGTCTTTTCAGGCCTATCCACAACAGCCCAATCTCCTGATTTAGCGTGCTTTTGTgggcgggg                                                                                                                                                                                                                                    >  1:695020/1‑296 (MQ=11)                                                                                tgCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGTCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGa                                                                                                                                                                                                                                                            <  2:444077/254‑1 (MQ=18)                                                                                     aCGCCAGACAGCGCAATAGCCGTGTTCAGAGTGATTTTTTAGGTCGGATAAGGCGTTTACG‑CCGCATCCTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAtn                                                                                                                                                                                                         <  2:259262/301‑2 (MQ=255)                                                                                       cccaacgcacgacatggCCTGATTTAGCGTGTTTTTGTAGGCCGGTTAAGGCGTTTACG‑CCTCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaa                                                                                                                                                                                                       <  2:871100/285‑1 (MQ=255)                                                                                             aaacgcAATATCCTGATTTAGCGTGATTTTGTAGTTCGGATAATGTATTTGCG‑CTGCATACAACACCAACGCCTGCTGCGACGCTTACCGTCTCTTATCAGGCCTACCGCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                                                                                 <  1:509177/298‑1 (MQ=11)                                                                                                cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATCCctgctg                                                                                                                                                                                              >  2:643990/1‑301 (MQ=19)                                                                                                                gtggTGATATTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                              <  2:888469/298‑1 (MQ=12)                                                                                                                  gatagactatatgcatcgctgagggcgttACGCCCGCATCCGCCAGCCACGCCATATGCGACGCTTTCCGGGGCTTGTCAGACATGTGCCAGACAACGCAATAGCCTAGTTTCGCGTGATGTTGGAGGTCGTATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                              <  1:350198/273‑1 (MQ=11)                                                                                                                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACGCCAGACAGCGCCATAGCCCGATTTCTCGTGATTTTGTAGGTCGCATAAGGCGTTTACGCCGCATCCGGCATCAATGCCTGAGTCTACGCtcgccccgtctt                                                                                                                                                                            >  2:851525/1‑290 (MQ=11)                                                                                                                   gTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGACTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGACTGCTGCGACGATTGCCGCGTCttt                                                                                                                                                                           >  2:193501/1‑300 (MQ=255)                                                                                                                        tttGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGGGTCTTATCAgg                                                                                                                                                                      >  2:69283/1‑301 (MQ=17)                                                                                                                              ggTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaat                                                                                                                                                                                                      >  2:376922/1‑262 (MQ=25)                                                                                                                               gccGGATAAGGCGTTTACG‑CCGAATCCGACATCAACGCCTGAGGGGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCACAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                                                                               <  2:227822/299‑4 (MQ=11)                                                                                                                               gTCGGATGCGGCATAAACG‑CCTTATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGCGATGTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGACTGATGCGACGCTTACAGCGTCttcttcggca                                                                                                                                                                    >  2:521823/1‑288 (MQ=11)                                                                                                                                tCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGTGTTTACGCTCCCGTCGACAGCAAGGCCGGATGCGACCCTTCCCTCGTCTTTGCAGGCCtttctt                                                                                                                                                              >  1:229392/1‑296 (MQ=11)                                                                                                                                tCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTCCTCCACCAATGCCTGATGCCACgtgtctcctgtcttatcgaaagagtccc                                                                                                                                                              >  2:71276/1‑274 (MQ=11)                                                                                                                                                                  aggcctgataagacgcggcaagcgTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATg                                                                                                                                                                                              <  2:214267/215‑1 (MQ=11)                                                                                                                                                                                                                                                       taagtcgttctggcggcttctgatttCAACCCCTGTTGCGCCGCTTGCCGCGTCTTTTCGGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                         <  1:789924/275‑1 (MQ=11)                                                                                                                                                                                                                                                                  cgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATATGGCGTTTACGCCGCATCCGACATCCATGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCGATCTTAACCGTTGGTTAAtgt                             >  1:446151/2‑299 (MQ=255)                                                                                                                                                                                                                                                                  cgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAGGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCGTTGGTTAAttt                             >  2:242191/2‑301 (MQ=255)                                                                                                                                                                                                                                                                   gCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCTCATCCGACCTCAATGCCTGATCCGACGCTTGCCACTTCGTATCAGGCCTATCTTAACCGGTGGTTGAtttt                            >  2:441330/1‑301 (MQ=11)                                                                                                                                                                                                                                                                   gCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCCATCCGACCTCCCTGCCTGCTCCGACCTTTGCCGCCCCTTCTCAGGCCTATCGTAACCGTTGGTTAAtttt                            >  1:647238/1‑301 (MQ=11)                                                                                                                                                                                                                                                                           cGACGTCAACGCCTGATGCGACGCTTGCCCGGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                    <  2:446274/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                    cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCTACGGTTGCCGCTGCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc           >  1:339902/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                             cGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAGCCGACATCAATGCCTGATGCCACGCTTGCCGCGTCTTATCCGGCCTATCTTAACTGTTGGTTAATTTTTCAAACTCTTCACACCAGGTATCaa   >  1:571203/1‑300 (MQ=14)                                                                                                                                                                                                                                                                                             cGACGCTTACCGCTTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa  <  2:398219/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                    CGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACG‑CCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAA  >  NC_000913/4295899‑4296480

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