Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,110,680 25702 bp→GA [yghD]pitB 20 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3110680 3136381 25702 73 [0] [0] 70 [yghD]–pitB [yghD],yghE,yghF,yghG,pppA,yghJ,glcA,glcB,glcG,glcF,glcE,glcD,glcC,yghO,insH‑9,yghQ,yghR,yghS,yghT,pitB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 31106790 (0.000)67 (0.960)
+GA
47/94 0.0 100% coding (448/537 nt) yghD putative type II secretion system M‑type protein
?NC_000913 3136382 = 0 (0.000)intergenic (‑11/+281) pitB/gss metal phosphate:H(+) symporter PitB/fused glutathionylspermidine amidase/glutathionylspermidine synthetase

CCAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110630‑3110728
                                                 |                                                 
ccAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCata                                                   <  1:3933379/50‑1 (MQ=255)
ccAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCata                                                   <  1:3568508/50‑1 (MQ=255)
 cAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                  <  1:948732‑M1/50‑2 (MQ=255)
 cAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                  <  1:3907109‑M1/50‑2 (MQ=255)
  aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAga                                                 <  1:148173‑M1/50‑3 (MQ=255)
  aGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAg                                                  >  1:114531‑M1/1‑48 (MQ=255)
   ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                                <  1:1079727‑M1/50‑4 (MQ=255)
   ggCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaa                                                >  1:5214048‑M1/1‑47 (MQ=255)
    gCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaac                                               <  1:1256885‑M1/50‑5 (MQ=255)
      ttCTCAGCAACACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             >  1:3571503‑M1/1‑44 (MQ=255)
       tctcAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             >  1:1495326‑M1/1‑43 (MQ=255)
       tctcAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgc                                             >  1:701592‑M1/1‑43 (MQ=255)
           agcagcACTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattt                                        <  1:4024431‑M1/50‑12 (MQ=255)
             cagcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttg                                      <  1:2604080‑M1/50‑14 (MQ=255)
              agcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttg                                      >  1:2992601‑M1/1‑36 (MQ=255)
              agcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttg                                      >  1:3423153‑M1/1‑36 (MQ=255)
               gcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaa                                    >  1:1885183‑M1/1‑35 (MQ=255)
               gcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaa                                    >  1:730414‑M1/1‑35 (MQ=255)
               gcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaa                                    <  1:2926216‑M1/50‑16 (MQ=255)
               gcaCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaa                                    >  1:2452031‑M1/1‑35 (MQ=255)
                caCTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaat                                   >  1:3751208‑M1/1‑34 (MQ=255)
                  cTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaattt                                 <  1:4383385‑M1/50‑19 (MQ=255)
                  cTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatt                                  >  1:3191874‑M1/1‑32 (MQ=255)
                  cTGACATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatt                                  >  1:3064552‑M1/1‑32 (MQ=255)
                     aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                              <  1:540241‑M1/50‑22 (MQ=255)
                     aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttac                               <  1:5024461‑M1/49‑21 (MQ=255)
                     aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttac                               >  1:3990504‑M1/1‑29 (MQ=255)
                     aCATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaattta                                <  1:2466864‑M1/48‑20 (MQ=255)
                      cATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                              >  1:2678506‑M1/1‑28 (MQ=255)
                      cATCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacg                              >  1:2658558‑M1/1‑28 (MQ=255)
                       aTCAATTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaa                            <  1:2312868‑M1/50‑24 (MQ=255)
                          aaTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacg                          >  1:4545385‑M1/1‑24 (MQ=255)
                          aaTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacg                          >  1:566319‑M1/1‑24 (MQ=255)
                           aTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacgga                        <  1:2998150‑M1/50‑28 (MQ=255)
                           aTTTGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacgga                        >  1:5013220‑M1/1‑23 (MQ=255)
                            tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       >  1:3151031‑M1/1‑22 (MQ=255)
                            tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       <  1:1953116‑M1/50‑29 (MQ=255)
                            tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       <  1:378618‑M1/50‑29 (MQ=255)
                            tttGTGTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggac                       <  1:1806314‑M1/50‑29 (MQ=255)
                              tgtgTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacg                      >  1:693006‑M1/1‑20 (MQ=255)
                              tgtgTCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacg                      >  1:695988‑M1/1‑20 (MQ=255)
                               gtgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgca                    >  1:1780405‑M1/1‑19 (MQ=255)
                               gtgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgca                    >  1:4554299‑M1/1‑19 (MQ=255)
                               gtgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgca                    >  1:2166625‑M1/1‑19 (MQ=255)
                                tgtCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcat                   >  1:5055269‑M1/1‑18 (MQ=255)
                                  tCACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcatt                  <  1:4194806‑M1/49‑34 (MQ=255)
                                   cACCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                >  1:743289‑M1/1‑15 (MQ=255)
                                    aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                <  1:4426230‑M1/49‑36 (MQ=255)
                                    aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                <  1:5513341‑M1/49‑36 (MQ=255)
                                    aCCCGCAGCGCATAgaacgccattttgaatttacgaacggacgcattat                >  1:1941216‑M1/1‑14 (MQ=255)
                                     cccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatca              <  1:3561058‑M1/50‑38 (MQ=255)
                                      ccGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcag             <  1:942953‑M1/50‑39 (MQ=255)
                                       cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt            <  1:477892‑M1/50‑40 (MQ=255)
                                       cGCAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagt            >  1:3109838‑M1/1‑11 (MQ=255)
                                         cAGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtga          <  1:4823540‑M1/50‑42 (MQ=255)
                                          aGCGCATAgaacgccattttgaatttacgaacggacgcattatcagtgac         >  1:3177520‑M1/1‑8 (MQ=255)
                                           gcgcATAgaacgccattttgaatttacgaacggacgcattatcagtgact        <  1:4878667‑M1/50‑44 (MQ=255)
                                            cgcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:3952468‑M1/1‑6 (MQ=255)
                                            cgcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:2537736‑M1/1‑6 (MQ=255)
                                            cgcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:5226974‑M1/1‑6 (MQ=255)
                                            cgcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:4906506‑M1/1‑6 (MQ=255)
                                             gcATAgaacgccattttgaatttacgaacggacgcattatcagtgacttt      <  1:3356822‑M1/50‑46 (MQ=255)
                                             gcATAgaacgccattttgaatttacgaacggacgcattatcagtgactt       >  1:1318424‑M1/1‑5 (MQ=255)
                                              cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttta     <  1:4216709‑M1/50‑47 (MQ=255)
                                              cATAgaacgccattttgaatttacgaacggacgcattatcagtgacttta     <  1:2277648‑M1/50‑47 (MQ=255)
                                               aTAgaacgccattttgaatttacgaacggacgcattatcagtgactttaa    >  1:3482907‑M1/1‑3 (MQ=255)
                                               aTAgaacgccattttgaatttacgaacggacgcattatcagtgactttaa    <  1:4406949‑M1/50‑48 (MQ=255)
                                               aTAgaacgccattttgaatttacgaacggacgcattatcagtgactttaa    <  1:2280738‑M1/50‑48 (MQ=255)
                                                tAgaacgccattttgaatttacgaacggacgcattatcagtgactttaac   <  1:514654‑M1/50‑49 (MQ=255)
                                                tAgaacgccattttgaatttacgaacggacgcattatcagtgactttaac   >  1:5411700‑M1/1‑2 (MQ=255)
                                                tAgaacgccattttgaatttacgaacggacgcattatcagtgactttaac   <  1:5024352‑M1/50‑49 (MQ=255)
                                                 agaacgccattttgaatttacgaacggacgcattatcagtgactttaacg  >  1:271762‑M1/1‑1 (MQ=255)
                                                 agaacgccattttgaatttacgaacggacgcattatcagtgactttaacg  >  1:1525793‑M1/1‑1 (MQ=255)
                                                 |                                                 
CCAGGCTTCTCAGCAGCACTGACATCAATTTGTGTCACCCGCAGCGCATATTTTTCATCCAGTGCGTTTAACCATTTCAGCAGGTCATTAAACACCACA  >  NC_000913/3110630‑3110728

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