Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,059	C→T	13.1%	intergenic (+266/+378)	gltP → / ← yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,059	0	C	T	13.1%	79.2 / 4.6	43	intergenic (+266/+378)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (16/21);  new base T (4/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACG‑CCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT  >  NC_000913/4295795‑4296267                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                  tGAGGATGGCGGGGGATATCACCCCGGCAATGGGTGAATGCTTGGTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCATGACAGCGAAATAGCCTGAGTTAGCGTGAGTTTGTAGGTCGGATAATGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                                                                                                         <  2:306737/275‑1 (MQ=11) ggaggattgcccgtgtgatcgccccGGAAATGTGTGAATGCCTGAGGCGACGCTTGCCGCGTCTGTTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAGCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                                                                                                         <  1:288130/253‑1 (MQ=11)        tGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGAGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATACACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGACTGATGCgcagatgg                                                                                                                                                                                                  >  1:115560/1‑268 (MQ=255)        tGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATAAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCACGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCTCATCCGACATCAATGCCTGATCCGACGCTTg                                                                                                                                                                                                  >  2:357154/1‑275 (MQ=255)              gtgatgccccccGGGCATGTGGGAAAGCCGTGTGGGACGCTTGCCCCGTCTTATCAGGCCTACGCCAGACAACAAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAGGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt                                                                                                                                                                                            <  2:30052/266‑1 (MQ=11)                   tccacccaggaattgtggctggcctcgtGCGACGCTTGCCGCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGGGTGATTTTGTGGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                       <  2:25096/248‑1 (MQ=11)                   tacgccgcatccgacatcAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                                                                                        <  1:96047‑M1/256‑1 (MQ=255)                    ccgcccggaaacggtggAATGCTGGATGCCACTCTTTCCGCGTCTTCTCAGGCCTACGCCAGAGAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                       <  2:275092/258‑1 (MQ=11)                     cACCCGGCAATGTGTGAATGCCTGATGCAACGCGTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                                      <  1:307992/274‑1 (MQ=255)                      ccccGGCAATGTGTGAATGCCTGATGCTACTCTTGCCCGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                                      <  2:356333/272‑1 (MQ=255)                         cgggaattgtgtgagGCCTGGTGCGAGGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATGGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                  >  2:355686‑M1/16‑274 (MQ=255)                         cGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                 >  1:356334/1‑275 (MQ=39)                                             aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTATGCCGCACACGACATCAACGCCGGATGCGACGCTTATCGC‑GTCTTATCAGTCCTACGCCAAACAGTGCAATCGCCGGCCTTAGCGTGATTTTGTAGTACGTATACGTCGTTTAGGCCGCACACGCAATCAATGCCTGACGCGATGCTTGCCGCGTCGTAACAAGCCTatcacttcctcggggtta                                                                                                                                                              >  2:249668/1‑257 (MQ=11)                                             aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGC‑GTCTTATCAGGCCTACGCAAAACGCCGCAATAGCTTGCTTAAGCGTGATTTTTTGGTCCGGATAAGCCGTTTACGCCCCATCCCCCACAATTGCCTGCTGCGACGCTTCCCGCGTCTTTCCCGGCCtcttttaaccttgggttac                                                                                                                                                              >  2:130944/1‑256 (MQ=11)                                               cggagcctggcggggttttATGAGGCATACTCGAGACAGCGAATTGGCCTTATTTCGGGGGAATTTGTGGGTCGGATAAGGCGTTTATGCCGAATCCGACATCAACGCCTGATGCGACGCTTGACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACAGCGCaa                                                                                                                                                                <  2:369055/253‑1 (MQ=11)                                                  ACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGCCATCAATGCCTGATGGGACGcttgccgcgtcttgtcaggcctatctttacattttgttaatttgt                                                                                                                                                         <  1:321966‑M1/274‑23 (MQ=255)                                                                      tGCCTAAGCCAGACGGCCGAATAGCCTAATTTGGCGGTATTTTGTAGGTCGGGTGAGGCGTTTATGCCGCATCCGCAATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt                                                                                                                                    >  2:274353‑M1/2‑275 (MQ=255)                                                                         CTACGCCAGATAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATACGGCCTACGCCAGACAGCGAACTAGCTTGCTTTAACGTGATTTTGTATGTCg                                                                                                                                  <  1:492532‑M1/274‑1 (MQ=255)                                                                             gacaggaagtggaaTAGCTTGAGTTGGCATGATTGTGGGGGTCGGGGAAGGCGTTTGTGCCGGAGCCGGCAGCGAGGCCGGATGGGACGCTTACCGC‑GTCTTATGAGGCCTAAGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATa                                                                                                                              <  2:30082/262‑1 (MQ=11)                                                                             gCCAGACGGCGCATAAGCCTGATTTAGCGTGATTTTATGGGTCGGATAAGGCGTTTATGCCGCATCCGGCATCAACGCCTGATGCGACGCTTAACGC‑GTCTTAGCAGGCCGAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa                                                                                                                             <  2:162137/275‑1 (MQ=255)                                                                                 gACAGCGCAATAGCCTGTTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACCTCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACCCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTGTGTAGGTCGGATAAGGc                                                                                                                          >  1:317955/1‑274 (MQ=16)                                                                                                ggttttAGCGGGGTTTTGTAGGGCGAGTAGGGGGTTGATGCCGCATCCGACATCAAGCCCTGATGGGAAGGTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATcc                                                                                                           <  2:41693/271‑1 (MQ=11)                                                                                                 ggtttggcggtactttGTAGGTCGGATAGGGCGTTTACGCCGCATCCGACATCACCGCCGGATGCGACGCTTACCGC‑GGCTTATCAGGCCTACGCCATACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCg                                                                                                          <  1:95540/261‑1 (MQ=11)                                                                                                        cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCACACATCGCACTAGCCTCATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCGCCGCATCCGCCCTCa                                                                                                    >  2:260521/1‑273 (MQ=11)                                                                                                            aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCt                                                                                              >  2:221526/1‑275 (MQ=35)                                                                                                                  tAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGACTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTAAGCCAGCCAGCGCAACAGGCTTCATTCGCTTCACTTTGTAGGTCCGAACAGTACTTTCTGCCGCATCCGaactgcactcctcactt                                                                                         >  2:315588/1‑258 (MQ=11)                                                                                                                      tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACg                                                                                     <  1:221526/274‑1 (MQ=35)                                                                                                                         gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTa                                                                                 >  1:317585/1‑275 (MQ=35)                                                                                                                                cttcatcgCCCCAACCGGCATCAACGCTGGGTGCGCCAC‑TAACGCGGTCTGTTCTGGCGTACCGCACAGAGCGAGATAGCCTGATTTTGCGTTATTTTGTAGGGCGGAGAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt                                                                           <  2:70606/267‑1 (MQ=11)                                                                                                                                cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGCCATCAAAGCCTGATGCGACGCTTACCGCGt                                                                           >  1:445305/1‑274 (MQ=16)                                                                                                                                                            ggacgcaacgcgtgaggc‑gTTTTATCAGGCCGAGCGCAGGCGGCGAAATAGCCTGATTTGGTGTGATTTTGTAGGGCGGATAGGGCGTTTCCGCCGCAACCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACG‑CCAGACAGCGCaa                                              <  2:104723/258‑1 (MQ=11)                                                                                                                                                                         aaCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCAATCAGACATCAACTCCTGCTGCGACCTTTTCCTCGTCTTATCAGGCCTACC‑CCCGACAGCGCCAAAGctttctttat                                 >  2:184504/1‑265 (MQ=11)                                                                                                                                                                                                   gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACG‑CCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg       >  1:255848/1‑275 (MQ=33)                                                                                                                                                                                                        cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                                      >  2:92521/1‑96 (MQ=33)                                                                                                                                                                                                        cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt                                                                                                                                                     >  2:81069/1‑129 (MQ=38)                                                                                                                                                                                                        cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGTCCGCACA‑CGCTATAGCCTGCTTTAGCTTGCTTCTGGAGGTCGGATAAGGCGtttcc  >  2:151815/1‑273 (MQ=11)                                                                                                                                                                                                            aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                  <  1:92521/96‑1 (MQ=33)                                                                                                                                                                                                            aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                 <  1:81069/129‑1 (MQ=38)                                                                                                                                                                                                                                                      tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACG‑CCAGACAgc                                                  >  2:93128/1‑181 (MQ=35)                                                                                                                                                                                                                                                          ccgcATCCTACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACG‑CCAGACAGCGCaa                                              <  1:93128/181‑1 (MQ=18)                                                                                                                                                                                                                                                              aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                                                                                                   >  2:69725/1‑77 (MQ=39)                                                                                                                                                                                                                                                                  gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg                                                                                                                                               <  1:69725/77‑1 (MQ=39)                                                                                                                                                                                                                                                                         aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                 <  2:322423/68‑1 (MQ=33)                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                  TGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACG‑CCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAT  >  NC_000913/4295795‑4296267

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