Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,197,137 25956 bp→CAGATAATTCAA [ymfD][ymgI] 42 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1197137 1223092 25956 208 [0] [0] 204 [ymfD]–[ymgI] [ymfD],ymfE,lit,intE,xisE,ymfH,ymfI,ymfJ,ymfK,ymfT,ymfL,ymfM,oweE,ymfN,aaaE,ymfR,beeE,jayE,ymfQ,ycfK,tfaP,tfaE,stfE,pinE,mcrA,icdC,iraM,ymgK,ymgL,ycgX,bluR,bluF,ycgZ,ymgA,ariR,ymgC,pdeG,ymgF,ycgH,ymgD,ymgG,[ymgI]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 11971360 (0.000)158 (1.410)
+CAGATAATTCAA
62/74 0.0 100% coding (397/666 nt) ymfD e14 prophage; putative SAM‑dependent methyltransferase
?NC_000913 1223093 = 0 (0.000)coding (71/174 nt) ymgI uncharacterized protein YmgI

Only 100 of 204 total aligned reads displayed.

GATCTGCTTAACTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATAAT  >  NC_000913/1223044‑1223142
                                                 |                                                 
ctcctgattttagtaatctcttgatgctaagaactatcagataattcaaa                                                   <  2:4352473‑M2/1‑1 (MQ=255)
 tcctgattttagtaatctcttgatgctaagaactatcagataattcaaAg                                                  >  1:2875739‑M2/49‑50 (MQ=255)
  cctgattttagtaatctcttgatgctaagaactatcagataattcaaAGc                                                 >  1:1765149‑M2/48‑50 (MQ=255)
  cctgattttagtaatctcttgatgctaagaactatcagataattcaaAGc                                                 >  1:2538144‑M2/48‑50 (MQ=255)
   ctgattttagtaatctcttgatgctaagaactatcagataattcaaAg                                                  <  2:5153717‑M2/2‑1 (MQ=255)
   ctgattttagtaatctcttgatgctaagaactatcagataattcaaAGc                                                 <  1:4608238‑M2/3‑1 (MQ=255)
     gattttagtaatctcttgatgctaagaactatcagataattcaaAGCGTg                                              >  1:5246597‑M2/45‑50 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGt                                               <  2:3070627‑M2/5‑1 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGTgg                                             <  2:4759797‑M2/7‑1 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGTg                                              <  1:5689109‑M2/6‑1 (MQ=255)
      attttagtaatctcttgatgctaagaactatcagataattcaaAGCGTg                                              <  1:869629‑M2/6‑1 (MQ=255)
       ttttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGc                                            >  2:5396530‑M2/43‑50 (MQ=255)
        tttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                           >  2:3964491‑M2/42‑50 (MQ=255)
        tttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                           >  1:2637931‑M2/42‑50 (MQ=255)
        tttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                           >  2:3119143‑M2/42‑50 (MQ=255)
        tttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                           >  1:5043016‑M2/42‑50 (MQ=255)
         ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtg                                          >  2:3321560‑M2/41‑50 (MQ=255)
         ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                           <  2:4054652‑M2/9‑1 (MQ=255)
         ttagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCt                                           <  2:5298093‑M2/9‑1 (MQ=255)
          tagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgt                                         >  2:5438999‑M2/40‑50 (MQ=255)
          tagtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgt                                         >  1:5602578‑M2/40‑50 (MQ=255)
            gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgtg                                        <  1:1885341‑M2/12‑1 (MQ=255)
            gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCtgtg                                        <  1:472034‑M2/12‑1 (MQ=255)
            gtaatctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtga                                       <  2:5221113‑M2/13‑1 (MQ=255)
               atctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                    <  2:296690‑M2/16‑1 (MQ=255)
               atctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatg                                     <  2:5649867‑M2/15‑1 (MQ=255)
               atctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatg                                     <  2:3581943‑M2/15‑1 (MQ=255)
                tctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                    >  1:1415061‑M2/34‑49 (MQ=255)
                tctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                    >  1:1534969‑M2/34‑49 (MQ=255)
                 ctcttgatgctaagaactatcagataattcaaAGCGTGGCTGtgatga                                    >  1:956768‑M2/33‑48 (MQ=255)
                  tcttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 >  1:2190752‑M2/32‑50 (MQ=255)
                    ttgttgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                                >  2:3616809‑M2/30‑49 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTgatgat                                   <  1:4843163‑M2/17‑1 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTgatgat                                   <  1:5235109‑M2/17‑1 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTgatgat                                   <  2:3892133‑M2/17‑1 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                                >  1:1307517‑M2/30‑49 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAgg                                >  2:2081849‑M2/30‑49 (MQ=255)
                    ttgatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAg                                 <  2:1226556‑M2/19‑1 (MQ=255)
                      gatgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGaa                              <  2:1716898‑M2/22‑1 (MQ=255)
                       atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAtt                            <  1:3328342‑M2/24‑1 (MQ=255)
                       atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAtt                            <  1:415711‑M2/24‑1 (MQ=255)
                       atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAtt                            <  2:4509198‑M2/24‑1 (MQ=255)
                       atgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAAtt                            <  2:5499304‑M2/24‑1 (MQ=255)
                        tgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTa                           >  1:4986484‑M2/26‑50 (MQ=255)
                        tgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTa                           <  1:3921163‑M2/25‑1 (MQ=255)
                        tgctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTa                           >  1:3710216‑M2/26‑50 (MQ=255)
                         gctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                          >  2:3363996‑M2/25‑50 (MQ=255)
                         gctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                          >  1:538984‑M2/25‑50 (MQ=255)
                          ctaagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTAt                          <  1:5178468‑M2/26‑1 (MQ=255)
                           taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  2:80289‑M2/23‑50 (MQ=255)
                           taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  2:3786685‑M2/23‑50 (MQ=255)
                           taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  1:3429178‑M2/23‑50 (MQ=255)
                           taagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        >  2:600918‑M2/23‑50 (MQ=255)
                            aagaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtt                       >  2:1525621‑M2/22‑50 (MQ=255)
                             agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttt                      <  2:1039357‑M2/30‑1 (MQ=255)
                             agaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGt                        <  1:3590761‑M2/28‑1 (MQ=255)
                              gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     <  1:1478601‑M2/31‑1 (MQ=255)
                              gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     <  1:2383769‑M2/31‑1 (MQ=255)
                              gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     <  2:4788049‑M2/31‑1 (MQ=255)
                              gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     <  1:3978736‑M2/31‑1 (MQ=255)
                              gaactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttt                     >  2:4769111‑M2/20‑50 (MQ=255)
                               aactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttttt                    >  2:4543522‑M2/19‑50 (MQ=255)
                               aactatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttttt                    <  2:2466990‑M2/32‑1 (MQ=255)
                                actatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttttt                   <  1:2239848‑M2/33‑1 (MQ=255)
                                actatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGtttttt                   <  1:4762018‑M2/33‑1 (MQ=255)
                                actatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGttttt                    >  2:2545546‑M2/18‑49 (MQ=255)
                                 ctatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTa                  >  2:3850281‑M2/17‑50 (MQ=255)
                                 ctatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTa                  <  1:4380509‑M2/34‑1 (MQ=255)
                                 ctatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTa                  >  2:1325197‑M2/17‑50 (MQ=255)
                                 ctatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTa                  <  2:2207335‑M2/34‑1 (MQ=255)
                                  tatcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTAc                 >  2:4454247‑M2/16‑50 (MQ=255)
                                  tatcagataattcaaAGCGTGGCTGGGATGATAGGAATTATGTTTTTTAc                 >  2:1972628‑M2/16‑50 (MQ=255)
                                   atcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACg                >  1:2115969‑M2/15‑50 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTAc                 <  1:3167866‑M2/35‑1 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTAc                 <  1:540578‑M2/35‑1 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACg                <  2:5054659‑M2/36‑1 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACg                <  2:183591‑M2/36‑1 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGt               >  1:4473343‑M2/14‑50 (MQ=255)
                                    tcagataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGt               >  1:2983555‑M2/14‑50 (MQ=255)
                                      agataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTTCGTg              <  1:3386879‑M2/38‑1 (MQ=255)
                                       gataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGaa            <  2:4406484‑M2/40‑1 (MQ=255)
                                       gataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGaa            <  1:729143‑M2/40‑1 (MQ=255)
                                       gataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGaa            >  1:4632748‑M2/11‑50 (MQ=255)
                                        ataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           <  2:4654809‑M2/41‑1 (MQ=255)
                                        ataattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           <  2:2870057‑M2/41‑1 (MQ=255)
                                          aattcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           >  2:4159388‑M2/8‑48 (MQ=255)
                                           attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           >  2:2360715‑M2/7‑47 (MQ=255)
                                           attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAAt           >  1:2934064‑M2/7‑47 (MQ=255)
                                           attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATgag        <  1:3723495‑M2/44‑1 (MQ=255)
                                           attcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATga         <  1:242468‑M2/43‑1 (MQ=255)
                                            ttcaaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATgaga       <  2:3683699‑M2/45‑1 (MQ=255)
                                              caaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaa      >  2:4591146‑M2/4‑49 (MQ=255)
                                               aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGataata    <  2:1451248‑M2/48‑1 (MQ=255)
                                               aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaata    <  1:1268268‑M2/48‑1 (MQ=255)
                                               aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaata    <  2:4893882‑M2/48‑1 (MQ=255)
                                               aaAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaata    <  1:1777455‑M2/48‑1 (MQ=255)
                                                aAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaataa   >  2:32943‑M2/2‑50 (MQ=255)
                                                aAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaataa   >  2:1279326‑M2/2‑50 (MQ=255)
                                                aAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaataa   <  1:1189410‑M2/49‑1 (MQ=255)
                                                 aGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGaataat  <  2:761838/50‑1 (MQ=255)
                                                 |                                                 
GATCTGCTTAACTAAAGAACGCTCGCTATTATTCAGATAATTCAAAATGAGCGTGGCTGTGATGATAGGAATTATGTTTTTTACGTGAATGAGAATAAT  >  NC_000913/1223044‑1223142

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