Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	2,406,557	(TGCGGCAGC)1→2	coding (85/939 nt)	lrhA ←	DNA‑binding transcriptional dual regulator LrhA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2406541 =	1 (0.020)	15 (0.360)	13/476	1.0	88.3%	coding (101/939 nt)	lrhA	DNA‑binding transcriptional dual regulator LrhA
	?	NC_000913	= 2406549	3 (0.070)					coding (93/939 nt)	lrhA	DNA‑binding transcriptional dual regulator LrhA

GATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGT  >  NC_000913/2406316‑2406799                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                           gattcatctgaagcaccgatagttaatacgccctgaagattactgaacattaatgatgagcaggcctcatcattaaaacgcaggattttcctggcgtaaccaagaagttgaatgccatgttcagttaacagtttgttgcgaccgtgacgagcgaacagttctttcccaacgagttggtccagacgctgcatttgctgacttacggcggactgagtacgacacacagcggcagctgcggcagctgcggcagc                                                                                                                                                                                                                                                    >  2:123964‑M1/235‑251 (MQ=255) gATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGGTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagc                                                                                                                                                                                                                                           >  2:123964‑M2/1‑234 (MQ=255)               accgatagttaatacgccctgaagattcctgaacattaatgatgagcaggcctcatcattaaaacgcaggattttcctggcgtaaccaagaagttgaatgccatgttcagttaacagtttgttgcgaccgtgacgagcgaacagttctttcccaacgagttgttccagacgctgcatttgctgacttacggcggactgagtacgacacacagcggcagctGCGGCAGCTGCGGCCGcaaaaggtttc                                                                                                                                                                                                                                          <  2:52711‑M1/27‑6 (MQ=255)               aCCGATAGTTAATACGCCCTGAAGATTCCTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggccgcaaaaggtttc                                                                                                                                                                                                                                 <  2:52711‑M2/247‑28 (MQ=255)                                cctgaagattactgaacattaatgatgagcaggcctcatcattaaaacgctggattttcctggcgtaaccaagaagttgaatgccatgttcagtgaacagtttgttgcgaccgtgacgagcgaacagttcttgcccaacgagttgttccagacgctgcatttgctgacttacggcggactgagtacgacacacagcggcagctGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAAt                                                                                                                                                                                                                     >  2:261114‑M1/204‑251 (MQ=255)                                ccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCTGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTGAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTGCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagcaaaagtgttcagatcggcaacagcaacaaat                                                                                                                                                                                                            >  2:261114‑M2/1‑203 (MQ=255)                                       attactgaacattaatgatgagcaggcctcatcattaaaacgcaggattttcctggcgtaaccaagaagttgaatgccatgttcagttaacagtttgttgcgaccgtgacgagcgaacagttctttcccaacgagttgttccagacgctgcatttgctgacttacggcggactgagtacgccacacagcggcaactGCGGCAGCTGCGGCAGCAAAAGtgttaagaacggccacagcaaccaattgtctcc                                                                                                                                                                                                              <  2:309833‑M1/55‑12 (MQ=255)                                       aTTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGCCACACA‑‑‑‑‑‑‑‑‑GCGGCAACTgcggcagctgcggcagcaaaagtgttaagaacggccacagcaaccaattgtctcc                                                                                                                                                                                                     <  2:309833‑M2/251‑56 (MQ=255)                                                                                gcatgattttcctggcttaaccaagaagttgaatgccattttcagttatcaggatgttgcgaccgtgacgagcgaacagtgctttcccaacgagttgttccagacgcttcatgtgctgacttacggcggactgagtacgacacacagcggcagctGCGTCAGCTGCGGCAGCAAAAGTGTTCAGATCGGAAACAGCAACAAAGTTTCTTAGCATATCGAGGTCGAGGTTAATTATCGGACGATTTGCACtt                                                                                                                                                                     >  1:38345‑M1/156‑251 (MQ=255)                                                                                gCATGATTTTCCTGGCTTAACCAAGAAGTTGAATGCCATTTTCAGTTATCAGGATGTTGCGACCGTGACGAGCGAACAGTGCTTTCCCAACGAGTTGTTCCAGACGCTTCATGTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcgtcagctgcggcagcaaaagtgttcagatcggaaacagcaacaaagtttcttagcatatcgaggtcgaggttaattatcggacgatttgcactt                                                                                                                                                            >  1:38345‑M2/1‑155 (MQ=255)                                                                                   ggattttcctggcgtaaccaagaagttgaatgccatgttcagttaacagtttgttgcgaccgtgacgagcgaacagttctttcccaacgagttgttccagacgctgcatttgctgacttacggcggactgagtacgacaaacagcggcagctGCGGGAGCTGCGGCAGAAAAAGAATTCAGATCGGCAAAAACCACCAATGTTCTCAGCAGATCGGGGTCGAGGTTAATTATAGGACAATTTGCACTTATc                                                                                                                                                                  <  2:244407‑M1/99‑1 (MQ=255)                                                                                   ggATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGACAAACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcgggagctgcggcagaaaaagaattcagatcggcaaaaaccaccaatgttctcagcagatcggggtcgaggttaattataggacaatttgcacttatc                                                                                                                                                         <  2:244407‑M2/251‑100 (MQ=255)                                                                                                                  gccatgttcagttaccagtttgttgcgaccgtgccgagcgaacagttctttcccaacgcgttgttccagacgctgcatttgctgacttacggcggactgagtacgacacacagcggcagctGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGCGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCCTACTg                                                                                                                                   <  1:264781‑M1/130‑1 (MQ=255)                                                                                                                  gCCATGTTCAGTTACCAGTTTGTTGCGACCGTGCCGAGCGAACAGTTCTTTCCCAACGCGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagcaaaagtgttcagatcggcaacagcaacaaatgttctcagcagatcgaggtcgcggttaattatcggacgatttgcacttatcatatattatcacttactggcggctcctactg                                                                                                                          <  1:264781‑M2/251‑131 (MQ=255)                                                                                                                       gTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGTACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGCATTACGGCGGACTGGGTACGACAAACAGCGGCCGCTGCGGCAGCTGCGGCAGCCAACGGGTTCAGATCGGCAACAGCCACAAATATTCTCAGCAGGGCGAGGTCGAGGTTAAATATAGGGCGGTTTGCTCTTCTCTTATATTATCACTTACTGGCGGCTccttctgcgctg                                                                                                                              >  2:210186/1‑240 (MQ=255)                                                                                                                            gttaacagtttgttgcgaccgtgacgagcgaacagttctttcccaacgagttgttccagacgctgcatttgctgacttacggcggaatgagtacgacccacagcggcagctGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTa                                                                                                                          <  1:411112‑M1/139‑1 (MQ=255)                                                                                                                            gTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGAATGAGTACGACCCACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagcaaaagtgttcagatcggcaacagcaacaaatgttctcagcagatcgaggtcgaggttaattatcggacgatttgcacttatcatatattatcacttactggcggctcatactgagctggtta                                                                                                                 <  1:411112‑M2/250‑140 (MQ=255)                                                                                                                                                        cgaacagttctttcccaacgagttgttccagacgctgcatttgctgacttacggcggactgagtacgacacacagcggcagctGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCCACAGCCAAAACTGTTCTCCGCAGATCGCGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTTTGCACACCCCAACAAGCAATTcc                                                                                             <  2:123629‑M1/168‑1 (MQ=255)                                                                                                                                                        cGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagcaaaagtgttcagatcggccacagccaaaactgttctccgcagatcgcggtcgaggttaattatcggacgatttgcacttatcatatattatcacttactggcggctcatactgagctggttaatgctttgcacaccccaacaagcaattcc                                                                                    <  2:123629‑M2/251‑169 (MQ=255)                                                                                                                                                                          agagttgttccagacgctgcatttgctgccttacggcggactgcgtacgacacacaacggccgctGCGGCAGCTGCGGCAGCAAAAGGGTTCAGAACGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGGTAATTATCGGCCGATTTGCACTTATCATATATTCTCACTTACTGGCGGCTCCTACTGAGCTGGTTAAAGCCGTGCACACACAAACAAACACATCCATTTTTAAATTGCCTccc                                                                           <  2:324973‑M1/186‑1 (MQ=255)                                                                                                                                                                          agaGTTGTTCCAGACGCTGCATTTGCTGCCTTACGGCGGACTGCGTACGACACACA‑‑‑‑‑‑‑‑‑ACGGCCGCTgcggcagctgcggcagcaaaagggttcagaacggcaacagcaacaaatgttctcagcagatcgaggtcgagggtaattatcggccgatttgcacttatcatatattctcacttactggcggctcctactgagctggttaaagccgtgcacacacaaacaaacacatccatttttaaattgcctccc                                                                  <  2:324973‑M2/250‑187 (MQ=255)                                                                                                                                                                                              atttgctgacttacggcggactgagtacgacacacagcggcagctGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAAcc                                                       >  2:214304‑M1/46‑251 (MQ=255)                                                                                                                                                                                              atttgctgacttacggcggactgagtacgacacacagcggcagctGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAAcc                                                       >  2:297015‑M1/46‑251 (MQ=255)                                                                                                                                                                                              aTTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagcaaaagtgttcagatcggcaacagcaacaaatgttctcagcagatcgaggtcgaggttaattatcggacgatttgcacttatcatatattatcacttactggcggctcatactgagctggttaatgctgtgcacacacaaacaagcaattccatttgtaatgtgcctccctggcagtttcatccggaacc                                              >  2:214304‑M2/1‑45 (MQ=255)                                                                                                                                                                                              aTTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagcaaaagtgttcagatcggcaacagcaacaaatgttctcagcagatcgaggtcgaggttaattatcggacgatttgcacttatcatatattatcacttactggcggctcatactgagctggttaatgctgtgcacacacaaacaagcaattccatttgtaatgtgcctccctggcagtttcatccggaacc                                              >  2:297015‑M2/1‑45 (MQ=255)                                                                                                                                                                                                  gctgacttacggcggactgagtacgacacacagcggcagctGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAAGTCCATTTGTAATGTGCCTCCCTGGCAGTTTCAGCCGGAACCCGGa                                                   >  2:191597‑M1/42‑251 (MQ=255)                                                                                                                                                                                                  gCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagcaaaagtgttcagatcggcaacagcaacaaatgttctcagcagatcgaggtcgaggttaattatcggacgatttgcacttatcatatattatcacttactggcggctcatactgagctggttaatgctgtgcacacacaaacaagcaagtccatttgtaatgtgcctccctggcagtttcagccggaacccgga                                          >  2:191597‑M2/1‑41 (MQ=255)                                                                                                                                                                                                                            cacacagcggcagctGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGc                         >  1:143694‑M1/16‑251 (MQ=255)                                                                                                                                                                                                                            cacaca‑‑‑‑‑‑‑‑‑GCGGCAGCTgcggcagctgcggcagcaaaagtgttcagatcggcaacagcaacaaatgttctcagcagatcgaggtcgaggttaattatcggacgatttgcacttatcatatattatcacttactggcggctcatactgagctggttaatgctgtgcacacacaaacaagcaattccatttgtaatgtgcctccctggcagtttcatccggaacccggacgaaagtaaaaatgcatatgagttgc                >  1:143694‑M2/1‑15 (MQ=255)                                                                                                                                                                                                                                          tGCGGCAGCTGCGGCAGCAAAAGAGTTCAGATCGGCAACAGCCACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCAATTATCATATATTATCACTTAATGTCGGCTCATACTGAGATGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGGTTCATCCGGAACCCGGACGAAAGTAAAAATACATATGAGTTGCACTAAAAAAGCGAc           <  1:70157‑M1/250‑1 (MQ=255)                                                                                                                                                                                                                                               cAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGcac                                                                                                               >  2:406849/1‑146 (MQ=255)                                                                                                                                                                                                                                               cAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGcac                                                                                                               <  1:406849/146‑1 (MQ=255)                                                                                                                                                                                                                                                   tgcggcagctgcggcagcaaaagagttcagatcggcaacagccacaaatgttctcagcagatcgaggtcgaggttaattatcggacgatttgcaattatcatatattatcacttaatgtcggctcatactgagatggttaatgctgtgcacacacaaacaagcaattccatttgtaatgtgcctccctggcaggttcatccggaacccggacgaaagtaaaaatacatatgagttgcactaaaaaagcgac  <  1:70157‑M2/251‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                           GATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAGTACGACACACA‑‑‑‑‑‑‑‑‑GCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGT  >  NC_000913/2406316‑2406799

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