Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,190	A→G	17.8%	intergenic (+396/+246)	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,190	0	A→G	17.8%	15.7 ‑0.8 ‑0.4	49	intergenic (+396/+246)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  new base (7/2):  ref base (23/16):  total (31/18)
Fisher's exact test strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.33e-01

TTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4295894‑4296482                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                           tttAGCTTGATGTTGTAGGTCGGATAAGGCTGTTATGCCGCCTCCGACATCAACGCCTGATGCGACGGTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAGCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc                                                                                                                                                                                                                                                                                                        <  1:382436/301‑1 (MQ=255)          gttttGTAGGTCGTATAAGGCTGTTATTCCGAATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                                <  1:574110/299‑1 (MQ=255)              tgtggggtggATATGGCGTTTAGACCGCAGCAGACATCAACGCCGGATGCGACGCTTCACGCGTCTTATTAGTGCTACGCAAGCCAGCGCAATAGCCTGATTTAGCGTGAGTTTGTAGGTCGGATAAGGCGTTTACGCCGCAGCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                                           <  2:593024/293‑1 (MQ=11)              tGTAGTGCGGATAAGGCGTTTATGCCGCACCCGACATCAACGCCTGATGCGACGCTTACCGCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGGGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                                           <  1:416174/301‑1 (MQ=16)                 aGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCCGACATCGCAATATCCTGATTTATCTTGATTTTGTCGGTCGGATAAGGCGTTTCAGCCGCATCGGCCGTCACCGCCTGATGCGCGGCTGACCGCTTCTTAtacggcgtc                                                                                                                                                                                                                                                                                         >  2:356058/1‑292 (MQ=11)                       gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAg                                                                                                                                                                                                                                                                                  >  1:591260/1‑300 (MQ=25)                        ttAAGGCGTTGAGGCCGCATCCGCCACCACCGCCTTATGCGACTTTTAACTCGTCTTATCAGGCCTAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAGCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGAc                                                                                                                                                                                                                                                                                <  2:152311/300‑1 (MQ=11)                        aaaaGGCGTGTATGTCCCATCCGATATCAACGCTTTATGCGACGCTTATCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGAc                                                                                                                                                                                                                                                                                <  1:68717/299‑1 (MQ=11)                              cGTGTATGCCGCATCCGACACAAAGGCCTGATGCGACGCTTAACGCGTCTGATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCTGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAgcgc                                                                                                                                                                                                                                                                           <  1:107940/300‑1 (MQ=255)                                          ttCCGACATCAACGCAGGAGTCGACGGTGAACGGTTTCTAGCAGGACTACGCCAGACAGCGGAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGAttt                                                                                                                                                                                                                                                              <  2:591235/300‑1 (MQ=11)                                           tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCGCATCCGACATCAAGCCCTGATGCGACGCTTACACCGTCTTATCAGG‑CTTCGCGCAGACAGCGCAAtgagcggatgta                                                                                                                                                                                                                                                             >  1:61122/1‑290 (MQ=11)                                                 aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGTGGTTTTGTAGGTCGGATAAGGCGTTGATGCCGCATCCGCCATCAACGCCTGATGTGACTCTTAACGCTCTTTCTCAGGCCTACGC‑CATCCCCCGCCATAGCCTGCTTTCACGTGa                                                                                                                                                                                                                                                       >  2:673710/1‑301 (MQ=11)                                                 aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATACGTCGTTTATGCCGCATCCGACATCAACGCCTGAGCCGACGCTTAACGCGTCTTATCAGTCCTACCC‑CAACCAgccctttaccctacacttccgct                                                                                                                                                                                                                                                        >  2:236657/1‑279 (MQ=11)                                                 aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGAGCAATAGCCTGATTTAGCGTgt                                                                                                                                                                                                                                                       >  1:99592/1‑300 (MQ=21)                                                                      ttcgctgtcgtatactgcgtACGCCAGACAGCTGAATGTCGTGTGTTAGCGTTATGTGGTCGGTCGGATAAGGCGTTTACGCCGCATCCGAAAGCACGGCCTGATGCGACGCTTGCTGCGTCTTATGAGGTACACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGCTTGGGTAGGTCGGATAAGGTGATTATGCCGCATCCGACATCAGCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                                                                                                                                                                                                  <  1:205305/282‑1 (MQ=11)                                                                           cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGCTAAGACGTTTATGCCGCCTCCGCCATCAAATCCTGATCCGACGCCTAACGCGTCTTATTAGGCCTACGC‑CAGCTAGCGCAAACGACTGATTCATCCTGATTTTGTAGGTAGGATAAGGCTttcaa                                                                                                                                                                                                                             >  1:373757/1‑298 (MQ=11)                                                                            ggtttATCAGGCCCACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATGAGGCGTTTac                                                                                                                                                                                                                             <  1:80676/297‑2 (MQ=17)                                                                                tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCTCGTCTTATCAGGCCTACGC‑CAGACAGCGCACTACCCTGACTTAGCTTGATTTTGTAGGTCGGATAAGCCGTTTACGCCg                                                                                                                                                                                                                         >  1:335587/1‑300 (MQ=255)                                                                                          aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTGGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGCATAAGGCGTCCATTCCCCCTCCGACATCAACGCCTGATGCGGAGCTTAACGCGCCTTATCATGCCAACAC‑CAGACCGCGCAATAGCCTGATTTAGCGTGCTTTTTTAGGTGCGCTAAGCCGATTagtgcgactccgccat                                                                                                                                                                                                               >  2:383928/1‑285 (MQ=11)                                                                                          aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CCGACAGCGCAATAGCCTGCTTTCGCGTGATCTTGTAGGTCGGATAAGGCGTTTCAGTCGCACCCGAccta                                                                                                                                                                                                              >  1:234195/1‑298 (MQ=255)                                                                                          aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTGTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCTGATGCTACGATTAACGCGTCTTCTCAGGCCTACGC‑CAGACAGCCCACTCACCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTACGCCGCATCCGACAt                                                                                                                                                                                                               >  1:342332/1‑300 (MQ=11)                                                                                          aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAACGTGATTTTGTAGGTCGGCTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTCACGAGTCTGATCAGGCCTACGC‑CAGACAGCGGTATAGCCTGATTTCACGTAATTTTTTCTGTCGACTAAGGGGTTTACGCGGACTCCGACAt                                                                                                                                                                                                               >  1:643245/1‑300 (MQ=11)                                                                                                tatagctctaaagtcTGATTTAGTGTGACGTTGTAGTTTGGATAAGGCGTTTACGCCGCATCCCACATCAACTCCTGATGCGACGCTCGCCGTGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAACCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                        <  1:673736/287‑1 (MQ=11)                                                                                                    gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTCATTTACCGTGATTTTCTCGGTCGGATACTGCGTTTACGCCCCATCCGACATCCACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTAAGC‑CATACACCGCCATATCCTGCTTTCCCGTTATTTTGTAGGTCGTCTACGCCTTTTCTCCCCCATCCGCCATctattcccgt                                                                                                                                                                                                     >  2:587180/1‑291 (MQ=11)                                                                                                     cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACACCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTCTGCCGCATCCGTCATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGC‑CAGCCCGCGCACTAGCCTGAGTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTCCGCCCCATTCGTCAtgcattccagttt                                                                                                                                                                                                   >  2:582732/1‑289 (MQ=11)                                                                                                            gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACCCCTGCTGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CCGCCAGCGCAATAGCCTGCTTTAGCTTGATTTTGTATTTCGGATAATGCGTTTCCGCCGCCTCCGCCTTCAATCCCTgctccgccc                                                                                                                                                                                              >  1:566485/1‑291 (MQ=11)                                                                                                            gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATAAGGCGTTTATGCCCCATCCGCCCTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAACAGCCTTATTTAGCGTGAGTTTGTATGTCGGATAAGCCGCTTACGcaacctacgaactcaatcccctattccgcc                                                                                                                                                                                              >  2:648309/1‑270 (MQ=11)                                                                                                            gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAAGCCTGATGCGACGATTAACGCGTCTTCTCAGGCCTACGC‑CAGACCGAGCAATAGCCTGATTTAGCGTGCTTTTGTAGGACGGATAAGACGCTTACGCCGCCTCCGACATCAATGCCTGCTGCCACg                                                                                                                                                                                              >  1:367806/1‑299 (MQ=11)                                                                                                                gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGTCCTCCGC‑CCGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCTTTTACGCCGCATCCGACCTCAATGCCTGATGCGAGGCttgc                                                                                                                                                                                         >  2:605640/1‑298 (MQ=11)                                                                                                                gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAGGGCTGTTACGCCCCATCAGACGTCAACgaccggcgcagctcttac                                                                                                                                                                                         >  2:80672/1‑283 (MQ=11)                                                                                                                           gtgtcgtagctttgtagaatagccagcgccgcAGCTCACACCAACCTCTGACGCCATGCTTGCCTCGTCTTTTCAGGCCTACCCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                              <  1:321302/274‑1 (MQ=11)                                                                                                                           aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATAAGGCCTACGC‑CAGCCAGCGCAATAGCCTGATTTATCGTACTTTTGTAGGTCTGATAAGGAGTTTACGCCGCATCCGACATCAATGCCTTATACGACGCTTGCCGCGTCGTATc                                                                                                                                                                              >  1:429638/1‑300 (MQ=11)                                                                                                                                gtaggccgctgaggctgttgccgccgcACGCGACAGCAACGCGGGGTGCGACGTTTGCCTCTGCTTATCAGGACTCCGCCAGAGCGGGCAATATCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTAGTCCGCATCAGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                         <  2:373747/279‑1 (MQ=11)                                                                                                                                     tCGGTTAATGCGTTTACGCCGCATCCGAAAGCAACGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                                                                                    <  1:605664/300‑4 (MQ=11)                                                                                                                                      cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAAt                                                                                                                                                                                                                                                                        >  1:12770/1‑199 (MQ=34)                                                                                                                                      cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAAt                                                                                                                                                                                                                                                                        <  2:12768/199‑1 (MQ=34)                                                                                                                                              gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAgg                                                                                                                                                                                                                                    >  2:136328/1‑227 (MQ=18)                                                                                                                                              gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAgg                                                                                                                                                                                                                                    <  1:136334/227‑1 (MQ=33)                                                                                                                                              gCGTTTACGCCGCATCCGACATCAACGACTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGCTACGAAGCTTAACGCGTCTTATCCGGCATAAGC‑CAGACAGCGCAATCGCCTGACTTAGCGTGATTTTGTAGTTCGGTTACGGCGTTTACGCCGCCTACGCCACCATGGCCTGCTGCGCCTCTGGCCTCGTCTTAGCAGACCTaacttaacagtttt                                                                                                                                                          >  1:707188/1‑288 (MQ=11)                                                                                                                                                     cgccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTACGCCGCATACGACATCAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                              >  2:321296/1‑274 (MQ=25)                                                                                                                                                                                                                                                           tAAGGCGTTTACGCCGCATGCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGC‑AAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGCTAAGACGTTGATGCCGCATCCGACATCAACACCTTATGCGGCGGTTAACGCGTCTTATCAGGCCTACGCCAGACGCCCAAATAGCCGGATGTAGCGTGAATTTGTAGGTCGGATCAGGCGATACA‑‑C‑CAGCATCCGACATCAATGCCTAGTGGGACGATTGCCGCGTGTTCTACGGCGTATGTTAAcc                                           >  1:739110/1‑300 (MQ=11)                                                                                                                                                                                                                                                            aggtgtttgtggagcgtaacGACATGATCGCCGGATGCGACGTTGGACGCGTCTTAGCAGGCCTACTG‑CAGCAAGCGCAAAAGCCTGATTTAGCGTGATTTTGTAGGGCGGAGAAGGCGTTGATGCGGCATCCGCCAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                           <  2:530762/280‑1 (MQ=11)                                                                                                                                                                                                                                                              ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATACGGCG‑TTTACGC‑C‑GCATCCGACATCCATGCCTGATGCGACGCTTGCCGCGTCTCATGAGGCCTATGTTAACCGTTg                                       >  2:257762/1‑301 (MQ=11)                                                                                                                                                                                                                                                                                 ccATCCAGGCCCGATGCGACGCTTGACGCAGCTTATGTGTCCGACGC‑AAGCCCACGCAAGGGCCTGATTTAGCTTAATTGTGGAAGTCGGATAAGGCGTGTATGCCGCATACTACATAAGCGCTTGAAGCGACGCTTAACGCTGCTGATCAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                        <  1:582756/296‑1 (MQ=11)                                                                                                                                                                                                                                                                                        cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGC‑CAGACAGCGCACTAGCCTGATTTCGCGTGATTTTGTCGGTCGGATCAGGCGTTTATGCCGCATCCGCCCTCAACGCCTGATGCTACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACCGCGCCGTCGCCTGATTTATCGTGCTTTTATCGGTCGGATACTGCG‑TTTACGC‑C‑GCCTCCGCCGTCACTGCCTGATGCGCCCCTTCCCTCATCTTCTCCGGCCTATCTTCACCTTTGGTTAATTTTTCAAACTCTTccctcc              >  1:450962/1‑295 (MQ=11)                                                                                                                                                                                                                                                                                               tGCGACGCTTGCCGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCTTCTTACCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGtcctt      >  1:38627/1‑297 (MQ=14)                                                                                                                                                                                                                                                                                                    cGCTTGCCGCGTCTTATCAGGCCTACGC‑CAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTATGCCGCATCCGACATCACCGCCTGATGCGCAGCTTCACGCGTCTTCTCCTGCCTACGCCCGACAGCCCCATCGCCTTCTTTAGCGTGATTTTGTATTTCGGATAAGGCG‑GTTA‑‑CGCATCATACGACATCTATGCCTTCTGCGCCGCTTACCCCGTCTTATCAGTCCTCTCTTGACCGTGGGGTACTTTTTCATCCTCTTCACACCCGTTCTCAACcc  >  2:684066/1‑300 (MQ=2)                                                                                                                                                                                                                                                                                                    cGCTTGCCGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc  >  1:614966/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                    cGCTTGCCGCGTCTTACCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCG‑TTTACGC‑C‑GCCTCCGACCTCAATGCATGATTCTACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTCGGTTCCTTTTTCCAACTCTTCCTACCCGCTacccaccc  >  2:101851/1‑293 (MQ=11)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                           TTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGC‑CAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4295894‑4296482

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.