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	1979278	776	53 [0]	[0] 53	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)	51 (0.900)	41/268	0.4	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

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913_3_pae_tpiA/1979133‑1979421                                                                                                                                                   |                                                                                                                                               aactgggtctttacacttatctaagagttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                              <  1:928445‑M2/3‑1 (MQ=255) aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                              <  2:799073‑M2/3‑1 (MQ=255)     gtgtctttacacttaactaagattgttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                                                                          <  2:669569‑M2/7‑1 (MQ=255)      tgtctttacacttatctaagttttttcataaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                                                                         <  2:7570‑M2/8‑1 (MQ=255)                      taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                                                                    >  1:340126‑M2/126‑138 (MQ=255)                      taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                                                                    <  2:340126‑M2/13‑1 (MQ=255)                       aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                        >  1:226955‑M2/125‑149 (MQ=255)                       aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                        >  1:12267‑M2/125‑149 (MQ=255)                       aagagtttttccaaattgacgctactgtactcgtcactacacgcacatacaacggagggtggctgcgattttcaataatgcgtgatgcatatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                        <  2:1028274‑M2/25‑1 (MQ=255)                          ttttttcctaagtcgacgcaactgtacgagtcactacacgcacatacaacggtggggggctgcgattttcagtaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATGCTTGGGTTTTGGGTGTa                                                                                                                     <  2:125294‑M2/28‑1 (MQ=255)                                  gaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                             >  1:741864‑M2/114‑149 (MQ=255)                                      tcgactcaactctacgcgtcactacacgcacatacaacggaggggtgctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                         <  1:11491‑M2/40‑1 (MQ=255)                                                         cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                      <  1:333229‑M2/59‑1 (MQ=255)                                                          actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTGTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                                                                     >  1:118230‑M2/90‑149 (MQ=255)                                                           ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                    >  1:838020‑M2/89‑149 (MQ=255)                                                              cacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                                                                 >  2:571825‑M2/86‑149 (MQ=255)                                                               acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                >  1:59845‑M2/85‑149 (MQ=255)                                                               acgcacatacaacggaggggggctgcgattttcaaccagccgagggggagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                <  1:459392‑M2/65‑1 (MQ=255)                                                                   acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTAAAAAAGTTTTTCGATACGCTAgcaagcaa                                                                            >  2:732564‑M2/81‑149 (MQ=255)                                                                       acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgGATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                        <  2:641524‑M2/73‑1 (MQ=255)                                                                        caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                       >  1:655613‑M2/76‑149 (MQ=255)                                                                          accgaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCCTACTTTCGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                     <  1:202460‑M2/76‑1 (MQ=255)                                                                             gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaccgatcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                  >  1:994596‑M2/71‑149 (MQ=255)                                                                              agggtggctgcgattttcaataatgcgtgattcagatctcacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                 <  1:436377‑M2/80‑1 (MQ=255)                                                                                gggggctgcgattttaaataattcgtgatgcagataacaaaaaacactcaattagttaagatgaatgTATAAGTCATACTTTTGTTATGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGAAAAAAAACAAAAATGAAACAGGa                                                               <  2:401936‑M2/82‑1 (MQ=255)                                                                                   ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTAGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                            <  1:913798‑M2/85‑1 (MQ=255)                                                                                     ctgtgattttcaataaggcgtgatgcagattacacaaaacattcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTTTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                          <  2:282658‑M2/87‑1 (MQ=255)                                                                                            tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAACCAAAAATGAAAAAGGAATCATCGAAATg                                                   >  2:158724‑M2/56‑149 (MQ=255)                                                                                              tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCGATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATAGAAATggg                                                 >  1:933957‑M2/54‑149 (MQ=255)                                                                                                 ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTCCGCTAGCAAGCAANAATGAAACAGGAATAANCGCAATGGGATg                                              >  1:172893‑M2/51‑149 (MQ=255)                                                                                                      gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTCCGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcg                                         >  2:191340‑M2/46‑149 (MQ=255)                                                                                                        gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                       >  1:328441‑M2/44‑149 (MQ=255)                                                                                                          gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATCTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                     <  1:970414‑M2/108‑1 (MQ=255)                                                                                                          gataaagatcacacaaaacactccattacttaaaataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAACGTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                     <  2:172893‑M2/108‑1 (MQ=255)                                                                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                   >  1:785032‑M2/40‑149 (MQ=255)                                                                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGATTAGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAGGGGATGTTGCGCACAGt                                   >  1:199919‑M2/40‑149 (MQ=255)                                                                                                             gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTc                                  <  2:328441‑M2/111‑1 (MQ=255)                                                                                                                 atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                              >  2:340495‑M2/35‑149 (MQ=255)                                                                                                                 atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTAaaaa                              >  2:1029689‑M2/35‑149 (MQ=255)                                                                                                                 atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAACAAAAATGAAACAGGAATAATCGAATTGGGATGTTGCGCACAGTCaaaa                              >  1:140100‑M2/35‑149 (MQ=255)                                                                                                                      acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACt                         <  1:340495‑M2/120‑1 (MQ=255)                                                                                                                        agaacactcaattacttaacataaaggTATAAGTCATACTTTTGTTTTGCGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAACAGCATAACAGAAATGGGATGTTGCGCAAAGTAAAGATAAATCa                       >  2:745627‑M2/28‑149 (MQ=255)                                                                                                                        aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                       >  1:249300‑M2/28‑149 (MQ=255)                                                                                                                           acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCg                    >  1:866805‑M2/25‑149 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGTCATAATCAAAATGGTATGTTGTGCACAGTCAAAATAACTCACCGAAatta              >  2:238987‑M2/19‑149 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata              >  1:233517‑M2/19‑149 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata              >  1:727663‑M2/19‑149 (MQ=255)                                                                                                                                    taattaacataaatgTATAAGGCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc           <  2:394828‑M2/134‑1 (MQ=255)                                                                                                                                       ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGGCAAAATAACTCACCGTAAATAatcatc        >  1:161313‑M2/13‑149 (MQ=255)                                                                                                                                       ttaacataaatgGATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc        <  2:198275‑M2/137‑1 (MQ=255)                                                                                                                                        taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCt       <  1:327732‑M2/138‑1 (MQ=255)                                                                                                                                            ataaatgTATAAGTCATACTTTTGTTCTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCta   <  2:77497‑M2/142‑1 (MQ=255)                                                                                                                                             taaatgTATAAGTCATACTTTTGTTTTGGGTGTACTACAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAATAACCGTAAATAATCATCTGCtat  >  2:490416‑M2/7‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913_3_pae_tpiA/1979133‑1979421

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: ‑