Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913_3_bme_pgi 1,978,503 Δ776 bp 100% insB1insA insB1, insA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_bme_pgi 1978503–1979268 1979278 11–776 62 [0] [0] 62 insB1–insA insB1,insA

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

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_bme_pgi/1979133‑1979427
                                                                                                                                                  |                                                                                                                                                    
aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                    <  1:631660‑M2/3‑1 (MQ=255)
  atgggtctgtggacttatcgaagttgtttcctatatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacataacactcaattacttaacataaatgTATaa                                                                                                                                                  <  2:93063‑M2/5‑1 (MQ=255)
         tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcggacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  1:744489‑M2/138‑149 (MQ=255)
         tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  2:866935‑M2/138‑149 (MQ=255)
         tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  1:450347‑M2/138‑149 (MQ=255)
         tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggagggctgcgattttcaatactcccccatgcacatcacacaaaacacacaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  2:317639‑M2/138‑149 (MQ=255)
                    ctaatatttttcttaaatcggcgcaactgtactcgtcactacacgcacatacaaagaaggatggatgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATActtttgttttgt                                                                                                                                <  2:680457‑M2/23‑2 (MQ=255)
                    ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacagaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                <  2:199684‑M2/23‑1 (MQ=255)
                     taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggcggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactaaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                               >  1:612784‑M2/126‑149 (MQ=255)
                     taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATACTTTTGTTTTggg                                                                                                                               >  1:463454‑M2/126‑149 (MQ=255)
                     taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgagtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTtgg                                                                                                                               >  2:232872‑M2/126‑149 (MQ=255)
                         atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgagtgatgcaaatcacacaaaaaactcaattacttaacataaatgAATAAGCCATAATTTTGTTTCGGGtata                                                                                                                           >  2:655016‑M2/122‑149 (MQ=255)
                         atttttcctaaatcgacgcaactgcactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                                                                           >  1:779255‑M2/122‑149 (MQ=255)
                         atttttcctaaatcaacgcaactgcactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                                                                           >  1:779004‑M2/122‑149 (MQ=255)
                                 taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                   <  1:441556‑M2/36‑1 (MQ=255)
                                 taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGATTTGGGTGTATTTCAATc                                                                                                                   >  1:632816‑M2/114‑149 (MQ=255)
                                     tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                               >  2:754777‑M2/110‑149 (MQ=255)
                                        acgcaactgtactcgtcactacactcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                            <  1:674125‑M2/43‑1 (MQ=255)
                                          tcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTTACACTTTCTTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                          <  1:892238‑M2/45‑1 (MQ=255)
                                          ggaactgtactcgacactaaacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                          <  2:626734‑M2/45‑1 (MQ=255)
                                            aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcagacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTATTAGATGTATTTCAATATGATAAAAAGt                                                                                                        >  2:284672‑M2/103‑149 (MQ=255)
                                            aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                                                                        >  1:756477‑M2/103‑149 (MQ=255)
                                            aactgtactcgtcactacacgcacatacaacggaggagggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTATATTTCAAtctgttaccaacg                                                                                                        >  2:444072‑M2/103‑143 (MQ=255)
                                                gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                    <  1:819274‑M2/51‑1 (MQ=255)
                                                gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                    <  1:218656‑M2/51‑1 (MQ=255)
                                                   ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                                                                 <  2:868084‑M2/54‑1 (MQ=255)
                                                    tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTCTGGCTGTATTACAATATATTAAAAAGCTtttctca                                                                                                >  2:806669‑M2/95‑146 (MQ=255)
                                                    tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagataacacaaaagcctcaattacttaacataaatgTATAAATCATACGTTTATTTTGGGTGTATGTCAATCTGTTAAAAAATTTTTcgct                                                                                                >  1:525097‑M2/95‑149 (MQ=255)
                                                       tccctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaactctcaattacttaacattaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                                                                             <  1:270452‑M2/58‑1 (MQ=255)
                                                                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTCCGCTAGCAGGCacaaa                                                                               >  2:510132‑M2/78‑149 (MQ=255)
                                                                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGCTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                               >  2:679316‑M2/78‑149 (MQ=255)
                                                                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                                                                  <  1:621065‑M2/53‑1 (MQ=255)
                                                                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                                                                  >  2:621065‑M2/77‑129 (MQ=255)
                                                                         acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                           >  2:124175‑M2/74‑149 (MQ=255)
                                                                          cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                          >  1:446426‑M2/73‑149 (MQ=255)
                                                                             aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                       <  1:113132‑M2/80‑1 (MQ=255)
                                                                                       cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGATTTTCGCTACTCTAGCCAGCCAAAATGGAAAAGCAATAatcca                                                             >  2:226552‑M2/60‑147 (MQ=255)
                                                                                         attttcaataatgcgtgatgcagatcaaaaaaatcattcaattacttaacataaatgTATAAGTTATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                           <  2:721474‑M2/92‑1 (MQ=255)
                                                                                           tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATATGTTAAAAAGTTTTTAGCTACGATAGCAAGCAAAAATGAAACAGGAATAATAGAAATg                                                         >  1:742661‑M2/56‑149 (MQ=255)
                                                                                              caataatgcatgatgaagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAACGAAACAGGAATAATCGAAATGGGa                                                      <  2:553761‑M2/97‑1 (MQ=255)
                                                                                                       gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGATGTATTTTAATCTGTTAAAAACTTTTTCGCTACGCTAGCAAGCAAAAATAACACAGGAATAATCGAAATGGTATGttgcgac                                             >  2:337074‑M2/44‑147 (MQ=255)
                                                                                                         gatgcagattatacaaaatactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAACAGTTTCTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  1:849432‑M2/108‑1 (MQ=255)
                                                                                                         gatgcagatcacacaaaatacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAATTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  1:393741‑M2/108‑1 (MQ=255)
                                                                                                         gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  2:549031‑M2/108‑1 (MQ=255)
                                                                                                             cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                       <  1:820522‑M2/112‑1 (MQ=255)
                                                                                                               gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTaaaa                                     >  1:86839‑M2/36‑149 (MQ=255)
                                                                                                                atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAAGTTACGCACAGTCaaaa                                    >  2:890852‑M2/35‑149 (MQ=255)
                                                                                                                 ttacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                                   <  1:32511‑M2/116‑1 (MQ=255)
                                                                                                                     acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACt                               >  2:785851‑M2/30‑149 (MQ=255)
                                                                                                                     aaaaaacactcaattacttaacataaatgTATAAGTTATACTTTTGTTTTAGGTGTATTTCAATCTGTTAAAAAGAATTTCGCGACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACt                               <  1:123227‑M2/120‑1 (MQ=255)
                                                                                                                        aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTCCGATAGGAAGCAAAAATTAAAAAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaacttca                            >  2:802244‑M2/27‑146 (MQ=255)
                                                                                                                              gcaattacataaaataaatgAAATAGTCATAGTTTTGTTTTCGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAGGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa                      <  2:93153‑M2/129‑1 (MQ=255)
                                                                                                                                aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACTCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAATCACCGTAaata                    >  1:514366‑M2/19‑149 (MQ=255)
                                                                                                                                 attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataa                   >  1:483245‑M2/18‑149 (MQ=255)
                                                                                                                                  ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                                <  2:168327‑M2/103‑1 (MQ=255)
                                                                                                                                  ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                                >  1:168327‑M2/17‑119 (MQ=255)
                                                                                                                                    acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAAAAatca                >  1:762385‑M2/15‑149 (MQ=255)
                                                                                                                                     ttgacaataaatgAATAAGTTATACTTTTGTTTTGGTTGTATTTCTATCTTTTAAAAATTTTTTTGCTACGCTAGAAAGCCAAAATGAAACAGGAAAAATCGAAAGGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat               <  1:450200‑M2/136‑1 (MQ=255)
                                                                                                                                     cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTACCCGTAAACAaacat               >  2:668448‑M2/14‑149 (MQ=255)
                                                                                                                                     cttaacataaatgTAAAAGTTACACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCATCGTAAATAatcat               <  2:446426‑M2/136‑1 (MQ=255)
                                                                                                                                           ataaatgTATAAGTCATACTTTTTTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGAAAAAATGAAAAAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAATTACAGTAAATAATCATCTGctc         >  2:534165‑M2/8‑148 (MQ=255)
                                                                                                                                                  tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATATGCTataaataa  >  1:863391/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_bme_pgi/1979133‑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: 
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.