Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NZ_LN832404	1,982,569	(ATTA)2→1	86.7%	intergenic (‑536/‑117)	flhD ← / → uspC	flagellar transcriptional regulator FlhD/universal stress protein UspC

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_LN832404	= 1982568	4 (0.050)	26 (0.320)	22/560	NT	86.7%	intergenic (‑535/‑121)	flhD/uspC	flagellar transcriptional regulator FlhD/universal stress protein UspC
	?	NZ_LN832404	1982573 =	4 (0.050)					intergenic (‑540/‑116)	flhD/uspC	flagellar transcriptional regulator FlhD/universal stress protein UspC

TTCAATCTGTTAAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATATTTACCTGTTTGACCGAGTTGGGATTGCGTCGTTTCTCCATTAGGAGTAAAGCTTTAATGTCACCTGAAGTTCACAGAATAAGGAACAGGTTATGAGCTATAGCAATATTCTTGTCGCTGTTGCGGTTACACCGGAAAGTCAGCAACTGCTGGCAAAAGCGGTATCTATCGCCAGGCCAGTAAAGGGACACATCAGTTTAATTACTCTCGCTTCCGACCCGGAAATGTACAATCAATTAGCTGCGCCGATGCTTGAAGATTTACGC  >  NZ_LN832404/1982435‑1982862                                                                                                                                      |                                                                                                                                                                                                                                                                                                       aagaaaggcaacaatcaacttaatgtTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaa                                                                                                                                       >  1:93541‑M1/13‑134 (MQ=255)     aaggcaacaatcaacttaATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtc                                                                                                                                   <  1:97088‑M1/287‑166 (MQ=255)       ggcaacaatcaactTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcag                                                                                                                                 <  2:33188‑M1/289‑168 (MQ=255)        gcaacaaTCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATATTTACCTGTTTGACCGAGTTGGGATTGCGTCGTTTCTCCATTAGGAGTAAAGCTTTAATGTCACCTGAAGTTCACAGCATAAGGAACAGGTTATGAGCTATAGCAATATTCTTGTCGCTGTTGCGGTTACACCGGAAAGt                                                                                                                                >  2:69285/6‑295 (MQ=255)         caacaaTCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATATTTACCTGTTTGACCGAGTTGGGATTGCGTCGTTTCTCCATTAGGAGTAAAGCTTTAATGTCACCTGAGGTTCACAGAATAAGGAACAGGTTATGAGCTATAGCAATATTCTTGTCGCTGTTGCGGTTACACCGGAAAGTc                                                                                                                               >  1:480415/5‑295 (MQ=255)           acaatcAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaac                                                                                                                             >  1:538613‑M1/3‑124 (MQ=255)                cAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctg                                                                                                                        <  2:538613‑M1/295‑177 (MQ=255)                        tgttaatgAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagc                                                                                                                >  1:90766‑M1/1‑111 (MQ=255)                             atgAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgtgtcgtttctccattaggagtaaagctttaatgtcacctgaagttaacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtat                                                                                                           <  2:512761‑M1/295‑190 (MQ=255)                                    gTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgc                                                                                                    >  1:283119‑M1/1‑99 (MQ=255)                                         aCATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgatggcaaaagcggtatctatcgccaggc                                                                                               >  1:675328‑M1/1‑94 (MQ=255)                                             atatTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttaccccggaaagtcagcacctgctggcaaaagcggtatctatcacaaggccaga                                                                                           >  2:449229‑M1/1‑90 (MQ=255)                                              tatTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcacctgctggcaaaagcggtatctaacgccaggccctga                                                                                          >  2:134128‑M1/1‑89 (MQ=255)                                              tatTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacacaggaaagtcagcaactgacggcaacagcggtatctatcgccaggccagta                                                                                          >  2:595217‑M1/1‑89 (MQ=255)                                                  aaCCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagg                                                                                      >  1:650757‑M1/1‑85 (MQ=255)                                                  aaCCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagg                                                                                      <  1:294145‑M1/295‑211 (MQ=255)                                                        ttGTTTAAAGCGGAGAATCATGAAACTAGTTACTTATAATTACATTTGAAAGTATATGAGCACCGGTAGTGATTGATTAtctttactaataatcagacgaatatttccctgtttgaccgagttgtgattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacaca                                                                                <  2:182033‑M1/295‑217 (MQ=255)                                                             tAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaacagcggtatctatcgcccggccagtaaagggacacatcagt                                                                           >  1:734056‑M1/1‑74 (MQ=255)                                                                             taacatagtttcttataatTACAATTGTAATTATATGTTCACCGGAGTTCATTGAATAtcttttctaataatcagactaatatttacctgtttgacctagttgggattgggtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacacatcagtttgattactctcgctt                                                           <  2:675328‑M1/287‑238 (MQ=255)                                                                              acacTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacacatcagtttaattactctcgcttc                                                          >  1:322333‑M1/1‑57 (MQ=255)                                                                               cacTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacacatcagtttaattactctcgcttcc                                                         >  1:644599‑M1/1‑56 (MQ=255)                                                                                acTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggaaaaagcggtatctatcgcccggccagtaaagggacacatcggtttaattactatcgcttccg                                                        >  2:225184‑M1/1‑55 (MQ=255)                                                                                                          aTTATATGATCACCGGGAGTGATTGATTAtctttactaataagcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgatctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacacatcagtttaattactctcgcttccgacccggaaatgtacaatcaattagct                              <  1:483054‑M1/295‑267 (MQ=255)                                                                                                              tatGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacacatcagtttaattactctcgcttccgacccggaaatgtacaatcaattagctgcgc                          >  1:595854‑M1/1‑25 (MQ=255)                                                                                                              tatGATCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggcaagtaaagggacccatcagtttaattactctcgcttccgacccggaaatgtacaatcaattagctgcgc                          >  2:211414‑M1/1‑25 (MQ=255)                                                                                                               atGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATATTTACCTGTTTGACCGAGTTGGGATTGCGTCGTTTCTCCATTAGGAGTAAAGCTTTAATGTCACCTGAAGTTCACAGAATAAGGAACAGGTTATGAGCTATAGCAATATTCTTGTCGCTGTTGCGGTTACACCGGAAAGTCAGCAACTGCTGGCAAAAGCGGTATCTATCGCCAGGCCAGTAAAGGGACACATCAGTTTAATTACTCTCGCTTCCGACCCGGAAATGTACAATCAATTAGCTg                         >  1:132216/1‑295 (MQ=255)                                                                                                                   tCACCGGGAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacacatcagtttaattactctcgcttccgacccggaaatgtacaatcaattagctgcgccgatg                     <  1:661416‑M1/295‑276 (MQ=255)                                                                                                                     aCCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATATTTACCTGTTTGACCGAGTTGGGATTGCGTCGTTTCTCCATTAGGAGTAAAGCTTTAATGTCACCTGAAGTTCACAGAATAAGGAACAGGTTATGAGCTATAGCAATATTCTTGTCGCTGTTGCGGTTACACCGGAAAGTCAGCAACTGCTGGCAAAAGCGGTATCTATCGCCAGGCCAGTAAAGGGCCACATCAGTTTAATTACTCTCGCTTCCGACCCGGAAACGTACAATCACTTAGCTGccccta                   >  2:88660/1‑293 (MQ=255)                                                                                                                         ggAGTGATTGATTAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctattgccaggccagtaaagggacacatcagtttaattactctcgcttccgacccggaaatgtacaatcaattagctgcgccgctgcttgaa               >  1:249473‑M1/1‑14 (MQ=255)                                                                                                                            gTGATTGATTAtttttattaataaataagctaatatttacctgtttgaccgagttgggattgggtagtttctccattaagagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacacatcagtttaattactctcgcttccgacccggaaatgtacaatcaattagctgcgccgatgcttgaagat            <  2:73072‑M1/295‑285 (MQ=255)                                                                                                                                 tgattAtctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtacatatcgccaggcacgtaaagggacacatcagtttaattactctcgcttcagacccggaacggcacaacaaatacgctgcgcggttcgttgaagagcaaca       >  2:314650‑M1/1‑6 (MQ=255)                                                                                                                                      atctttactaataatcagactaatatttacctgtttgaccgagttgggattgcgtcgtttctccattaggagtaaagctttaatgtcacctgaagttcacagaataaggaacaggttatgagctatagcaatattcttgtcgctgttgcggttacaccggaaagtcagcaactgctggcaaaagcggtatctatcgccaggccagtaaagggacacatcagtttaattactctcgcttccgacccggaaatgtacaatcaattagctgcgccgatgcttgaagatttacgcagcg  <  2:322333‑M1/295‑295 (MQ=255)                                                                                                                                      |                                                                                                                                                                                                                                                                                                       TTCAATCTGTTAAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACACTATTTTCTTATAATTACATTTGAAATTATATGATCACCGGGAGTGATTGATTAATTATCTTTACTAATAATCAGACTAATATTTACCTGTTTGACCGAGTTGGGATTGCGTCGTTTCTCCATTAGGAGTAAAGCTTTAATGTCACCTGAAGTTCACAGAATAAGGAACAGGTTATGAGCTATAGCAATATTCTTGTCGCTGTTGCGGTTACACCGGAAAGTCAGCAACTGCTGGCAAAAGCGGTATCTATCGCCAGGCCAGTAAAGGGACACATCAGTTTAATTACTCTCGCTTCCGACCCGGAAATGTACAATCAATTAGCTGCGCCGATGCTTGAAGATTTACGC  >  NZ_LN832404/1982435‑1982862

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