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	51 [0]	[0] 52	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)	49 (0.870)	40/274	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

ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_pae_tpiA/1979135‑1979427                                                                                                                                                 |                                                                                                                                                     atgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                                  <  1:507405‑M2/5‑1 (MQ=255)             acttatctaagatttttcctaaatcgacgcaactgtactcgtcactccacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                                                                      >  1:440957‑M2/133‑149 (MQ=255)                 atctaagatttttcctaaatcgacggaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                                  <  2:182979‑M2/21‑1 (MQ=255)                   ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                >  1:751479‑M2/127‑149 (MQ=255)                   ctaaaattttacttaattcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  2:161443‑M2/23‑1 (MQ=255)                               ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggagtggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaaaataaatgTATAAGGCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    >  2:363421‑M2/115‑149 (MQ=255)                               ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    <  2:985319‑M2/35‑1 (MQ=255)                               ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaagggggctgcgattttcaataatgcttgatgcagatcacgcaaaacactcaattacttaacataaatgTATAACGCATATTTTCACTTTGGGTGTATTTCAAt                                                                                                                    >  2:205142‑M2/115‑149 (MQ=255)                               ctaaatcgacgcaactgtactcgtcactacaagcacatacaacggagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    <  1:734478‑M2/35‑1 (MQ=255)                                   atcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                                <  2:342595‑M2/39‑1 (MQ=255)                                         gcatttgtactcgtcactcaacgcacatagaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                          <  2:803099‑M2/45‑1 (MQ=255)                                         ccaactgtactcaccactacacgcacatacaacagaggggggctgcgactatccacaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                          <  2:896240‑M2/45‑1 (MQ=255)                                              tgtactcgtcactacacgcacatacaacggaggggggctgcgatgttcaataatacgtgatgcagatcacacaaaacactcaattacttaacaaaaatgAAAAAATCATACTATTGTTATAGGTGTAATTCAAACTGTTAAAAAGtttt                                                                                                     >  1:43199‑M2/100‑149 (MQ=255)                                               gtactcgtaactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                    <  1:753121‑M2/51‑1 (MQ=255)                                                 actcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                                                                  >  2:778880‑M2/97‑149 (MQ=255)                                                   tcgtcactacgcgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTATAAAGTTTttcgac                                                                                                >  1:1017083‑M2/95‑147 (MQ=255)                                                   tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacccaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTCTTGGCTGTATTTCAATCTGATAAAAAATTATTcgct                                                                                                >  1:231667‑M2/95‑149 (MQ=255)                                                   tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                                >  1:494788‑M2/95‑149 (MQ=255)                                                   tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGAGTATTTCAATCTGTTAAAAAGTTTttcgtt                                                                                                >  1:704019‑M2/95‑147 (MQ=255)                                                       tactacacgcacaatcaacggaggggggcttcgattttaaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                            <  2:769704‑M2/59‑1 (MQ=255)                                                              cgcacatccaacggaggtgtgttgcgattttcaataatgcgttatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                                     <  1:352560‑M2/66‑1 (MQ=255)                                                                cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGctagcaagaa                                                                                   >  1:822733‑M2/82‑147 (MQ=255)                                                                         cggaggggagctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAATAACGTTTGCGCTACGCTAGCAAGCAAAAATGaaa                                                                          <  2:224695‑M2/77‑1 (MQ=255)                                                                           gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                        <  2:974968‑M2/79‑1 (MQ=255)                                                                                    tgttattttcaataatgcgtgattcagatctcataaaacattcaaatacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                               <  1:233153‑M2/88‑1 (MQ=255)                                                                                    tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                               >  2:515422‑M2/62‑149 (MQ=255)                                                                                             caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgAATACGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTACGCTAGCAAGCCAAAATGAAACAGGAATAATCGAAATGGGa                                                      >  2:118351‑M2/53‑149 (MQ=255)                                                                                              aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                                     >  2:277175‑M2/52‑149 (MQ=255)                                                                                                    gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGatgttctg                                               >  2:421995‑M2/46‑146 (MQ=255)                                                                                                        gatgcagatcacacaaaacactaaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCTCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  2:122638‑M2/108‑1 (MQ=255)                                                                                                            cagatcacacaaaacactcaattacttaacataaatgTATAAATCATACTTGTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                       >  2:526977‑M2/38‑149 (MQ=255)                                                                                                             agatcacaaaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaa                                      <  1:351623‑M2/113‑1 (MQ=255)                                                                                                              gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATTGAAATGGGATGTTGCGCACAGTCaaa                                     >  2:48818‑M2/36‑149 (MQ=255)                                                                                                                tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                                   >  1:593038‑M2/34‑149 (MQ=255)                                                                                                                 cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATa                                  >  2:1035125‑M2/33‑149 (MQ=255)                                                                                                                  acaaaaaacactctattactaaacataaatgTATAAGTCATACTTTTCTTTTGGGTGTACTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                                 <  2:199251‑M2/118‑1 (MQ=255)                                                                                                                    acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTGAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAaca                               >  1:198972‑M2/30‑148 (MQ=255)                                                                                                                      aaaacactcaattacttaacaaaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                             <  2:649680‑M2/122‑1 (MQ=255)                                                                                                                       aaacactcaatgacttagcataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                            <  1:161036‑M2/123‑1 (MQ=255)                                                                                                                        aatcatcaattatttaacatatactTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCGGTTAAAAAGTGTTTCGCTGCGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAcc                           <  1:983198‑M2/124‑1 (MQ=255)                                                                                                                        aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACtcaac                           >  1:1004921‑M2/26‑147 (MQ=255)                                                                                                                        aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAcc                           >  2:561402‑M2/26‑149 (MQ=255)                                                                                                                        aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAcc                           >  2:766408‑M2/26‑149 (MQ=255)                                                                                                                            ctcaattacttaacattaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa                       <  1:531192‑M2/128‑1 (MQ=255)                                                                                                                             tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa                      >  2:1019691‑M2/21‑149 (MQ=255)                                                                                                                               aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCTAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTCaata                    >  2:52322‑M2/19‑149 (MQ=255)                                                                                                                               aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGCCAAAACAACTCACCGTAaata                    >  1:643128‑M2/19‑149 (MQ=255)                                                                                                                                attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataa                   >  1:252574‑M2/18‑149 (MQ=255)                                                                                                                                        acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATACTCATCTGc           >  1:380447‑M2/10‑149 (MQ=255)                                                                                                                                              atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa     <  2:252574‑M2/146‑1 (MQ=255)                                                                                                                                              atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa     >  1:226586‑M2/4‑149 (MQ=255)                                                                                                                                                 tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACTCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  <  1:363421/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_pae_tpiA/1979135‑1979427

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