Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_pae_tpiA	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_pae_tpiA	1978503–1979269	1979278	10–776	41 [0]	[0] 41	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 1978502	0 (0.000)	39 (0.850)	35/274	0.3	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913_3_pae_tpiA	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_pae_tpiA/1979133‑1979420                                                                                                                                                   |                                                                                                                                              aaatggtgctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                             <  2:68266‑M2/3‑1 (MQ=255) aaatgggtctttacacttatctaagatttttcctaaatagacgcaactgtactcgtcactaacagcacacacaacggagaggggctgcgattttcaataatgcgtgatgcacatcaaacaaaacacccaattacttaacatcaatgtct                                                                                                                                             >  1:375102‑M2/147‑147 (MQ=255)  aatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                            >  2:285526‑M2/146‑149 (MQ=255)      ggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                                                                        <  1:122584‑M2/8‑1 (MQ=255)       gtctttacacttatccaagagttttcctaaatcgacgcaactgtactcgtcactacacgaacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                                                                       <  2:158951‑M2/9‑1 (MQ=255)               acttatctaagatttttcctaaatcgacgcaactgtactcgtcactaccagcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattgcttaacataaatgTATAAGTCATACttttg                                                                                                                               <  1:374322‑M2/17‑1 (MQ=255)                   atctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                           <  1:268485‑M2/21‑1 (MQ=255)                             tttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                 <  1:749453‑M2/31‑1 (MQ=255)                                    aaccgacgcaactgtactcggcaccacacggaaatacaacggagggcggcggcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                          <  1:175664‑M2/38‑1 (MQ=255)                                         acgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                     >  1:394158‑M2/107‑149 (MQ=255)                                                tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                              >  1:313591‑M2/100‑149 (MQ=255)                                                     tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                         >  1:96441‑M2/95‑149 (MQ=255)                                                               acgcgcatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAATTCATACTTTTGTTTTGGGTGTATTTCTATCTGTTAAACAGTTATTAGCTAGTCTagcaa                                                                               >  1:272881‑M2/85‑149 (MQ=255)                                                               acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                               >  1:387892‑M2/85‑149 (MQ=255)                                                                       acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                       >  1:744070‑M2/77‑149 (MQ=255)                                                                         aacggaggggcgctgctattttaaataaggagtgaggcagatcacaaaaaacactcaagtacttaacaaaaatgTATAAGTCATTCTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                     <  2:444082‑M2/75‑1 (MQ=255)                                                                           cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                         <  1:102184‑M2/39‑1 (MQ=255)                                                                           cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                         >  2:102184‑M2/73‑111 (MQ=255)                                                                               ggggggctgcgattttcaataatccgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAACAATGAAACAgg                                                               <  1:521499‑M2/81‑1 (MQ=255)                                                                                gggggcggcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                              <  1:479760‑M2/82‑1 (MQ=255)                                                                                    gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                                                          <  2:96441‑M2/86‑1 (MQ=255)                                                                                    gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                                                          <  2:387892‑M2/86‑1 (MQ=255)                                                                                    gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGATAAAAAGTTTTTCGCTACGCTAGCAACCAAAAATGAAACAGGaataa                                                          >  2:728451‑M2/64‑149 (MQ=255)                                                                                     ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                         <  1:568926‑M2/87‑1 (MQ=255)                                                                                         gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAACATGAAAAAGAAATAATAGaa                                                     >  2:344863‑M2/59‑149 (MQ=255)                                                                                          attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                    >  1:220642‑M2/58‑149 (MQ=255)                                                                                                   aatgcgtgattcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                           <  1:640004‑M2/101‑1 (MQ=255)                                                                                                       cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcgc                                       >  1:583726‑M2/45‑149 (MQ=255)                                                                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                  >  1:585281‑M2/40‑149 (MQ=255)                                                                                                               agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaa                               <  1:715579‑M2/113‑1 (MQ=255)                                                                                                                    acaaaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                          <  2:356543‑M2/118‑1 (MQ=255)                                                                                                                      gcaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCCCACAGTCAAAATAACt                        >  2:462709‑M2/30‑149 (MQ=255)                                                                                                                        aaaacaatcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                      <  2:641394‑M2/122‑1 (MQ=255)                                                                                                                             actcaattacttaacataaatgTATAAGTGATACTTTTGTTTTGGCTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGTGATGTTGCGCACAGTACAAATACCTCCCCGTa                 >  2:773472‑M2/23‑149 (MQ=255)                                                                                                                             actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGATGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGTAATGTTGCGCACAGTCAAAATAATTCACCGTa                 >  2:309306‑M2/23‑149 (MQ=255)                                                                                                                                       ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGCATTTAGCTTCTGTAGCGTGCTAAAATGAAACAAGAGTAATTGAAATGGGTTGTTGAGCACAGTCAAAATAACTCACCGTAAATAatcatc       >  2:775800‑M2/13‑149 (MQ=255)                                                                                                                                        taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCt      <  2:744070‑M2/138‑1 (MQ=255)                                                                                                                                         tatataaatgTATAAGTCATACTTTTGTTTTGGGTTTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAACGAAACAAGTATGTTCTAGATGGGATTTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg     <  2:393441‑M2/139‑1 (MQ=255)                                                                                                                                         aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGATACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAACTAATCATCTg     >  1:22857‑M2/11‑149 (MQ=255)                                                                                                                                         aacataaatgTATAAGTCATACTATCGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg     <  1:644581‑M2/139‑1 (MQ=255)                                                                                                                                            ataaatgTATAAGTCATACTGTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCta  <  2:411408‑M2/142‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_pae_tpiA/1979133‑1979420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑