Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,048,805 Δ48,085 bp [gfcB]ycdU 45 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1048805 1096889 48085 23 [0] [2] 27 [gfcB]–ycdU [gfcB],gfcA,insA4,insB4,cspH,cspG,ymcF,ymcE,gnsA,yccM,torS,torT,torR,torC,torA,torD,cbpM,cbpA,yccE,agp,yccJ,wrbA,ymdF,rutG,rutF,rutE,rutD,rutC,rutB,rutA,rutR,putA,ymdG,putP,efeO,efeB,phoH,pgaD,pgaC,pgaB,pgaA,dgcT,insF4,insE4,ycdU

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 10488040 (0.000)23 (0.500) 21/480 0.9 100% coding (529/645 nt) gfcB lipoprotein GfcB
?NC_000913 1096890 = 0 (0.000)intergenic (+61/+675) ycdU/serX uncharacterized protein YcdU/tRNA‑Ser

CCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATT  >  NC_000913/1048565‑1049017
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                     
ccACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgact                                                                                                                                                                                                             <  1:292546‑M1/250‑11 (MQ=255)
       gCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGagtgtggactgagtaacc                                                                                                                                                                                                     <  2:55280‑M1/251‑19 (MQ=255)
           aCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagac                                                                                                                                                                                                 <  1:317650‑M1/251‑23 (MQ=255)
                  gtgtCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCCAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCCCTGattgtggactgagtaaccagacaggtgat                                                                                                                                                                                          >  2:664‑M1/1‑222 (MQ=255)
                      cTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCCAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgc                                                                                                                                                                                      >  1:308644‑M1/1‑218 (MQ=255)
                      cTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGACGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgc                                                                                                                                                                                      >  2:115422‑M1/1‑218 (MQ=255)
                                        aTAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTATCATTACGTAAAGTACGCTGGCCATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccct                                                                                                                                                                    >  1:285342‑M1/1‑200 (MQ=255)
                                                                               cATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTGCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatata                                                                                                                              <  2:1793‑M1/250‑90 (MQ=255)
                                                                                aTTACGTAAAGTTCGCTGGCAATGAAATACGACTGTAATTTATTCAGGGTGCGGGTTTAAGAAACGTCGTTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtggtcgggggcacccggtgtgtgcaatatacg                                                                                                                            <  2:249672‑M1/251‑92 (MQ=255)
                                                                                                                    aaTTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttcttt                                                                                        <  1:456893‑M1/251‑128 (MQ=255)
                                                                                                                             aTGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcccgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatac                                                                                                                             >  2:12701‑M1/1‑115 (MQ=255)
                                                                                                                             aTGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatataa                                                                                                                             <  1:12701‑M1/205‑91 (MQ=255)
                                                                                                                               ggTGCTGGTTTAAGACACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGTCTTTGACGAACCTGGCCAGTAGCCGAATCGACCCGGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcttccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcg                                                                              <  1:249018‑M1/250‑138 (MQ=255)
                                                                                                                                   cGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttc                                                                                                        >  1:20703‑M1/1‑109 (MQ=255)
                                                                                                                                           aGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGCATCGAACCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatctacgaaaaacaagcccgtactttcgtacgagctcttctttaaatatggcggggggtggggggg                                                                 >  2:458637‑M1/1‑101 (MQ=255)
                                                                                                                                                 gTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgtggactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggtgcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgaggtggggattgactc                                                           <  2:401230‑M1/251‑157 (MQ=255)
                                                                                                                                                  tCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggggagggggggggtgagtcg                                                          >  1:37636‑M1/1‑94 (MQ=255)
                                                                                                                                                   cATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtttgcaatatacgaaaacaaagcccgtactttcgtccgagctcttctttaaatagggcggttgggggggggagtgctgcg                                                         >  2:258667‑M1/1‑93 (MQ=255)
                                                                                                                                                       tCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgagggggggattgactcgcttc                                                      <  1:115422‑M1/250‑162 (MQ=255)
                                                                                                                                                                  aaCCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCAAGCAGCAGAATCGACCCAGGATGCGTTTTGCCACTGcttgttgcctgagtaaccagacagttgatgtgcacgatttccactcaccctaacagacgtggccgggggcacccggggagtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgcgggggggattgactcacttctctcccgcctcg                                          >  1:222739‑M1/1‑78 (MQ=255)
                                                                                                                                                                                                  gAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgcgggggggggtgacttgctttcttcgcccctcggggcgcccgcccactccggtcacgggctgtcc          >  1:269442‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                        ggCCAGTAGCCGAATCGACCCAGAATGTGGTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgagggggggattgactcgcttcgctcgccctgcgggcatcccgctcactgcgttctcggtctgtccaactgg    <  2:233744‑M1/251‑212 (MQ=255)
                                                                                                                                                                                                          ccAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgagggggggattgactcgcttcgctcgccctgcgggcagcccgctcactgcgttcacggtctgtccaactggct  <  1:107865‑M1/251‑214 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                     
CCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATT  >  NC_000913/1048565‑1049017

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