Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,812,597	C→T	9.4%	intergenic (‑38/‑134)	yicR ← / → dfp	UPF0758 family protein/fused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,812,597	0	C→T	9.4%	2.3 ‑0.4 ‑0.0	49	intergenic (‑38/‑134)	yicR/dfp	UPF0758 family protein/fused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  new base (3/2):  ref base (20/23):  total (24/25)
Fisher's exact test strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.57e-01

CGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCACG‑CCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA  >  NC_000913/3812297‑3812888                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                    ctcgcCAGTCAAGCAACCCCCTTTAACTGGTAAAATTT‑GTCATGTCCAATTCCATGAACGCCACTAAATTGTTCATAGTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCGTCAGCATCGCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGTTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGcc                                                                                                                                                                                                                                                                                                     <  2:383912/299‑1 (MQ=255)   cgacAGTTCAACAATCCCTTTTAACTGGGCAAATTTCGCCAGT‑CAGATTCCAGGAACGACACTCAATTGTTTATATTAAGAGGTTAACAAGCCATAAAGACAGCCGAAATGCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGCTACGCAGAAAGAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCGGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCg                                                                                                                                                                                                                                                                                                    <  2:392778/297‑1 (MQ=255)           cAGCAATCCCCTTTAACTGGGCAAATTTCGCCACG‑CCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGACATGTTTTTCACCTTCAACTCCT‑TTGTGGTGCTCGCAGCCTGTCATAGCCGGTTCGTTaa                                                                                                                                                                                                                                                                                           >  2:35360/1‑301 (MQ=255)                aTCCCCTTTAACTGGGCAAATTTCGCCACG‑CCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCACACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTCCTCGATTCCTGTCTTAGCCGTTTCGTTAAATCGa                                                                                                                                                                                                                                                                                      >  1:420581/1‑301 (MQ=255)                   cccTTTAACTGGGCAAATTTCGCCACG‑CCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGCCGCTAATACCAAACGTCAGCATTTTTTCGCCCGGAATCAACAGATGTGAATTGTTTTTCAACTTCACTTCCT‑TTGTGGTCCTGCCATCCTGTCCTAGCCTTTTCGTTAAATCTACgg                                                                                                                                                                                                                                                                                   >  1:40290/1‑301 (MQ=255)                   cccTTTAACTGGGCAAATTTCGCCACG‑CCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGAAATCCGTTAAGGCGCTAATACCAAACTTCAGCCTTTTTTCGCGCTGCATCAACAGCTGTGACTTGTTTTTAACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGACGTTTCGTTTAATCGACgg                                                                                                                                                                                                                                                                                   >  2:234721/1‑301 (MQ=255)                   cccTTTAACTGGGCAAATTTCGCCACG‑CCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCCCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTATTGTGCTCGCATCCTCTCATAGCCGTTTTGTTAAATccccgg                                                                                                                                                                                                                                                                                   >  2:229905/1‑296 (MQ=255)                                 aaaATTCACCACG‑CCAATTCCATGAACGCCACTAAAATGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATcc                                                                                                                                                                                                                                                                     <  2:155990/301‑1 (MQ=255)                                                              ccACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCTCTAATACCAAACTTCAGCATTTTTTCGCGCGGACTCAACAGCTGTTAATTGTTTTTCACCTTCACCTCCT‑TTTTGGTGCTCGCATCCTGTCATATCCGTTTCTTCACCTTGACGGCCCTTTTTAACTCATGGAAAGCGCCTCGCCAAGTTCATCTCa                                                                                                                                                                                                                                         >  2:237489/1‑301 (MQ=255)                                                              ccACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGTCATCCACAGCTGTGCATTGTTTTTCACCTTCACCTCCT‑TTGTTGTGCTCTCCTCCTGTCATAGCCGTTTCGTTAAATCTCCGGCCCGTTTTCACTCCTTGAACACGCCTCCCCAATTGAATCTCa                                                                                                                                                                                                                                         >  1:318384/1‑301 (MQ=255)                                                              ccACGAAAGTGTTCATATTCATAGGTTAACAAGCCATAAAGAGGCCCGAAATTATCCAGCAGCTCGTGTGCCAGGGTTAAGACATCTGTACCGCGCGTTCCGGTACGAAGAAATAGCGCCATCAGCTCGACATCCGTAAATGCGCTAATACCAAACTTCAGCATTTTGTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAAtctc                                                                                                                                                                                                                                          <  2:420590/300‑1 (MQ=255)                                                                aatAAAGTGTTCATATTCAGAGGTTAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCCGTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAgg                                                                                                                                                                                                                                       <  1:229900/299‑1 (MQ=255)                                                                       ttttctgctacacagGTTAACAAGCCAAAAAGAGCGCCTCGATTCGCCCGCATCTCTTTTGCCAGGGTTAATACATCGTGACAGAGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACAGCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGGGCTGCATCAACAGCTGTGAATTGTTTTTCACCTTCACATCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGc                                                                                                                                                                                                                                <  1:167753/288‑1 (MQ=255)                                                                                    tggTTAAAAAGCCATAAAAAGTGCCTAAATTCTTGTGCATCTCTTTTGCCAGGTTTAATACAGCTTTACCGCTCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTAAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTAGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAg                                                                                                                                                                                                                   <  2:318392/300‑1 (MQ=255)                                                                                       ttAACAAGCCATAAAGAGAGCCGAAATTCTCCAACATCTCTTTTGGCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTGAAGGCGCTAATACCAAACTTCAGCATTGTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGGGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTAt                                                                                                                                                                                                                <  1:237484/301‑1 (MQ=255)                                                                                            tatactagcaagggtccGAAATTTCTCACCATATCTCTTGCAAGGTTTAATACCTCCTGACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACAGCCGTCAAGTCGCTAATCCCAAACTTCATCATTTTTTCGCGCGTCACTAACATCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTTTGGTGCGCGCATCCTCGCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGa                                                                                                                                                                                                           <  2:178406/286‑1 (MQ=255)                                                                                              gCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTCAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCAACAACAGCTGTGAATTGTTTTTCACCTTCACCCCCT‑TTGTGGTGCTCGCATCCTGTAATAGCCGTTTCGTTAAATCGAACGCCAGTTTTCAATCCTGGACAGTCCCTCTCCAAGTGAATCTCAGGTTAAGGccagaaaccacacatattgtgtta                                                                                                                                                                                                         >  2:34882/1‑278 (MQ=255)                                                                                               ccATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACACCAAAAGTATTTTGATaa                                                                                                                                                                                                        >  1:363452/1‑301 (MQ=255)                                                                                                        gagCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAACGTGAATCTCAGGTGACGGCCCGACCACCAACGTATTTTGa                                                                                                                                                                                                           >  1:178405/1‑289 (MQ=255)                                                                                                                       atcagttctTTTGCCAGGGTGAATACATCTTGCCTGCGTGTTCCGGGACGCAGAAAGAGCGCCCGCAGCTCGACATCCGTTAAGGCGCTAAGAACAAACTTCAGCATTGTGTAGCGCGGCATCAACAGCTGTGAATTGTGTTTAACATTCACCTCCT‑TTGTGGTGCGAGCATCCTGGCATAGCCGTGTCGTTAAAGCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGTCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGTTGTCAccc                                                                                                                                                                                <  1:365017/295‑1 (MQ=255)                                                                                                                              ctTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTAAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACATACT‑TTGTGGTGCTCGAATCCTGTCATAGCCGTGTCGTTAAATCGCCGCCAAGTTTTCAATCCTGGAAAGGGCCTCGCAAAGGGAATCTACGGTGAAGCAACGACAACAACAGAATTGTGGTACACACGCCAACTTATGGAGTCAccaacaataa                                                                                                                                                                         >  2:365678/1‑293 (MQ=255)                                                                                                                                  tGCCAGGGTTAATACATCTGTACCGCGCGGTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGaaaaa                                                                                                                                                                     <  1:234716/301‑1 (MQ=255)                                                                                                                                      cgggTTAATACCGTTTTATCGACCGTTCCGGTACGCAGAAATAGCGCCAGCCTCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAatcat                                                                                                                                                                 <  2:40290/300‑1 (MQ=255)                                                                                                                                       gggTTAATACATCTTTACCGCGCGTTCGGGTACGCAGCAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCACACTTCAGCATTTTTTAGCGAGTCATCAACAGCTGTGAATTGTTTTTCACCTTCCCCTCCT‑TTGTGGTGCTCGCATCATGTCATAGCTGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGCACGTCCCTCGCCACGTGAATCTCCGTGCAAGGCCCGACAACAACAGTATTCTGATAAAATCCCCCACTTCTGGTGTCACCCtcactgaacactcata                                                                                                                                                                >  1:234459/1‑286 (MQ=255)                                                                                                                                        gttgaaTACATCTTTACCGCCCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACGTCAGCATTTTTTCTCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGa                                                                                                                                                                                                           <  1:198122/253‑1 (MQ=255)                                                                                                                                          ttAATACACCTGTACCGCGAGGTCCGGTCAGCAGACACAGCTCCAGCTACTCGACCGACGTGAAGGCGCGCATACCCAACTTTAGCATTTTTGCGCGCCGCATCAACAACTGTTAAGTTTTTTTTCCCTTCACCTACT‑TTGTGGTGCTCGCAGACTTTCATAGCCGTTTCGTTAAAACGACGGCCAGTTGTCAATCCTGGAAAGCGTCTCGCAAAGTGAATCTCAGGTGAAGTCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATg                                                                                                                                                             <  2:234464/301‑1 (MQ=255)                                                                                                                                           tAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCAGTTTCGTTAAATCGACTGCCCGTTTTCAATCCTGGAACTCGCCTCGCACACTGCATTCTATGTCAAGCCACCACACCCAAAGTATTTTTATCAAACACCCAACTCCTTgtgtacacctcccggacaaaccctaataa                                                                                                                                                            >  2:158107/1‑276 (MQ=255)                                                                                                                                                cATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCCTCCTGTCATAGCCGTTTCGTTACATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCCGGTGACGGCACGACAACAACAGTATTGTGATAAACTCGCCACCTTCTGGTGTCACCCTCCAGGAAAAatcacactgagc                                                                                                                                                          >  1:384102/1‑290 (MQ=255)                                                                                                                                                  tCTTTTCCCCGCGTTCCGCTACGCAGCAATAGCTGCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCATCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGc                                                                                                                                                     <  1:35360/301‑1 (MQ=255)                                                                                                                                                  tCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTTCTCGCATCCTGTCATAGCCGTTTCGTTAACTCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCGCGACAACGAAAGTCTTGTGATAAAATCGCAAACTTCTGGGGTCAACCTACAGGAAAAATAACCATGAGACTgag                                                                                                                                                     >  1:394190/1‑299 (MQ=255)                                                                                                                                                                                           cTCGAAATCCGGTACGGCGCGAAGACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCGTTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTAAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGGATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATtgctg                                                                                                            <  1:303302/301‑1 (MQ=255)                                                                                                                                                                                           aacagcATCCGTTAAGGCGCTAATACCAAACTTCACCAGTTTGTCGCGCGGGATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCTGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATtgctg                                                                                                            <  1:83701/296‑1 (MQ=255)                                                                                                                                                                                                cATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCTGTTTCGTTAAATCGACGGCCCGTTTTCAATCCTGGAAAGAGCCTCGCAAAGTGAATCCCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACCTCTGGTGTCACCCTCCAGGCAAACTCATCATGAGCCTTGCCGGTAAAACAATCCTTCACGGCGTTAACGGCGGCAGtgctc                                                                                                            >  2:83701/1‑295 (MQ=255)                                                                                                                                                                                                   ccGTTAAGGCGCTAATACTAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAACGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTCAAACAATCGTTCTCGCCGTTAGCGGGGGTGTTGCTGCCTATaaa                                                                                                    >  1:165437/1‑301 (MQ=255)                                                                                                                                                                                                      ttAAGTAGCTAAGACCAATCTTCCGCAGTTCTTCGCGCCGCATCAACAGCTTTTAATTGGTTTTCTCATTCATCGCCT‑TCGTGGCACTCGCCTCCTGTCATAGCCGTTTCGTTAAATCTACGGCCAGGTTTCAATCCTGGAAAGCGTCTCGCAAAGTTAATCTCAGGTTAATGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCGGGTGTCACCCTACAGGAAAAATCATCATGCGCCTGGCCGGTAAAAAAATCGTTATCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAcc                                                                                                 <  2:394199/301‑1 (MQ=255)                                                                                                                                                                                                           gcgcTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTAGTTAAATCGACGGCCAGGTTTCAATCCTGGACAGCGCCTCGCCAAGTGAATCTCAGGTGACGGCCAGACAACAAACGTGTTGTGATAAAATCGCCACCTTCTTGTGTCACCCTACCGTAAAAAGCATCATGAGCATTGCCGGTACAAAACACTTTCTCGACGGTACAGGCGGTATATCTTCC‑ATACAAccccctga                                                                                           >  2:235372/1‑297 (MQ=255)                                                                                                                                                                                                             gcTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAatc                                                                                                                                                                   >  1:347742/1‑228 (MQ=255)                                                                                                                                                                                                             gcTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAatc                                                                                                                                                                   <  2:347751/228‑1 (MQ=255)                                                                                                                                                                                                             cgcaattgcAAACTTCAGCATTTGTTCGCGCGGCAGCAACAGCTTGTAATTGTTTTGCACCGTCACCCCCG‑TTGTGGTTCTAGCAGCCTGTAATAGCCGTTTCTTTAAATCGACGGCCAGTTTGCAATCCGTGAAATCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAAc                                                                                          <  2:363461/293‑1 (MQ=255)                                                                                                                                                                                                              cTACTACCTCACTTCATCCTTGTTGTGCGCTCCATCAACCGCTGTTAATTGTGTTCCACCTTCACCTCCT‑TTTTGTTGATCGCATCTTGTCCTACCCGTTTCGTTAAATGGACTGCCAGTTCTCAATTCTGGAAATCTCCTCGCAAAGTGAATCTCAGGTTAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGTTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTATCGGCGGTATTGCTGCCTATAAAACCCCTGAAc                                                                                          <  2:384111/300‑1 (MQ=255)                                                                                                                                                                                                                                        cctcgGCATTAAGAGCGGTTAATTTTTTTTCCCCTTCACCTCCTGTTGT‑GTGCTCGCATCCGGTGAGAGCCGTTTCGTTAAATCGACGTCCAGTTGTCAAGCCGGGAAAGCGCCTCGCAAATTGAATCGCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGTTGTCACCCTACAGGAAAAACCATCATGAGCCTTGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCgggg                                                               <  1:235367/298‑1 (MQ=255)                                                                                                                                                                                                                                         ccgggCTTCATCAGCTGTGAATTGTTTTCTACCTTCCCCTCTT‑TTTGTGTGCACGCATCCCGTCAGCGCCGTTTCGTTAAATCGCCGGCCAGGTTTCAAGCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGc                                                              <  1:125563/298‑1 (MQ=255)                                                                                                                                                                                                                                             gCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTTCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAc                                                                                                  <  2:69581/261‑1 (MQ=255)                                                                                                                                                                                                                                                tCAACAGCTGTGAATTGTTTTTCACCTGCACCTCCT‑TTGGTGTGCTCGTATCCTGACATAGCCGTTTCGTTAAATCGACGGCCAGTTCTCAATCCTGGAAAGCGACTCTCAAAGTGAATTTCAGGTGAAGGCACTACAACAAAAGTATTGTGATAAAATCGCCAACTTCGGGTGTCACCCTACAGGAAAAAGTATCATGAGCCTGGCCGGTAAAAAACTCGTTCTCTGCGTTAGCTGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCTGGGCCGACGTc                                                       <  2:149042/301‑1 (MQ=255)                                                                                                                                                                                                                                                      gCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGTCAAGCGCCTCGCCAAGTGAATCTCAGGTGAAGGCACGACAAAAAAAGTATTGTGAGAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAACATCAGCATGAACCTGGCAGGTAAAAAAATCGTTCTCGGCGTTAGCCGCCGTATTCCTTCCTCTACAACCCCTGCACTGCTGCTGCGttttccgcatccagtcgccgtcgtcccgtta                                                 >  1:228162/1‑273 (MQ=255)                                                                                                                                                                                                                                                                gTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCCCGGGGCCGACGCCCGCGTAGCCATTAcc                                        >  1:21572/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          ttcacctCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGAGAATCTCAGGTGCAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTTTCACCCTCCAGGAAAAATCATCATGAGCCTAGCCGGTAAAAAAATCGTTCTCGGAGTTAGCGGCGGTATTGCGGCATATAACGACCCTGACATGTTGAGTCGTTTGAGCGCTCGCGGGGCCGACGGCCTGTTAGCCATGACCGTAGCggaaga                             >  2:148811/1‑297 (MQ=255)                                                                                                                                                                                                                                                                                 cct‑TTGTGGTGCTCGCATCCTGTCATAGCTGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGacaac                                                                                                                                                                                                                        >  2:242745/1‑107 (MQ=255)                                                                                                                                                                                                                                                                                 cct‑TTGTGGTGCTCGCATCCTGTCATAGCTGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGacaac                                                                                                                                                                                                                        <  1:242740/107‑1 (MQ=255)                                                                                                                                                                                                                                                                                                      tGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCAGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACACCTGGACTTGTACGTCGTTTGCGGGCTCGCGGGGCCGACGTCCGCATAGCCATGACCTACTCGGCCAACGCCTTTATCAGCCCACtacacttgcg  >  2:388219/1‑292 (MQ=255)                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                    CGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCACG‑CCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCGCGCGTTCCGGTACGCAGAAATAGCGCCAGCAGCTCGACATCCGTTAAGGCGCTAATACCAAACTTCAGCATTTTTTCGCGCGGCATCAACAGCTGTGAATTGTTTTTCACCTTCACCTCCT‑TTGTGGTGCTCGCATCCTGTCATAGCCGTTTCGTTAAATCGACGGCCAGTTTTCAATCCTGGAAAGCGCCTCGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA  >  NC_000913/3812297‑3812888

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ 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.