Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,286	C→T	19.7%	intergenic (+492/+150)	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,286	0	C→T	19.7%	13.1 0.0 0.0	51	intergenic (+492/+150)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (5/5):  ref base (19/19):  total (25/26)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCC  >  NC_000913/4295989‑4296570                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                       tatagctctaaagtcTGATTTAGTGTGACGTTGTAGTTTGGATAAGGCGTTTACGCCGCATCCCACATCAACTCCTGATGCGACGCTCGCCGTGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAACCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                                                                                                                     <  1:673736/287‑1 (MQ=11)             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACCCCTGCTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCCGCCAGCGCAATAGCCTGCTTTAGCTTGATTTTGTATTTCGGATAATGCGTTTCCGCCGCCTCCGCCTTCAATCCCTgctccgccc                                                                                                                                                                                                                                                                                           >  1:566485/1‑291 (MQ=11)             gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAAGCCTGATGCGACGATTAACGCGTCTTCTCAGGCCTACGCCAGACCGAGCAATAGCCTGATTTAGCGTGCTTTTGTAGGACGGATAAGACGCTTACGCCGCCTCCGACATCAATGCCTGCTGCCACg                                                                                                                                                                                                                                                                                           >  1:367806/1‑299 (MQ=11)                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGTCCTCCGCCCGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCTTTTACGCCGCATCCGACCTCAATGCCTGATGCGAGGCttgc                                                                                                                                                                                                                                                                                      >  2:605640/1‑298 (MQ=11)                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAGGGCTGTTACGCCCCATCAGACGTCAACgaccggcgcagctcttac                                                                                                                                                                                                                                                                                      >  2:80672/1‑283 (MQ=11)                            gtgtcgtagctttgtagaatagccagcgccgcAGCTCACACCAACCTCTGACGCCATGCTTGCCTCGTCTTTTCAGGCCTACCCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                           <  1:321302/274‑1 (MQ=11)                            aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATAAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTATCGTACTTTTGTAGGTCTGATAAGGAGTTTACGCCGCATCCGACATCAATGCCTTATACGACGCTTGCCGCGTCGTATc                                                                                                                                                                                                                                                                           >  1:429638/1‑300 (MQ=11)                                 gtaggccgctgaggctgttgccgccgcACGCGACAGCAACGCGGGGTGCGACGTTTGCCTCTGCTTATCAGGACTCCGCCAGAGCGGGCAATATCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTAGTCCGCATCAGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                      <  2:373747/279‑1 (MQ=11)                                      tCGGTTAATGCGTTTACGCCGCATCCGAAAGCAACGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                                                                                                                                                                                 <  1:605664/300‑4 (MQ=11)                                               gCGTTTACGCCGCATCCGACATCAACGACTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGCTACGAAGCTTAACGCGTCTTATCCGGCATAAGCCAGACAGCGCAATCGCCTGACTTAGCGTGATTTTGTAGTTCGGTTACGGCGTTTACGCCGCCTACGCCACCATGGCCTGCTGCGCCTCTGGCCTCGTCTTAGCAGACCTaacttaacagtttt                                                                                                                                                                                                                                                       >  1:707188/1‑288 (MQ=11)                                                      cgccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTACGCCGCATACGACATCAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                           >  2:321296/1‑274 (MQ=25)                                                                                                                                                            tAAGGCGTTTACGCCGCATGCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGCTAAGACGTTGATGCCGCATCCGACATCAACACCTTATGCGGCGGTTAACGCGTCTTATCAGGCCTACGCCAGACGCCCAAATAGCCGGATGTAGCGTGAATTTGTAGGTCGGATCAGGCGATA‑‑C‑‑A‑‑C‑CAGCATCCGACATCAATGCCTAGTGGGACGATTGCCGCGTGTTCTACGGCGTATGTTAAcc                                                                                                                                    >  1:739110/1‑300 (MQ=11)                                                                                                                                                             aggtgtttgtggagcgtaacGACATGATCGCCGGATGCGACGTTGGACGCGTCTTAGCAGGCCTACTGCAGCAAGCGCAAAAGCCTGATTTAGCGTGATTTTGTAGGGCGGAGAAGGCGTTGATGCGGCATCCGCCAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                                    <  2:530762/280‑1 (MQ=11)                                                                                                                                                               ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATACGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCCATGCCTGATGCGACGCTTGCCGCGTCTCATGAGGCCTATGTTAACCGTTg                                                                                                                                >  2:257762/1‑301 (MQ=11)                                                                                                                                                                                  ccATCCAGGCCCGATGCGACGCTTGACGCAGCTTATGTGTCCGACGCAAGCCCACGCAAGGGCCTGATTTAGCTTAATTGTGGAAGTCGGATAAGGCGTGTATGCCGCATACTACATAAGCGCTTGAAGCGACGCTTAACGCTGCTGATCAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                                                                                 <  1:582756/296‑1 (MQ=11)                                                                                                                                                                                       ccgtcctggggcgcagctggctgtgtgtcaTCAGGACTAGGCCAGACAGCGCAAAAGGTTGATGTAGACGGAGTGTGTCGGCAGGATAAGGCGTGTATGCCGCAGCCGACAGCAACGCCTCATGCGACGCTTAACGCGTCGTATCAGGCCTACGCCAGACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                                                                                        <  2:429627/272‑1 (MQ=14)                                                                                                                                                                                         cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTCGCGTGATTTTGTCGGTCGGATCAGGCGTTTATGCCGCATCCGCCCTCAACGCCTGATGCTACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACCGCGCCGTCGCCTGATTTATCGTGCTTTTATCGGTCGGATACTGCG‑TT‑‑T‑‑ACGC‑C‑GCCTCCGCCGTCACTGCCTGATGCGCCCCTTCCCTCATCTTCTCCGGCCTATCTTCACCTTTGGTTAATTTTTCAAACTCTTccctcc                                                                                                       >  1:450962/1‑295 (MQ=11)                                                                                                                                                                                             cgaggcgcagctggcggctttttATCAGGCCTACGCCAGAACGCGCAAAGGCCTGGTTTAGCGTGAATGTGTAGGGCGGATAAGGAGTTTATGCCGCATCCGACATCAACGGCTGATGCGACGGTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                                                  <  2:38625/280‑1 (MQ=14)                                                                                                                                                                                                tGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCTTCTTACCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGtcctt                                                                                               >  1:38627/1‑297 (MQ=14)                                                                                                                                                                                                    acgctggccgcgTCTTACAAGGCCGACTCAAGACAGCTAAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGTGTTCATGCCGCTTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                            <  1:648335/292‑1 (MQ=14)                                                                                                                                                                                                     ccctgtcgtcgTTTTCTCAGCTCTCCGCCATAACGCGTAATCGCCTTATTTAGCGTGATTTTGTTGGTCGGATAAGGCGTTTATGCCGCACCCGTCATCAACGCTTGATTCGAAGCTTACCGTGTCTTATCAGCCCTCCGCCCGACAGCTCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGATTAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCTACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCCCGGTATCAAAcc                                                                                           <  1:356066/291‑1 (MQ=255)                                                                                                                                                                                                     cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTATGCCGCATCCGACATCACCGCCTGATGCGCAGCTTCACGCGTCTTCTCCTGCCTACGCCCGACAGCCCCATCGCCTTCTTTAGCGTGATTTTGTATTTCGGATAAGGCG‑GT‑‑T‑‑A‑‑CGCATCATACGACATCTATGCCTTCTGCGCCGCTTACCCCGTCTTATCAGTCCTCTCTTGACCGTGGGGTACTTTTTCATCCTCTTCACACCCGTTCTCAACcc                                                                                           >  2:684066/1‑300 (MQ=2)                                                                                                                                                                                                     cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                           >  1:614966/1‑300 (MQ=255)                                                                                                                                                                                                     cGCTTGCCGCGTCTTACCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCCTCCGACCTCAATGCATGATTCTACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTCGGTTCCTTTTTCCAACTCTTCCTACCCGCTacccaccc                                                                                           >  2:101851/1‑293 (MQ=11)                                                                                                                                                                                                              ttggttATCAGGCCGTCGCAGGACACCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTTGGATCAGGCGTTTATGCCCCATCCGACATCAACGACGGACGCGACGCTTAACGCTTCTGATCAGTCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTAATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                  <  1:263273/296‑1 (MQ=17)                                                                                                                                                                                                              gtgctgaacAGTCCTACGCGATACATCGCAGTGACATGGTGGAGCGTGGTTGTGGAGGTCGCATAATGCGTGTCTGCCGCATCCGCCATCTATGCCGGGTGCGTAGCTTAACGCGTCTTATCAGGTCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATAACGCG‑TT‑‑T‑‑ACGC‑C‑TCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACACTTCACACCCGGTATCAAACCCTTCCATAc                                                                                  <  1:446629/292‑1 (MQ=17)                                                                                                                                                                                                              gggtttaaccggcgcaggCCAGACAGCGAAAAAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCACCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                  <  1:236663/283‑1 (MQ=17)                                                                                                                                                                                                              gagcTTAGCAGGCCTAAGCCATACAGCAGAATAGCCTGATTTAGCGTGAGTTTGTAGGCAGGATAAGGCGTGTATGCCGCAACCGACATGAACGCCTGATGCGACGCGTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                  <  1:460341/297‑1 (MQ=17)                                                                                                                                                                                                              cgtctgatgaggctgtcggcagcacgcgcAATAGGCTGATTTAGCGTTATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCGGATGCGACTCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                  <  2:300214/275‑1 (MQ=17)                                                                                                                                                                                                                   tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCATGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTAATTTTGTAGGTCGGATAAGGCG‑CT‑‑TGCG‑‑C‑C‑GCATCCGACCTCAATGCCTGATGCGACATTTGCCGCGTCTTATCAGCCCTATCTTACCAGTTGGTTTATCTTTCCAGCCCCTCAACCCCGGTATCAGACCATTCCATTCAGCTc                                                                             >  2:108430/1‑300 (MQ=11)                                                                                                                                                                                                                             aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGCTTTCGCGTGATTTTGTAGGTCGGATACGGCG‑TT‑‑TCCG‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACCCCCGATATCAAACCCTTCCCTACCGCTCCGTTTCA‑Aca                                                                  >  1:565694/1‑299 (MQ=14)                                                                                                                                                                                                                             aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATACGGCG‑TT‑‑T‑‑ACGC‑C‑CCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTCTCAGGCCTATCTTACCCGTTGGTTACTTTTTCACACTCTTCACACCCGGTACCCAACCCTTCCATCCACTCCCGGTTCA‑ACCa                                                                 >  2:102714/1‑301 (MQ=14)                                                                                                                                                                                                                                 acggcAGCGAAATCGCCGTGTTTTGGGTGATTTGGGAGGGCGGATAAGGCGTTTAGGCCGAATCCGAAATCAACGCCGGAGGGGACGCTTAACGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGc                                                             <  2:335581/297‑1 (MQ=18)                                                                                                                                                                                                                                       gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTATGTCGGGTACGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGCTGCGCGGCTTGCCGCGTCTTGTCAGGCATATCTTAACCGTTGGTTAATTTTTCCACCTCTTCACACCCGGTATCAAACACTTCCACACAGCACCGGTCCA‑ACCCATGCAccga                                                        >  2:694041/1‑299 (MQ=14)                                                                                                                                                                                                                                        cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTCCGACAAACAGCGCAATAGCCTGAATTAGCGTGATTTTTTAGGTCGAATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCAGCCGACCTCACTGCCTGCTGAGACGATTGCCGCGTCTTATCAGGCCTCTATTCACCGTTGTGTACTTTTTCCCACTCTTCACACCCGATCTACAACCCTGCCAtcacgaccaccttaccccccttccccgcc                                                        >  2:492174/1‑272 (MQ=11)                                                                                                                                                                                                                                             tggcctggtttgacggtaTTTTGTAGGTCGGATAAGTGGGTGAGGCCGCATCCGACACCAACGCCTTATGGGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGc                                                 <  1:257768/284‑1 (MQ=18)                                                                                                                                                                                                                                             tacctcgatgtatcttGAGTTTGTAGTTCGGAGAAGGCTTTTCTTCCCCAGCCGACTGCAACGCCGTATGCTACGCTTAACCCGTCTTATCAGGCTGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTTGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACAGCAATGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGc                                                 <  1:213936/286‑1 (MQ=18)                                                                                                                                                                                                                                               gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGTTTTAGCGTGATGTTGTAGGTCGGGTAAAGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCGGGCCTATCTTAACCGTTGGTTAATTTGTCAAACTCTTCACACCCGGTGTcca                                                                                              >  1:190918/1‑253 (MQ=14)                                                                                                                                                                                                                                                ccTGATTGAGTGTGAGTTTGGAGGTGTGGTAAGGCGTTTAGTCCGCATCCGACATCAAAGCCTGATGCGACGCTTAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCt                                                <  2:126663/299‑1 (MQ=18)                                                                                                                                                                                                                                                   taTTTAGCGTGATTTTGGAGGTCGGAGAATGCGTTTATTCCTCATCCGACATCAACTCCTGATGCTGCGCCTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGTCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGttc                                           <  1:419299/300‑1 (MQ=18)                                                                                                                                                                                                                                                                tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCCCAGGTTCA‑ACCAGTGCAGCCCCTTTTGCTTGTTCTTCTCGATaaaaa                              >  2:306614/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 tggttgggctgATAAGGCGTTTAAGCCGCATCCGACATCACAGCCTGATGCGCCGCTTAACTGGTTGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAcc                             <  2:342325/291‑1 (MQ=21)                                                                                                                                                                                                                                                                           gATAAGGCGTTTATGCCGCATCCGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCAGTCCTACGCCAGACCGCCCAATAGCCTGATTTAGCGTGATTTTGTAGTGCGGATACGGCG‑TTTCT‑‑G‑‑C‑C‑GCATCCGACATCAATGCCTGATTCGCCGCTTGCCGCGTCTTATCACGCCTATCTTAACCGTTGGTTACTGTTTCACACTCTTCACACCCGGTCTCAACTCCTCCCATACACCAC‑TGTTCACACCTGTGCAGCGCCTTTTCCTTGTTCTTCTCTATAAAACCTTCTTCTCCg                   >  2:549795/1‑301 (MQ=17)                                                                                                                                                                                                                                                                                 gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACCGCTCAGGTTCA‑ACCATTTCACCGCCTTTTGCTTTTTCTTCTCGATAAAACCTTCTCCACCTTTTaaa             >  1:587919/1‑300 (MQ=18)                                                                                                                                                                                                                                                                                  cgttcaagcccgtgcgtcaatacaCGCCTGATGCGACGCTTAACGCGTTTTATGTGGCCTAAGACAGACAGCGCAATAGCCTGTTTTAGCGTGCTTGTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAg             <  1:733439/277‑1 (MQ=25)                                                                                                                                                                                                                                                                                    tgtATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑G‑‑ACGC‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAc          <  1:574436/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                    gttacgcccactcagACATCAACGCCTGAAGGGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCAGCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTGCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa           <  2:565671/286‑1 (MQ=25)                                                                                                                                                                                                                                                                                       aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                               <  2:506508/230‑1 (MQ=255)                                                                                                                                                                                                                                                                                       aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                               >  1:506531/1‑230 (MQ=255)                                                                                                                                                                                                                                                                                          ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCGACGCCAGACAGCGCAATGGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCG‑TT‑‑TTCG‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatcat     <  2:184061/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                             cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCACTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTcc  >  1:446622/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                       GACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TT‑‑T‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCA‑ACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCC  >  NC_000913/4295989‑4296570

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.