Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,190	0	A	G	32.4%	27.0 / 20.1	36	intergenic (+396/+246)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (5/6);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4295899‑4296482                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                           cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTc                                                                                                                                                                                                                                                                                                     <  1:346377/299‑1 (MQ=35)           gatgtcggatgcggcataaacgccttaTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTTTACGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATCAAGCGTGATTTTGTAGGTCGGATAGTGCGTTTGTGCCGCATCCGACATCACCGCCTGATGTGACGCTTAACGCTTCTTATCAGGCCtt                                                                                                                                                                                                                                                                                         >  1:133823/27‑300 (MQ=11)                   aTAAGGGGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATAAGACGCGTTAAGCGTCGCATCAGGcgttgatgtcg                                                                                                                                                                                                                                                                                  <  1:468295/300‑11 (MQ=11)                   aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCCACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACCCCCGa                                                                                                                                                                                                                                                                                  >  1:94137/1‑300 (MQ=25)                        tttccgatttcGCATCCGACATCAACGCCTGATACGACGCTTAACGCGTCTTATCGTGCCTACGCCAGACAGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAgcac                                                                                                                                                                                                                                                                            <  1:359999/291‑3 (MQ=255)                         aatttCTGCCGCATCCACAATCAACGCCTAATGCGGCGCTTAACGCGGCTGATCAGGCCTACGCGAAACAGCGTAATAGCCTAATTTAGCGTGATTTTGTAGGTCGGATAGGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAgcgc                                                                                                                                                                                                                                                                            <  1:50234/298‑1 (MQ=11)                                  cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATCAACCCCTGATGCGACGCCTGCCGCGTCTTATCAGGCCTCCGCCAGACAGCGATACAGCATGATTTAGCTTGACTGTGTTGGTCGGATAAGGAGTTCATGTCGCCTCTGACATTAACGCCTGATGCTACGCCTTTCGCGCCTTCACAAGCTTCCGCCATACacagcactatgcata                                                                                                                                                                                                                                                                  >  2:359992/1‑287 (MQ=11)                                       ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGACGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGTCATATTTa                                                                                                                                                                                                                                                              >  1:415729/1‑300 (MQ=25)                                              caagtgcctggtgtaCGCGTAAGGCATCTTATGGGACCACCGCCAGACAGCACAATAGCGAGAGTTAGCGGGATTGTGTAGGTCGGATAGGGCGTTTACGCCGCATCCAACAGCAACGCCTAATGCGACCCTTGCCGCGTCTGATCAGGCCGACGCAAGACAGCGCAATAGCCTGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                           <  1:59131/282‑1 (MQ=11)                                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGATTCGACGATTAAGGCCTCGTATCAGGCCTACCCCACACACCGCATAACCCTCATTTCGCGTGttt                                                                                                                                                                                                                                                      >  2:191089/1‑301 (MQ=11)                                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATGCCAGACAGCGCAATAGCCTGATTTACCGTGAttt                                                                                                                                                                                                                                                     >  1:373359/1‑301 (MQ=21)                                                    cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCCACTCCTGATGCGACAATTAACGCGTCTTATCAGGCCTTCGTCAGACAACGCACTAGCCTTAttgtacgaattttactag                                                                                                                                                                                                                                                >  2:431753/1‑285 (MQ=11)                                                      tgtgatagctcaacagagTCTCATAAGGCCAACGCCAGAAATCGTAATAGGCTGGTTTAGAGTACATTTGTAGGTCGGATAAGGCGTTTAAGCGCCATCCGACATCAACGCCTAATGCGACGCTTGCGGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAgg                                                                                                                                                                                                                                               <  2:227140/283‑1 (MQ=11)                                                        tGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCg                                                                                                                                                                                                                                            >  1:73354/1‑301 (MQ=34)                                                          cGACGCTTAAGTGGTCTTATCAGTCCTGCGCGAAAAAGGTCATTAGCTTGATTTAGCGTGATTTTGTGGGTCGGATGGGGCGTTTACGCCGTATCAGACATAAACGCCTGATGCTACGCTTGCCTCGTCTTATCAAGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGATCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACACTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCgg                                                                                                                                                                                                                                           <  1:240062/300‑1 (MQ=11)                                                          cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCATACGACAGACAGCGCAATAGCCTGGTTTAGCGTGATTTTGTAGGGCggg                                                                                                                                                                                                                                          >  2:140281/1‑300 (MQ=19)                                                                 tAACGCGCCGTATCAAGACTACGACAGACAGCGAAATAGCCTGAATGAGCGTGATTTTGTAGGTCGGATACGGCATTTACTCCGCATTCGACATCAACGCCTGATGCGACGGTTGCCGCGTCTGATCAGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCATGATTTTGTAGGTCGGATAAGGc                                                                                                                                                                                                                                    <  1:337402/300‑1 (MQ=11)                                                                          ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCAGCCGACATCAACGCCTGCTGCTTCGATTACCGCGTCTTATCAGGCCTCCCCCAGACACCCCAATAGCCTGCACTCGCGTGATATTCTAGCTCggtatcgggttttaagccg                                                                                                                                                                                                                          >  2:346446/1‑284 (MQ=11)                                                                                 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)                                                                                                                                                                                                                                                                                     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)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                           CGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4295899‑4296482

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