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 20563 5175 66 [0] [0] 65 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)65 (0.800) 39/96 0.2 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

ACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTG  >  NC_000913/20515‑20610
                                                 |                                              
agggtgcgttgaagatatgcgagcacctgtaaagtggcgggggtcactct                                                >  1:4026343‑M2/50‑50 (MQ=255)
 gggtgcgttgaagatatgcgagcacctgtaaagtggcggggatcactct                                                <  1:4736483‑M2/1‑1 (MQ=255)
 gggtgcgttgaagatatgcgagcacctgtaaagtggcggggatcactct                                                <  1:1706497‑M2/1‑1 (MQ=255)
 gggtgcgttgaagatatgcgagcacctgtaaagtggcggggatcactct                                                <  1:1547633‑M2/1‑1 (MQ=255)
 gggtgcgttgaagatatgcgagcacctgtaaagtggcggggatcactcTa                                               <  1:1604615‑M2/2‑1 (MQ=255)
   gtgcgttgaagatatgcgagcacctgtaaagtggcgggggtccatcTAcc                                             >  1:3144092‑M2/47‑50 (MQ=255)
    tgcgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCt                                            <  1:1577342‑M2/5‑1 (MQ=255)
      cgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTc                                           <  1:5057675‑M2/6‑1 (MQ=255)
      cgttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCa                                          <  1:2683323‑M2/7‑1 (MQ=255)
        ttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAt                                        <  1:750319‑M2/9‑1 (MQ=255)
        ttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAt                                        <  1:5423329‑M2/9‑1 (MQ=255)
          gaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAtg                                       <  1:1482622‑M2/10‑1 (MQ=255)
          gaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAtg                                       <  1:140447‑M2/10‑1 (MQ=255)
           aagatatgcgagcacctgtaaagtggcgggggtcactcTACCCCAAtgtg                                     >  1:1033101‑M2/39‑50 (MQ=255)
           aagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAtgt                                      <  1:2874831‑M2/11‑1 (MQ=255)
            agatatgcgagcacctgtaaagtggcggggatcactcTACCTCAATgtgt                                    <  1:5782164‑M2/13‑1 (MQ=255)
            agatatgcgagcacctgtaaagtggcggggatcactcTACCTCAATgtgt                                    <  1:769481‑M2/13‑1 (MQ=255)
              atatgcgagcacctgtaaagtggcgggggtcactcCACCTCAATGTGTAt                                  >  1:15130‑M2/36‑50 (MQ=255)
              atatgcgagcacctgtaaagtggcgggggtcactcCACCCCAATGTGTAt                                  >  1:6237554‑M2/36‑50 (MQ=255)
              atatgcgagcacctgtaaagtggcgggggtcacccTACCTCAATTTGTa                                   >  1:648951‑M2/36‑49 (MQ=255)
              atatgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTAt                                  <  1:5460272‑M2/15‑1 (MQ=255)
              atatgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTAt                                  <  1:3252139‑M2/15‑1 (MQ=255)
                atgcgagcacctgtaaagtggcgggggtcactcCACCTCCATTTGtatta                                >  1:6096996‑M2/34‑48 (MQ=255)
                 tgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATcac                               <  1:4699724‑M2/18‑1 (MQ=255)
                   cgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACaa                             <  1:5871566‑M2/20‑1 (MQ=255)
                     agcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAata                           <  1:567636‑M2/22‑1 (MQ=255)
                     agcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAata                           <  1:5294300‑M2/22‑1 (MQ=255)
                      gcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAatat                          <  1:5744586‑M2/23‑1 (MQ=255)
                       cacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAATATc                         <  1:1611951‑M2/24‑1 (MQ=255)
                         cctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCa                       <  1:3316283‑M2/26‑1 (MQ=255)
                         cctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCa                       <  1:5171374‑M2/26‑1 (MQ=255)
                             taaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCata                     <  1:3442023‑M2/28‑1 (MQ=255)
                             taaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCATAtt                   <  1:3857415‑M2/30‑1 (MQ=255)
                               aagtggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCt                 <  1:3701119‑M2/32‑1 (MQ=255)
                               aagtggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCt                 <  1:3119886‑M2/32‑1 (MQ=255)
                                 gtggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCttt               <  1:582004‑M2/34‑1 (MQ=255)
                                  tggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtg              <  1:54597‑M2/35‑1 (MQ=255)
                                  tggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtg              >  1:3909393‑M2/16‑50 (MQ=255)
                                    gcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtgtg            <  1:2052661‑M2/37‑1 (MQ=255)
                                    gcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtgtg            <  1:5115108‑M2/37‑1 (MQ=255)
                                    gcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtgt             <  1:3703107‑M2/36‑1 (MQ=255)
                                    gcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtgt             <  1:4485380‑M2/36‑1 (MQ=255)
                                      ggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTggg          >  1:4249238‑M2/12‑50 (MQ=255)
                                      ggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTggg          >  1:417246‑M2/12‑50 (MQ=255)
                                       gggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg         >  1:546329‑M2/11‑50 (MQ=255)
                                        ggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg         >  1:596003‑M2/10‑49 (MQ=255)
                                        ggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg         >  1:1335083‑M2/10‑49 (MQ=255)
                                         gatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGa       <  1:49212‑M2/42‑1 (MQ=255)
                                         gatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGa       <  1:4902764‑M2/42‑1 (MQ=255)
                                         gatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGa       <  1:5277081‑M2/42‑1 (MQ=255)
                                          atcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAg      >  1:5463933‑M2/8‑50 (MQ=255)
                                          atcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAg      >  1:2126813‑M2/8‑50 (MQ=255)
                                           tcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGt     >  1:4703782‑M2/7‑50 (MQ=255)
                                           tcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGt     <  1:2131046‑M2/44‑1 (MQ=255)
                                           tcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGt     >  1:1655083‑M2/7‑50 (MQ=255)
                                            cactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTc    >  1:1745730‑M2/6‑50 (MQ=255)
                                            cactcTACCTCAATGAGTATCACAATATCCATATTCTTTGTGGGGGAGTc    >  1:855402‑M2/6‑50 (MQ=255)
                                             actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTc    <  1:4921605‑M2/45‑1 (MQ=255)
                                             actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTc    <  1:1340366‑M2/45‑1 (MQ=255)
                                              ctcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  <  1:4046351‑M2/47‑1 (MQ=255)
                                              ctcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  >  1:3452974‑M2/4‑50 (MQ=255)
                                              ctcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  <  1:2199418‑M2/47‑1 (MQ=255)
                                               tcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  >  1:407980‑M2/3‑49 (MQ=255)
                                               tcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  >  1:2091577‑M2/3‑49 (MQ=255)
                                               tcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg  >  1:1432281‑M2/3‑49 (MQ=255)
                                                 |                                              
ACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTG  >  NC_000913/20515‑20610

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