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 50 [0] [0] 49 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)45 (0.900) 30/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

CCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGG  >  NC_000913/20516‑20611
                                                |                                               
gggtgcgttgaagatatgcgagcacctgtaaagtggcggggatcactct                                                 <  1:1953573‑M2/1‑1 (MQ=255)
     cgttgaagatatgcgagcacctgtaaagtggcgggggtcactctacccc                                            >  1:52573‑M2/44‑47 (MQ=255)
      gttgaagatatgcgagcacctgtaaagtggcgggggtcacccCACctcct                                          >  1:1240364‑M2/43‑48 (MQ=255)
      gttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCaa                                          <  1:3900528‑M2/8‑1 (MQ=255)
      gttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCaa                                          <  1:2379865‑M2/8‑1 (MQ=255)
       ttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAt                                         <  1:99072‑M2/9‑1 (MQ=255)
       ttgaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAt                                         <  1:7367‑M2/9‑1 (MQ=255)
        tgaagatatgcgagcacctgtaaagtggcgggggtcactcTACCTCAAtg                                        >  1:1262740‑M2/41‑50 (MQ=255)
         gaagatatgcgagcacctgtaaagtggcgggggttactcTTCCTCCAtgt                                       >  1:1006709‑M2/40‑50 (MQ=255)
         gaagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAtg                                        <  1:2179943‑M2/10‑1 (MQ=255)
          aagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAtgtg                                      <  1:3509947‑M2/12‑1 (MQ=255)
          aagatatgcgagcacctgtaaagtggcggggatcactcTACCTCAAtgtg                                      <  1:2707079‑M2/12‑1 (MQ=255)
           agatatgcgagcacctgtaaagtggcgggggtcactcCACCTCCAtgtg                                      >  1:3640134‑M2/38‑49 (MQ=255)
           agatatgcgagcacctgtaaagtggcggggatcactcTACCTCAATgtgt                                     <  1:745226‑M2/13‑1 (MQ=255)
           agatatgcgagcacctgtaaagtggcggggatcactcTACCTCAATgtgt                                     <  1:2739631‑M2/13‑1 (MQ=255)
           agatatgcgagcacctgtaaagtggcggggatcactcTACCTCAATgtgt                                     <  1:1732126‑M2/13‑1 (MQ=255)
           agatatgcgagcacctgtaaagtggcggggatcactcTACCTCAATgtgt                                     <  1:1390436‑M2/13‑1 (MQ=255)
            gatatgcgagcacctgtaaagtggcgggggtcaaccTACCTCAAAGTGTa                                    >  1:1436153‑M2/37‑50 (MQ=255)
             atatgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTAt                                   <  1:897340‑M2/15‑1 (MQ=255)
              tatgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATc                                  <  1:1145161‑M2/16‑1 (MQ=255)
              tatgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATc                                  <  1:685296‑M2/16‑1 (MQ=255)
                tgcgagcacctgtaaagtggcggggattacacCACCTCAATGTGtatta                                 >  1:139144‑M2/33‑47 (MQ=255)
                tgcgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATca                                 <  1:1941870‑M2/17‑1 (MQ=255)
                  cgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATcaca                               <  1:2205860‑M2/19‑1 (MQ=255)
                  cgagcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACaa                              <  1:3057689‑M2/20‑1 (MQ=255)
                    agcacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAata                            <  1:1429566‑M2/22‑1 (MQ=255)
                      cacctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAATATc                          <  1:1629297‑M2/24‑1 (MQ=255)
                        cctgtaaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCa                        <  1:3243566‑M2/26‑1 (MQ=255)
                            taaagtggcggggatcactcTACCTCAATGTGTATCACAATATCCata                      <  1:1199590‑M2/28‑1 (MQ=255)
                              aagtggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTc                   <  1:1373827‑M2/31‑1 (MQ=255)
                              aagtggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTc                   <  1:2996193‑M2/31‑1 (MQ=255)
                                gtggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCttt                <  1:3145514‑M2/34‑1 (MQ=255)
                                  ggcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtgt              >  1:1101201‑M2/15‑50 (MQ=255)
                                   gcggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTtgtg             >  1:3494708‑M2/14‑50 (MQ=255)
                                     ggggatcactcTACCTCAATGTGTATCACATTATCCATATTCTTTGTggg           >  1:1354537‑M2/12‑50 (MQ=255)
                                     ggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTggg           >  1:1708503‑M2/12‑50 (MQ=255)
                                     ggggatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTggg           >  1:2453659‑M2/12‑50 (MQ=255)
                                        gatcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg          >  1:1112647‑M2/9‑48 (MQ=255)
                                          tcactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGt      <  1:2975548‑M2/44‑1 (MQ=255)
                                           cactcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTc     >  1:208463‑M2/6‑50 (MQ=255)
                                            actcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTc     <  1:745221‑M2/45‑1 (MQ=255)
                                             ctcTACCTCAATGTGTATCACAATATCCATATTATTTGTGGGGGAGTCTg   <  1:3770063‑M2/47‑1 (MQ=255)
                                              tcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg   >  1:3502003‑M2/3‑49 (MQ=255)
                                              tcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg   >  1:35571‑M2/3‑49 (MQ=255)
                                              tcTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg   >  1:798557‑M2/3‑49 (MQ=255)
                                               cTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg   >  1:1034914‑M2/2‑48 (MQ=255)
                                               cTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg   >  1:1301412‑M2/2‑48 (MQ=255)
                                                tacctCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTgg  <  1:2215096/48‑1 (MQ=255)
                                                tacctCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTg   >  1:1043020/1‑47 (MQ=255)
                                                |                                               
CCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGG  >  NC_000913/20516‑20611

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