Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,169,293	A→G	6.8%	E69E (GAA→GAG)	tufA →	protein chain elongation factor EF‑Tu

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_007779	4,169,293	0	A→G	6.8%	8.6 0.0 ‑0.1	88	E69E (GAA→GAG)	tufA	protein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand):  new base (3/3):  ref base (46/36):  total (49/39)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.55e-01

TGGCTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGA  >  NC_007779/4169193‑4169393                                                                                                     |                                                                                                     tGGCTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGaa                                                                                                      <  1:6175404/101‑1 (MQ=32) tGGCTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGGGGTATCACCATCAACACTTCTCACGTTGaa                                                                                                      >  1:6772318/1‑101 (MQ=17)  ggCTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAAt                                                                                                     <  1:2786408/101‑1 (MQ=32)  ggCTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAAt                                                                                                     >  1:8454788/1‑101 (MQ=32)  ggCTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAAt                                                                                                     <  1:5152701/101‑1 (MQ=32)    cTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATAc                                                                                                   >  1:13883199/1‑101 (MQ=32)    cTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATAc                                                                                                   <  1:6510177/101‑1 (MQ=32)     tAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACg                                                                                                  >  1:3293640/1‑101 (MQ=32)     tAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACg                                                                                                  <  1:4680416/101‑1 (MQ=32)     tAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACg                                                                                                  <  1:1632541/101‑1 (MQ=32)      aaaaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGa                                                                                                 >  1:12952525/1‑101 (MQ=32)       aaaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGac                                                                                                >  1:12125292/1‑101 (MQ=32)       aaaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGac                                                                                                <  1:2640184/101‑1 (MQ=32)        aaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCTCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATATGaca                                                                                               >  1:11343514/1‑101 (MQ=12)        aaCCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGaca                                                                                               >  1:12729776/1‑101 (MQ=32)          ccTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcc                                                                                             >  1:7807751/1‑101 (MQ=32)           cTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAccc                                                                                            >  1:3179660/1‑101 (MQ=32)           cTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAccc                                                                                            >  1:1991133/1‑101 (MQ=32)            tACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcccc                                                                                           <  1:2326048/101‑1 (MQ=32)            tACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcccc                                                                                           <  1:8969490/101‑1 (MQ=32)            tACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcccc                                                                                           >  1:1396018/1‑101 (MQ=32)            tACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcccc                                                                                           <  1:4060122/101‑1 (MQ=32)             aCGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCg                                                                                          <  1:191818/101‑1 (MQ=32)               ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                        >  1:14283031/1‑101 (MQ=33)               ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                        >  1:13210753/1‑101 (MQ=33)               ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                        >  1:9186230/1‑101 (MQ=32)               ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                        >  1:12511769/1‑101 (MQ=32)               ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                        <  1:9117250/101‑1 (MQ=32)                gcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAcc                                                                                       <  1:6682759/101‑1 (MQ=32)                gcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAcc                                                                                       <  1:11493450/101‑1 (MQ=32)                gcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAcc                                                                                       >  1:11430739/1‑101 (MQ=33)                gcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAcc                                                                                       >  1:1763466/1‑101 (MQ=32)                 cggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAccc                                                                                      >  1:12426084/1‑101 (MQ=32)                  ggcgctgctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                     >  1:13920680/4‑101 (MQ=17)                  ggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                     >  1:13279962/1‑101 (MQ=32)                  ggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                     >  1:7064689/1‑101 (MQ=32)                  ggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                     <  1:6519725/101‑1 (MQ=32)                  ggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                     <  1:10146955/101‑1 (MQ=32)                   gTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGt                                                                                    >  1:10220643/1‑101 (MQ=32)                   gTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGt                                                                                    <  1:8324732/101‑1 (MQ=32)                   gTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGt                                                                                    <  1:5799949/101‑1 (MQ=32)                   gTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGt                                                                                    >  1:762001/1‑101 (MQ=32)                   gTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGt                                                                                    >  1:12587459/1‑101 (MQ=32)                    tgctgcTCGTGCATTCGACCAGATCGATAGCGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTc                                                                                   >  1:2229942/1‑101 (MQ=17)                      ctgctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCAc                                                                                 >  1:13845162/1‑101 (MQ=32)                       tgctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACt                                                                                >  1:11576055/1‑101 (MQ=32)                        gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                               >  1:5022443/1‑101 (MQ=32)                        gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                               <  1:1911560/101‑1 (MQ=32)                        gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                               >  1:11117897/1‑101 (MQ=32)                         ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                              <  1:230775/101‑1 (MQ=32)                           cGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGc                                                                            >  1:4832651/1‑101 (MQ=32)                            gTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGca                                                                           <  1:13135643/101‑1 (MQ=31)                            gTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGca                                                                           >  1:3639293/1‑101 (MQ=31)                                                                                     acgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTgctgct                  >  1:1698867/9‑101 (MQ=2)                                                                                     acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                  >  1:7475059/1‑101 (MQ=11)                                                                                     acacTTCTCACGTTGAACACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                  <  1:10111412/101‑1 (MQ=11)                                                                                      cgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTc                 <  1:12111192/94‑1 (MQ=2)                                                                                      cgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTc                 <  1:5831399/94‑1 (MQ=2)                                                                                       acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                >  1:385698/1‑101 (MQ=32)                                                                                       acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                <  1:2833706/101‑1 (MQ=32)                                                                                       acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                <  1:8830491/101‑1 (MQ=32)                                                                                       acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                >  1:11999868/1‑101 (MQ=32)                                                                                       acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                >  1:14050496/1‑101 (MQ=32)                                                                                       acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                <  1:14223804/101‑1 (MQ=32)                                                                                        cTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAg               >  1:12875481/1‑101 (MQ=32)                                                                                        cTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAg               <  1:1342930/101‑1 (MQ=32)                                                                                         ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa              >  1:7866586/1‑101 (MQ=32)                                                                                         ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa              <  1:13869726/101‑1 (MQ=32)                                                                                         ctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGa              >  1:6129257/5‑101 (MQ=2)                                                                                            tcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATgg           >  1:14267230/1‑101 (MQ=32)                                                                                            tcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATgg           <  1:10521368/101‑1 (MQ=32)                                                                                            tcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATgg           <  1:9296784/101‑1 (MQ=32)                                                                                             cACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGa          >  1:235828/1‑101 (MQ=2)                                                                                             cACGTTGAGTACAAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGa          <  1:8182838/101‑1 (MQ=2)                                                                                             cACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGa          <  1:5474552/101‑1 (MQ=32)                                                                                             cACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGa          >  1:12565821/1‑101 (MQ=32)                                                                                              aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc         <  1:838724/101‑1 (MQ=32)                                                                                              aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc         <  1:6937773/101‑1 (MQ=32)                                                                                               cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg        >  1:9241137/1‑101 (MQ=32)                                                                                                gTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACgg       >  1:7226557/1‑101 (MQ=32)                                                                                                gTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACgg       <  1:12380567/101‑1 (MQ=32)                                                                                                 ttGAATACGACCCCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc      <  1:660654/101‑1 (MQ=17)                                                                                                 ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc      >  1:2028469/1‑101 (MQ=32)                                                                                                 ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGCGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc      >  1:12264688/1‑101 (MQ=17)                                                                                                  tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg     <  1:5982064/101‑1 (MQ=33)                                                                                                   gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc    >  1:7305586/1‑101 (MQ=32)                                                                                                   gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc    <  1:11541244/101‑1 (MQ=33)                                                                                                   gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc    <  1:2095219/101‑1 (MQ=33)                                                                                                   gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc    >  1:9225183/1‑101 (MQ=32)                                                                                                   gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc    >  1:3425338/1‑101 (MQ=32)                                                                                                    aaTACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGcgcg   <  1:6157057/101‑1 (MQ=33)                                                                                                    aaTACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGcgcg   >  1:216080/1‑101 (MQ=32)                                                                                                     aTACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGa  >  1:12727271/1‑101 (MQ=32)                                                                                                     |                                                                                                     TGGCTAAAACCTACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGA  >  NC_007779/4169193‑4169393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ 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.