Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,907,165	IS150 (–) +4 bp :: +CA	intergenic (‑105/+58)	rlmA ← / ← cspC	23S rRNA m(1)G745 methyltransferase/transcription antiterminator and regulator of mRNA stability CspC

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1907165 =	0 (0.000)	16 (0.610) +TG	16/486	0.5	100%	intergenic (‑105/+61)	rlmA/cspC	23S rRNA m(1)G745 methyltransferase/transcription antiterminator and regulator of mRNA stability CspC
	?	NC_000913	3720633 =	NA (NA)					noncoding (1/1443 nt)	IS150	repeat region
*	?	NC_000913	= 1907168	0 (0.000)	10 (0.380)	10/486	1.0	100%	intergenic (‑108/+58)	rlmA/cspC	23S rRNA m(1)G745 methyltransferase/transcription antiterminator and regulator of mRNA stability CspC
	?	NC_000913	= 3722075	NA (NA)					noncoding (1443/1443 nt)	IS150	repeat region

CCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGT  >  NC_000913/1906916‑1907401                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             ccGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCa                                                                                                                                                                                                                                             <  2:67655/251‑1 (MQ=255)    tctggatcagcactaacgcctttaggtttgggaattaaacctttttctccatgtttttcatagagggcaacccatgtcctgacctgggttcgggggacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTc                                                                                                                                                                                                                                          <  1:52061‑M1/5‑1 (MQ=255)                      gCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACGGATGTCGCTGGGGAAAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTtgtactgaccccaaaaag                                                                                                                                                                                                                         <  1:178633‑M1/250‑19 (MQ=255)                                    aattaaacctttttctccatgtttttcatagagggcaacccatgtcctgacctgggttcgggggacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTAtt                                                                                                                                                                                                          >  1:176265‑M1/215‑251 (MQ=255)                                                ttctccatgtttttcatagagggcaacccatgtcctgacctgggttcgggggacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGCAGCGTATTCACCCTTc                                                                                                                                                                                                  <  2:121033‑M1/45‑1 (MQ=255)                                                    ccatgtttttcatagagggcaacccatgtcctgacctgggttctggggacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCa                                                                                                                                                                                            >  1:90612‑M1/199‑249 (MQ=255)                                                                  gagggcaacccatgtcctgacctgggttcgggggacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGAt                                                                                                                                                                            >  1:121033‑M1/185‑251 (MQ=255)                                                                   agggcaacccatgtcctgacctgggttcgggggacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATa                                                                                                                                                                           >  1:114629‑M1/184‑251 (MQ=255)                                                                     ggcaacccatgtcctgacctgggttcgggtgacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGc                                                                                                                                                                         >  1:214805‑M1/182‑251 (MQ=255)                                                                        aacccatgtcctgacctgggttcgggggacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTGATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGt                                                                                                                                                                      >  1:160491‑M1/179‑251 (MQ=255)                                                                              tgtcctgacctgggttcgggggaaaccacaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCCCACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGtt                                                                                                                                                                <  2:269377‑M1/79‑1 (MQ=255)                                                                                      cctgggttcggggtacaccaaaacgtgccgagatgatcctgtaaccatcatcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTg                                                                                                                                                        >  2:213289‑M1/165‑251 (MQ=255)                                                                                                    tttCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaatt                                                                                                                                           <  1:275851‑M1/250‑97 (MQ=255)                                                                                                                      gATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtat                                                                                                                         <  2:206073‑M1/250‑115 (MQ=255)                                                                                                                              aGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTtgtactgaccccaaaaagttggacagttaaacacgaggcat                                                                                                                                                                                                  <  1:130359‑M1/169‑42 (MQ=255)                                                                                                                              aGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTtgtactgaccccaaaaagttggacagttaaacacgaggcat                                                                                                                                                                                                  >  2:130359‑M1/1‑128 (MQ=255)                                                                                                                                       atcagttgtgaagtagtgattcacgacttcaaggcgcttttcaaaagggtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTc                                                                                                       >  1:269377‑M1/116‑251 (MQ=255)                                                                                                                                                aaagtagggattcacgacttcaagtcgcttttcaaaagggtattttggctttgacagagtaggggctattccttttcagcgtccaacaaaatgggtgcagtacacaCATTCGTCTAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCGGTTACTGTAACAGCTGCCGGACCTTTTTGGCCGTCCTGAATTTCGAACTCTACGTCCTGACCTTCAGCCAGAGt                                                                                              >  2:163962‑M1/107‑251 (MQ=255)                                                                                                                                                     agtgattcacgacttcaaggcgcttttcaaaagtgtattttggctttgacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTg                                                                                          <  2:146364‑M1/149‑1 (MQ=255)                                                                                                                                                           ggCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaatatattccgtaacagcatccttcagttcgcttat                                                                                   <  1:112403‑M1/251‑153 (MQ=255)                                                                                                                                                                  tACTAAAGGCCGCAAGTGTAACGCGACTCAAGCGGCGGGCCAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgctattcaattggagacagaccttttaatttcaggctaattcttctgctgttgtagtattcaatatattccgtaacagcatccttcagttcgcttatattactg                                                                            >  2:175141‑M1/1‑92 (MQ=255)                                                                                                                                                                      aaaGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAACAGCCGGACGCCAGCAACGAAGATAAAGCATAGTAAAAAGCCTCGCATTtgaactgaccccacaaagttggccagttaaacaagaggcatataggactgcttccgatattcaattggagccaaaccttttaatttcaggctaattcttatgctgttgtagaattcaatatattccgtaacagcatccttcagttcgcttatattactgaaat                                                                        >  1:178922‑M1/1‑88 (MQ=255)                                                                                                                                                                                                   ggcggGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTACAAAGCCTCGCATTtgtactgaccccaacaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaatatattccgtaacagcatccttcagttcgcttatattactgaactcatcaagataaaaacactccgacttta                                             >  2:125570‑M1/1‑59 (MQ=255)                                                                                                                                                                                                    gacatattaggggctattccatttcatcgtccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgctg                                          <  1:244674‑M1/197‑1 (MQ=255)                                                                                                                                                                                                                                  ccaacaaaatgggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCt            <  1:87204‑M1/227‑1 (MQ=255)                                                                                                                                                                                                                                          aTAGTAAAAAGCCTCTCATTtgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaattcttctgctgttgtagtattcaatatattccgtaacagcatccttcagttcgcttatattactgaactcatcaagataaaaacactccgactttaaggttccaaagaaacactccaccacagcattatccagaca     >  2:178305‑M1/1‑20 (MQ=255)                                                                                                                                                                                                                                            gggtgcagtacacaCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTACTGAAGCCAAAACCTTTAGACTCGt  <  2:14362‑M1/237‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             CCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAGTGGACAGGAAAAAGACATGACAACTCCGGCAGGATACTAAAGGCCGCAAGTGTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGT  >  NC_000913/1906916‑1907401

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ 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.