Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	CP009273	3,915,133	Δ4 bp	100%	coding (276‑279/816 nt)	atpB ←	F0 sector of membrane‑bound ATP synthase, subunit a

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 3915132	1 (0.050)	34 (1.810)	32/260	NT	97.1%	coding (280/816 nt)	atpB	F0 sector of membrane‑bound ATP synthase, subunit a
	?	CP009273	3915137 =	1 (0.050)					coding (275/816 nt)	atpB	F0 sector of membrane‑bound ATP synthase, subunit a

TGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGACCCAGTACATGTTCAGCAATGTACGGCAGCAGGTCGATAGGCAGTAAATCCATCAGGTTCATCAGGAATACCCAGACGAAGATCGTCAGGGCCAGCGGAGCAATCAGCTTGCTTTTGCCATGGTACATGTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAA  >  CP009273/3914889‑3915356                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                            ttcgtgaagccgccgatgcctttcattttgatgctgtagaacagaatcaggataaatacgcccagtgccatagacagcgttacgttcacgtccgcagacggaaccacacgcagtgcaggcagacccagtacatgttcagcaatgtacggcagcaggtcgataggcagtaaatccataaggttcatcaggaatacccagacgaagtacgtcagggccagcggagcaataagcttgattttgacagggtaTCATACACGCTCCCATTAAAAAAGcc                                                                                                                                                                                                    >  2:384123‑M2/249‑274 (MQ=255)                                     tagaacagaatcaggataaatacgcccagtgccatagacagcgttacgttcacgtccgcagacggaaccacacgcagtgcaggcagacccagtacatgttcagcaatgtacggcagcaggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCttt                                                                                                                                                                                                                         >  1:323255‑M2/213‑217 (MQ=255)                                             aatcaggataaatacgcccagtgccatagacagcgttacgttcacgtccgcagacggaaccacacgcagtgcaggcagacccagtacatgttcagcaatgtacggcagcaggtcgataggcagtaaatccatcaggttcatcaggactacccagccgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCCATCGCGGTCTGAAACTtacctggccacccgc                                                                                                                                                        >  1:103121‑M2/205‑267 (MQ=255)                                                                   gccttagacatcgttacgttcacgtcagcagccggaaccacacgcagtgcaggcagacccattacatgttcagcaatgttcggcagcaggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGgctacgcta                                                                                                                                 <  1:384129‑M2/93‑1 (MQ=255)                                                                                             cgagcacgagaacacacgaagagcagtcagacccagtacatggtcagaaatggacggcagcaggtcgataggcaggaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGAc                                                                                                        <  2:103117‑M2/118‑1 (MQ=255)                                                                                                                   gcaggcagacccagtacatgttcagcaatgtacggcagcaggtcgatagtcagtaaatccatcaggttcatctggaatacacagacgaagatcgtcagggccagcggagcaatcagcttgattttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATaaaa                                                                                                                       >  1:25093‑M2/135‑237 (MQ=255)                                                                                                                           acccagtacatgttcagcaatgtacggcagcaggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGa                                                                                                                                                                                                  >  2:322742‑M2/127‑154 (MQ=255)                                                                                                                               agtacatgttcagcaatgtacggcagtaggtcgagaggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGNCACACCGCNGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCGAGAAGAACATGGAGACAAtatta                                                                     >  1:274539‑M2/123‑274 (MQ=255)                                                                                                                               agtacatgttcagcaatgtacggcagcaggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCAc                                                                                                                                                                                              <  1:322748‑M2/32‑1 (MQ=255)                                                                                                                                                 tacggcagcaggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACt                                                                                                                                                                       >  2:349970‑M2/105‑159 (MQ=255)                                                                                                                                                    ggcagcaggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCttt                                                                                                                                                                                                                         >  2:272213‑M2/102‑106 (MQ=255)                                                                                                                                                     gcagcaggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTAcc                                                                                                                                                                   <  1:349976‑M2/59‑1 (MQ=255)                                                                                                                                                        gcaggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACg                                                                                                                                                                                                                     <  1:272218‑M2/9‑1 (MQ=255)                                                                                                                                                         caggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCa                                                                                                                                                                                                               >  1:363603‑M2/97‑111 (MQ=255)                                                                                                                                                           ggtcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCGAGAAAAACATGGAGTCAata                                                                        >  1:170631‑M2/95‑244 (MQ=255)                                                                                                                                                             tcgataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTaa                                                                                                                                                                                                           <  2:363597‑M2/19‑1 (MQ=255)                                                                                                                                                               gataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCGAGAAAAACATGGAGTCAATattga                                                                    <  2:170626‑M2/154‑1 (MQ=255)                                                                                                                                                                ataggcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtat                                                                                                                                                                                                                             >  1:91504‑M2/90‑90 (MQ=255)                                                                                                                                                                    gcagtaaatccatcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCttt                                                                                                                                                                                                                         <  2:91500‑M2/5‑1 (MQ=255)                                                                                                                                                                                tcaggttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCGAGAAGAACATGGAGTCAATattgattg                                                                 >  2:116238‑M2/74‑230 (MQ=255)                                                                                                                                                                                  aggttcatcaggaatacccagactaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTg                                                                                                                                          >  1:154879‑M2/72‑155 (MQ=255)                                                                                                                                                                                    gttcatcaggaatacccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCa                                                             <  1:116242‑M2/161‑1 (MQ=255)                                                                                                                                                                                      tcatcaggcatacccagactaagagcgtcagaggcagcggagcaagcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGACTTATTGgcta                                                                                                                                      <  2:154874‑M2/88‑1 (MQ=255)                                                                                                                                                                                                  cccagacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCttt                                                                                                                                                                                                                         >  1:334079‑M2/56‑60 (MQ=255)                                                                                                                                                                                                   agagaaaaagatggtcagggccaggngagcaatcactttgcttttgccttggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGGTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGaa  <  2:274534‑M2/220‑1 (MQ=255)                                                                                                                                                                                                      gacgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACg                                                                                                                                                                                                                     <  2:334073‑M2/9‑1 (MQ=255)                                                                                                                                                                                                       acgaagatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGaac                                                                                                               >  1:381839‑M2/51‑161 (MQ=255)                                                                                                                                                                                                           agatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAacaaca                                                                                                           <  2:381833‑M2/115‑1 (MQ=255)                                                                                                                                                                                                            gatcgtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGc                                                                                                                                                                                                                    >  2:294744‑M2/46‑55 (MQ=255)                                                                                                                                                                                                                gtcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTAcc                                                                                                                                                                                                                <  1:294750‑M2/14‑1 (MQ=255)                                                                                                                                                                                                                 tcagggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTg                                                                                                                                                                            >  2:85530‑M2/41‑90 (MQ=255)                                                                                                                                                                                                                     ggccagcggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAAc                                                                                                                                                                        <  1:85533‑M2/54‑1 (MQ=255)                                                                                                                                                                                                                           gCAATCAGCTTGCTTTTGCCATGGTA‑‑‑‑TCTCTCATCCTACCATCAACAAAGCCGATCAc                                                                                                                                                                                              <  1:278130/58‑1 (MQ=255)                                                                                                                                                                                                                           cggagcaatcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGa                                                                                                                                                                                                  >  2:278125‑M2/31‑58 (MQ=255)                                                                                                                                                                                                                                   tcagcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCttt                                                                                                                                             >  2:52101‑M2/23‑103 (MQ=255)                                                                                                                                                                                                                                     agcttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCaccacc                                                                                             >  2:125378‑M2/21‑149 (MQ=255)                                                                                                                                                                                                                                       cttgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTgg                                                                                                                                         <  1:52104‑M2/85‑1 (MQ=255)                                                                                                                                                                                                                                         tgcttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCgaga                                                                                         <  1:125382‑M2/133‑1 (MQ=255)                                                                                                                                                                                                                                            ttttgccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGgctac                                                                                                                                     >  1:8752‑M2/14‑102 (MQ=255)                                                                                                                                                                                                                                                gccatggtaTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGgctacgcta                                                                                                                                 <  2:8752‑M2/93‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                            TGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGACCCAGTACATGTTCAGCAATGTACGGCAGCAGGTCGATAGGCAGTAAATCCATCAGGTTCATCAGGAATACCCAGACGAAGATCGTCAGGGCCAGCGGAGCAATCAGCTTGCTTTTGCCATGGTACATGTCTTTCACGCTACCATTAACAAAGCCGATCACCAGCTCAATCGCGGTCTGAAACTTACCTGGCACACCGCTGGTCGCCTTTTTGGCTACGCTACGGAATAAAACCAGGAACAACAGACCCAGCACCACCGAGAAGAACATGGAGTCAATATTGATTGTCCAGAAGGTGGCTGGGGGGTTTTGTGGATCCACCAGCGAGAATGTACGCAGGTCCAGCTGAA  >  CP009273/3914889‑3915356

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