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

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	4,296,268	0	T→C	29.3%	14.9 ‑1.4 ‑0.5	35	intergenic (+474/+168)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (8/2):  ref base (11/10):  total (19/16)
Fisher's exact test strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.43e-01

GGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACAT  >  NC_000913/4295978‑4296564                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                             ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGCACACCAGCGCCATAGCCTGATTTAGCTTGATGTTGTAGGTCGGATCAGACGTTTACGCCGCATcc                                                                                                                                                                                                                                                                                                     >  2:217666/1‑300 (MQ=11) ggCCTACGCCAGACAGCGCAATAAACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCTACATCAACCACTGGTTCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACTGCCGCATCACCTACTTTAGCGTGCATTTCTAAGTCGTATCACGCGGTTACGACGccccc                                                                                                                                                                                                                                                                                                     >  2:193924/1‑296 (MQ=11)      aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCCAATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTCGCGTGATTTTTTAGGTCGGATACAGCTTTTACGCCGCATCCGACAt                                                                                                                                                                                                                                                                                                >  2:109495/1‑300 (MQ=255)                gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTTATGCCGCATCCGACATCAACGCCTGATGCGGAGCTTAACGCGACTTATCAGGCAGACGCAAGACAGAGCAATAGCCGCATCTATAGTGATTTGGTATGTCGGCGAAGACGTTGACGCCGCATCCGaactcactgaccacg                                                                                                                                                                                                                                                                                     >  2:186216/1‑287 (MQ=11)                   cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCTGATGCGACGCTTCACGCGTCTTATTAGGCCTACGCCACACACCGCAACAACCTGATCTCACGTTATTTTTTAGGGCCGCTACGGCCTTTACGCCGCATCCGACCTCACAGCCttcttc                                                                                                                                                                                                                                                                                   >  1:409053/1‑295 (MQ=11)                                       aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCCGACAGCGCCATCGCCTGCTTTCGCGTGCTTTTGTCGGTCGACTACGGCGTTTACGCCGCCTCCGACATCAATGCCTCATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                                                                                                               >  2:233704/1‑300 (MQ=11)                                       aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAGGCTTAACGCATCTTATCAGGCCTACGCCCGACCGCACAATAATCTGATTTAGCGTGCTCTTGTCGGTGGGTTCATGCGTTACCGCCGCGTCCGACTTCAACGCCTGAGTCGGCGCTTGCCGCGTTTTacat                                                                                                                                                                                                                                                              >  2:368268/1‑298 (MQ=11)                                                      tAAGGCGTTTACGCCGGATTAGAGACCAAGGGCGGAGGGGAGGCTGGCCGCGTGGTATCAGGCCGACGACAGAAAGCGCAATAGACGGCTTGAGCGTGATTGTGTAGGTCGGAGAAGGCGTTTATGACGCATCAGAAAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                <  1:312151/300‑1 (MQ=1)                                                                                                                                                                     gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTACCC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaca                                                                                                                              >  1:137159/1‑299 (MQ=2)                                                                                                                                                                       tagggCGTTTATGCCGGAGAGGACATAAACGGCTGGGGCGAGGGTGTGGGCGTCTGATCAGGACGACGCCAGACGGCGCAATAGCCGGATTTGGCGTGATTTTGGAGGTCGGATAAGGCGTTTAGTCCGCAGCCGACATAAACGCCTGATGCGACGCTTAATGCGTCTTATCAGGCCTACGCCAGACAGGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                             <  1:416331/297‑1 (MQ=11)                                                                                                                                                                          ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGACAAGGCGTTTATGCCGCGTCCGACATAAACGCCTGATGAGACGCTTAACGCGTCTTATCAGGACTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAGTTTGTAGGTCGGGTGAGGCGTTTAAGC‑C‑GCATCCGACATTAATGCGTGATGCGACGCTTGCAGCGTCTTGTCAGGCCTATCTTAACCGTTg                                                                                                                         >  2:161542/1‑301 (MQ=11)                                                                                                                                                                                      gCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTC‑‑CGCGGCATCCGACATCAATGCCTGATGAGACGCTTGCCGGGTCTTATAAGGCCTAGCTTAAACGTTGGTTGAtttttaa                                                                                                             >  2:93286/1‑299 (MQ=11)                                                                                                                                                                                               aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGCTCAGGCGTTTACGC‑C‑GCATCCGACATCCATGCCTGATGCGCCCCTTTCCGCGTCTTTTCATCCCTATTTACACCGTGCGTTCATTTTTTAACctttct                                                                                                     >  1:400475/1‑296 (MQ=1)                                                                                                                                                                                                   cggtctgtggggactgttgacggtacTTATCAGGCCGACGCCAGAAAGCGCAAGAGCCGGTGTTAGCGGGATTTTGTAGGTCGGATAAGGCGTTTATGACGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTATTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                <  1:25112/276‑1 (MQ=14)                                                                                                                                                                                                   aCGCCTGATGTCAAGCTGTCCGCGACTTATCAGGCCAACGCAAATCATCGCAACATCCTGATTTTGGGTTCTTTGGTCGGTCGCATAAGGCTTTTATGCCGCATCCGACATCATCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                <  2:137160/301‑1 (MQ=14)                                                                                                                                                                                                     gcctgtgtcgacgctggcagcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                <  2:34545/280‑1 (MQ=14)                                                                                                                                                                                                          aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTac                                                                                          >  2:87386/1‑299 (MQ=255)                                                                                                                                                                                                               actggtgcagtgTCTTAACTGGCCTAGGTCAAACAGGACACAAGCCTGATTTAGAGTGATTTTGAAGGTCGGAGAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                     <  2:409062/289‑1 (MQ=14)                                                                                                                                                                                                                         gttCTGATCAGGCCTACGCCAGACAGCGCAATAGTCTGGTGTAGCGGGAGTTTGTAGGTCGGATAAGGCGTTTATGCCGCGTCCGACATGAACGCCTGATGCGACGCTTAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                           <  1:295894/298‑1 (MQ=17)                                                                                                                                                                                                                         ggTCTTAACAGGGCTAGGCCGGACAGGGAAATGGCCTGATTTAGGGGTAGTTTGGAGGTCGGATAAGGCGTGTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                           <  1:180010/299‑1 (MQ=17)                                                                                                                                                                                                                         ggTCTCATCAGGGCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAGGGCGTGTAGGCCGCATCCGACATCAAAGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                           <  1:87386/299‑1 (MQ=17)                                                                                                                                                                                                                                   ggCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATAAACGACTGATGCGACGATTACCGCGTATTATCAGGCCTACGACAGAAAGCGCAATAGCCTGATTTAGCGTGCTTTTGTATGACGGATAAGGAGTGGACGC‑C‑GCATCCTACATCAATGACTGATGCGACGCTTTCAGCGGCTTATCAGGCCTATCTTAACCACTGTTTAATTTTTCACAAACTTACCACCAGTTCTGAAACCCTTcacaacggatagtgtta                                                                >  2:387125/1‑285 (MQ=2)                                                                                                                                                                                                                                       tagcccaatcAGCGCAATACCCTGATTGAGCGTGAGTTTGTAGGTCGGATAATGCGTTTATGCCGCATCCGACATCAAAGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAcc                                                            <  2:27262/292‑1 (MQ=17)                                                                                                                                                                                                                                        aagcgagagagcggtagggggtgagtgagggTGAGTTTGTGGGTCGGGGAAGGCGTTGAGGCCGCAGCGGGCATCAACGGCTGATGCGACGCTTAACGCGTCTTATCAGGCGTACGCGAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAcc                                                            <  1:384845/270‑1 (MQ=17)                                                                                                                                                                                                                                        aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATTCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACCGCCCGTGTTAAACCa                                                           >  2:177691/1‑301 (MQ=255)                                                                                                                                                                                                                                         cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCACACCCTTCCATCTAGTTCCGTTTCAACCat                                                          >  2:320936/1‑300 (MQ=17)                                                                                                                                                                                                                                                  gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGACAAGGCGTTTACGC‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACATCTTCACACCCGGTATCAAACCCTTCAATACCGCTCGGTTCCAACCCATGCacctc                                                  >  1:279583/1‑296 (MQ=14)                                                                                                                                                                                                                                                         tgCCTAATTTAGGGTGATTGTGTGGGTGGGATAAGGCGTTGATGCCGCAGCCGCAATCACCTCCTGATGCGACGCTTAACGCGTTTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCt                                          <  1:368259/300‑1 (MQ=18)                                                                                                                                                                                                                                                          gtcTGACTTCGCGTTAGTTTGTATGTGGGACAAGGCGTTTACGGCGCAGCCGACATTAACGCCGGATGCGACGCTTAACTGGGCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCt                                          <  1:109495/298‑1 (MQ=18)                                                                                                                                                                                                                                                              gATTTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACCTCACTGCCTGTTGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCTTTGGTTAATTTTTCAAACTCTTCACACCCGGTACCCACCCCTTCCCTACACCTCCGGTTCAACCAGTGCAGCCCCTTTTGCTTGtc                                      >  1:84520/1‑299 (MQ=17)                                                                                                                                                                                                                                                              gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGAATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTCGGTTAATTTTTCCAACTCTTCACACCCGGTATCAAACCCTTCCATTCATCTCAGGTTCAACCAAAGCACACCCCTTTGCTTTttt                                     >  1:403291/1‑300 (MQ=17)                                                                                                                                                                                                                                                                           tgtggaggtgcgATAAGGCGTTTATGCCGCCTCCGCAATCAACGCCTGATGCGACGCTTAACGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAc                        <  1:191468/290‑1 (MQ=21)                                                                                                                                                                                                                                                                                      gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTCTCAAACCCTTCCCTACCGCTCCGGtt                                                                 >  1:172097/1‑249 (MQ=255)                                                                                                                                                                                                                                                                                             cGTTTATGCCGCATTCGTCATCAACGCCTGATGCGGCTCTTAACGCGTCTTATCAGGCCTACGCCAGACTGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAg       <  2:128477/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                 tATGCCGCATCCGACAACAACGCCTGATGCGACTCTTAACGCGTCGTATCAGGCCTAAGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcat  <  1:161541/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                             GGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACAT  >  NC_000913/4295978‑4296564

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.