Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 15,389 Δ5,175 bp insL‑1insA‑1 insL‑1, hokC, mokC, sokC, nhaA, nhaR, insB‑1, insA‑1

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 15389–15438 20563 5126–5175 72 [0] [0] 74 insL‑1–insA‑1 insL‑1,hokC,mokC,sokC,nhaA,nhaR,insB‑1,insA‑1

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 153880 (0.000)72 (0.780) 46/96 0.1 100% noncoding (2/1345 nt) IS186 repeat region
?NC_000913 20564 = 0 (0.000)intergenic (‑56/+251) insA‑1/rpsT IS1 protein InsA/30S ribosomal subunit protein S20

CCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTG  >  NC_000913/20516‑20610
                                                |                                              
gggtgcgttgaagatatgcgagcacctgtaaagtggcggggatcactcTa                                               <  1:6191113‑M2/2‑1 (MQ=255)
  gtgcgttgaagatatgcgagcacctgtaaagtggcggggatcactcTAcc                                             <  1:6923673‑M2/4‑1 (MQ=255)
   tgcgttgaagatatgcgagcacctgtaaagtggcggggatcactcTaccc                                            >  1:2271572‑M2/46‑49 (MQ=255)
   tgcgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCt                                            <  1:6837201‑M2/5‑1 (MQ=255)
   tgcgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCt                                            <  1:3511114‑M2/5‑1 (MQ=255)
     cgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTc                                           <  1:4610687‑M2/6‑1 (MQ=255)
     cgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCa                                          <  1:1551513‑M2/7‑1 (MQ=255)
     cgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCa                                          <  1:6979935‑M2/7‑1 (MQ=255)
     cgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCa                                          <  1:6839638‑M2/7‑1 (MQ=255)
     cgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCa                                          <  1:5096426‑M2/7‑1 (MQ=255)
      gttgaagatatgcgagcacctgtaaagtggcgggggtcactcTACCtcc                                          >  1:7065756‑M2/43‑48 (MQ=255)
      gttgaagatatgcgagcacctgtaaagtggcgggggtcactcTACCTCaa                                         >  1:6110144‑M2/43‑50 (MQ=255)
      gttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCaa                                         <  1:98815‑M2/8‑1 (MQ=255)
       ttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAt                                        <  1:5742611‑M2/9‑1 (MQ=255)
         gaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAtg                                       <  1:3433968‑M2/10‑1 (MQ=255)
         gaagatatgcgagcacctgtaaagtggcggggaatactcCACCTCAAtgt                                      >  1:4541484‑M2/40‑50 (MQ=255)
          aagatatgcgagcacctgtaaagtggcgggggtcactcTACCCCCAtgtg                                     >  1:7008809‑M2/39‑50 (MQ=255)
           agatatgcgagcacctgtaaagtggcgggggtcactcTCCCTCCATttgt                                    >  1:4086081‑M2/38‑50 (MQ=255)
           agatatgcgagcacctgtaaagtggcggggatcactcTACCTCAATgtgt                                    <  1:5800320‑M2/13‑1 (MQ=255)
             atatgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTAt                                  <  1:2023085‑M2/15‑1 (MQ=255)
             atatgcgagcacctgtaaagtggcggcgatcactcTACCTCAATGTGTAt                                  >  1:5033579‑M2/36‑50 (MQ=255)
              tatgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATc                                 <  1:6326413‑M2/16‑1 (MQ=255)
              tatgcgagcacctgtaaagtggcggggaccactcTACCTCAATGTGTAt                                  >  1:4642619‑M2/35‑49 (MQ=255)
                tgcgagcacctgtaaagtggcgggggtcactcCTCCTCAATGTGTATcac                               >  1:4460321‑M2/33‑50 (MQ=255)
                tgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATcac                               <  1:6435133‑M2/18‑1 (MQ=255)
                tgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATcac                               <  1:4130457‑M2/18‑1 (MQ=255)
                tgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATcac                               <  1:973319‑M2/18‑1 (MQ=255)
                  cgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATcaca                              <  1:4007357‑M2/19‑1 (MQ=255)
                  cgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACaa                             <  1:4070211‑M2/20‑1 (MQ=255)
                   gagcacctgtaaagtggcgggggtcactcCACCTCCATTTGTATCACaa                             >  1:6544500‑M2/30‑49 (MQ=255)
                   gagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAat                            <  1:3339928‑M2/21‑1 (MQ=255)
                   gagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAat                            <  1:3209470‑M2/21‑1 (MQ=255)
                       acctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAATATcc                        <  1:5491246‑M2/25‑1 (MQ=255)
                        cctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCa                       <  1:2367770‑M2/26‑1 (MQ=255)
                            taaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCatat                    <  1:3696578‑M2/29‑1 (MQ=255)
                            taaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCata                     <  1:3281573‑M2/28‑1 (MQ=255)
                            taaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCATAtt                   <  1:6831960‑M2/30‑1 (MQ=255)
                               agtggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCtt                <  1:3217143‑M2/33‑1 (MQ=255)
                                gtggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCttt               <  1:3749003‑M2/34‑1 (MQ=255)
                                 tggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtg              >  1:2695670‑M2/16‑50 (MQ=255)
                                 tggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtg              >  1:2714998‑M2/16‑50 (MQ=255)
                                 tggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtg              >  1:5776918‑M2/16‑50 (MQ=255)
                                   gcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtgt             <  1:5567329‑M2/36‑1 (MQ=255)
                                   gcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtgt             <  1:4624788‑M2/36‑1 (MQ=255)
                                    cggggatcactgTACCTCAATGTGTATCACAATATCCATATTCTTTGTgg           <  1:6737644‑M2/38‑1 (MQ=255)
                                    cggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgg           >  1:1857155‑M2/13‑50 (MQ=255)
                                    cggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgg           >  1:1518707‑M2/13‑50 (MQ=255)
                                    cggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgg           >  1:1508141‑M2/13‑50 (MQ=255)
                                    cggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgg           >  1:5462145‑M2/13‑50 (MQ=255)
                                     ggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTggg          >  1:3311392‑M2/12‑50 (MQ=255)
                                     ggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTggg          >  1:342803‑M2/12‑50 (MQ=255)
                                     ggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTggg          <  1:1595020‑M2/39‑1 (MQ=255)
                                      gggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg         >  1:2132145‑M2/11‑50 (MQ=255)
                                      gggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg         <  1:2856533‑M2/40‑1 (MQ=255)
                                       ggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg         >  1:1903075‑M2/10‑49 (MQ=255)
                                        gatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg         >  1:3065676‑M2/9‑48 (MQ=255)
                                         atcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAg      >  1:6719639‑M2/8‑50 (MQ=255)
                                          tcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGt     >  1:1908930‑M2/7‑50 (MQ=255)
                                          tcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGt     >  1:4034527‑M2/7‑50 (MQ=255)
                                          tcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGt     >  1:4079724‑M2/7‑50 (MQ=255)
                                          tcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGt     >  1:4581036‑M2/7‑50 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTc    <  1:2150968‑M2/45‑1 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTc    <  1:5615820‑M2/45‑1 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt   >  1:1889288‑M2/5‑50 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt   <  1:6910858‑M2/46‑1 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt   >  1:2874344‑M2/5‑50 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt   <  1:4859190‑M2/46‑1 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt   >  1:5936182‑M2/5‑50 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt   >  1:5485884‑M2/5‑50 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt   >  1:5695664‑M2/5‑50 (MQ=255)
                                              tcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  >  1:1284279‑M2/3‑49 (MQ=255)
                                              tcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  >  1:55408‑M2/3‑49 (MQ=255)
                                                tacctCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  >  1:3787756/1‑47 (MQ=255)
                                                tacctCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  >  1:621973/1‑47 (MQ=255)
                                                |                                              
CCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTG  >  NC_000913/20516‑20610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: