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 105 [0] [0] 100 [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)80 (1.280)
+CAGATAATTCAA
51/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 105 total aligned reads displayed.

ATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAA  >  NC_000913/1197087‑1197185
                                                 |                                                 
aTAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   <  1:3363154/50‑1 (MQ=255)
 tAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:1147367/1‑49 (MQ=255)
 tAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:4643844/1‑49 (MQ=255)
 tAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:3142596/1‑49 (MQ=255)
  aaGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:5442413/1‑48 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTAt                                                   >  1:3440336/1‑47 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcag                                                <  1:1204204‑M1/50‑4 (MQ=255)
   agagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATca                                                 <  1:2602067‑M1/49‑3 (MQ=255)
     agTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagat                                              <  1:6315183‑M1/50‑6 (MQ=255)
     agTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcaga                                               >  1:359130‑M1/1‑45 (MQ=255)
      gTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagat                                              >  1:1811100‑M1/1‑44 (MQ=255)
       tctcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataa                                            <  1:3830662‑M1/50‑8 (MQ=255)
       tctcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataa                                            <  1:1924973‑M1/50‑8 (MQ=255)
         tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:3743543‑M1/1‑41 (MQ=255)
         tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:2662134‑M1/1‑41 (MQ=255)
          cACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:1979756‑M1/1‑40 (MQ=255)
          cACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:1454327‑M1/1‑40 (MQ=255)
          cACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          >  1:3008677‑M1/1‑40 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattc                                         >  1:2242669‑M1/1‑39 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattc                                         <  1:4933555‑M1/49‑11 (MQ=255)
           aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataatt                                          <  1:1507361‑M1/48‑10 (MQ=255)
            cctcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaa                                       <  1:4236453‑M1/50‑13 (MQ=255)
             ctcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaa                                      <  1:4710632‑M1/50‑14 (MQ=255)
             ctcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattca                                        >  1:1311064‑M1/1‑37 (MQ=255)
             ctcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattca                                        >  1:3187990‑M1/1‑37 (MQ=255)
              tcctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaag                                     >  1:6154771‑M1/1‑36 (MQ=255)
                ctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcg                                   >  1:1271249‑M1/1‑34 (MQ=255)
                ctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcg                                   <  1:5207414‑M1/50‑17 (MQ=255)
                ctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagc                                    <  1:1058647‑M1/49‑16 (MQ=255)
                ctGATTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaag                                     <  1:4906623‑M1/48‑15 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtg                                 <  1:674984‑M1/49‑19 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtg                                 <  1:3964760‑M1/49‑19 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtg                                 <  1:88085‑M1/49‑19 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgt                                  <  1:3604212‑M1/48‑18 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgt                                  <  1:6382420‑M1/48‑18 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgt                                  <  1:5728758‑M1/48‑18 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgt                                  <  1:3522543‑M1/48‑18 (MQ=255)
                   aTTTTAGTAATCTCTTGATGCTAAGAACGATcagataattcaaagcgtg                                 <  1:4791378‑M1/49‑19 (MQ=255)
                    ttttAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggc                               >  1:4968125‑M1/1‑30 (MQ=255)
                     tttAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggct                              >  1:3029448‑M1/1‑29 (MQ=255)
                     tttAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggct                              >  1:6339288‑M1/1‑29 (MQ=255)
                       tAGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgt                            >  1:4981964‑M1/1‑27 (MQ=255)
                        aGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgt                            <  1:223109‑M1/49‑24 (MQ=255)
                        aGTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgt                            <  1:5795860‑M1/49‑24 (MQ=255)
                         gTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtga                          <  1:6223955‑M1/50‑26 (MQ=255)
                         gTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtga                          >  1:542936‑M1/1‑25 (MQ=255)
                         gTAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtg                           <  1:4337203‑M1/49‑25 (MQ=255)
                          tAATCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgat                         >  1:5246840‑M1/1‑24 (MQ=255)
                           aaTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgat                         >  1:4274565‑M1/1‑23 (MQ=255)
                           aaTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgat                         >  1:1679042‑M1/1‑23 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:5806847‑M1/1‑22 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       <  1:3767919‑M1/50‑29 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       <  1:2118253‑M1/50‑29 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       <  1:3110062‑M1/50‑29 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:1678813‑M1/1‑22 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatg                        <  1:228808‑M1/49‑28 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatg                        <  1:4886284‑M1/49‑28 (MQ=255)
                            aTCTCTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatg                        <  1:3078670‑M1/49‑28 (MQ=255)
                             tctcTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:2924336‑M1/1‑21 (MQ=255)
                             tctcTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:1450080‑M1/1‑21 (MQ=255)
                             tctcTTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:4707136‑M1/1‑21 (MQ=255)
                              ctctTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatga                       >  1:2509035‑M1/1‑20 (MQ=255)
                                ctTGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    >  1:6291627‑M1/1‑18 (MQ=255)
                                 ttGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    <  1:1817079‑M1/48‑32 (MQ=255)
                                 ttGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    <  1:805270‑M1/48‑32 (MQ=255)
                                 ttGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatag                    <  1:2835716‑M1/48‑32 (MQ=255)
                                  tGATGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaa                 >  1:3961494‑M1/1‑16 (MQ=255)
                                   gATGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaat                >  1:4094623‑M1/1‑15 (MQ=255)
                                   gATGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaat                >  1:2388752‑M1/1‑15 (MQ=255)
                                   gATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatagga                  <  1:644841‑M1/48‑34 (MQ=255)
                                   gATGCTAAGAACTATcagataattcaaagcgtggctgtgatgatagg                   <  1:2202152‑M1/47‑33 (MQ=255)
                                    aTGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaatt               <  1:4094796‑M1/50‑37 (MQ=255)
                                    aTGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaat                >  1:3596411‑M1/1‑14 (MQ=255)
                                     tGCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaatta              >  1:977478‑M1/1‑13 (MQ=255)
                                      gCTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattat             >  1:1641301‑M1/1‑12 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattatg            <  1:37891‑M1/50‑40 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattat             <  1:4981985‑M1/49‑39 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattat             <  1:246770‑M1/49‑39 (MQ=255)
                                       cTAAGAACTATcagataattcaaagcgtggctgtgatgataggaattat             <  1:6194644‑M1/49‑39 (MQ=255)
                                        tAAGAACTATcagataattcaaagcgtggctgtgatgataggaattatgt           >  1:3143416‑M1/1‑10 (MQ=255)
                                         aaGAACTATcagataattcaaagcgtggctgtgatgataggaattatgtt          >  1:1586113‑M1/1‑9 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgttt         <  1:4511425‑M1/50‑43 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgttt         >  1:4864450‑M1/1‑8 (MQ=255)
                                          aGAACTATcagataattcaaagcgtggctgtgatgataggaattatgttt         <  1:2856873‑M1/50‑43 (MQ=255)
                                           gAACTATcagataattcaaagcgtggctgtgatgataggaattatgtttt        <  1:2110825‑M1/50‑44 (MQ=255)
                                           gAACTATcagataattcaaagcgtggctgtgatgataggaattatgtttt        >  1:4861327‑M1/1‑7 (MQ=255)
                                           gAACTATcagataattcaaagcgtggctgtgatgataggaattatgtttt        <  1:716719‑M1/50‑44 (MQ=255)
                                            aaCTATcagataattcaaagcgtggctgtgatgataggaattatgttttt       <  1:6153998‑M1/50‑45 (MQ=255)
                                             aCTATcagataattcaaagcgtggctgtgatgctaggaattatgttttt       >  1:5132655‑M1/1‑5 (MQ=255)
                                             aCTATcagataattcaaagcgtggctgtgatgataggaattatgtttttt      <  1:2538841‑M1/50‑46 (MQ=255)
                                             aCTATcagataattcaaagcgtggctgtgatgataggaattatgtttttt      <  1:3536012‑M1/50‑46 (MQ=255)
                                             aCTATcagataattcaaagcgtggctgtgatgataggaattatgttttt       >  1:762597‑M1/1‑5 (MQ=255)
                                              cTATcagataattcaaagcgtggctgtgatgataggaattatgtttttta     >  1:3371060‑M1/1‑4 (MQ=255)
                                              cTATcagataattcaaagcgtggctgtgatgataggaattatgtttttta     >  1:3717450‑M1/1‑4 (MQ=255)
                                               tATcagataattcaaagcgtggctgtgatgataggaattatgttttttac    >  1:970105‑M1/1‑3 (MQ=255)
                                                aTcagataattcaaagcgtggctgtgatgataggaattatgttttttacg   >  1:1525135‑M1/1‑2 (MQ=255)
                                                 tcagataattcaaagcgtggctgtgatgataggaattatgttttttacgt  <  1:689571‑M1/50‑50 (MQ=255)
                                                 tcagataattcaaagcgtggctgtgatgataggaattatgttttttacg   <  1:413842‑M1/49‑49 (MQ=255)
                                                 tcagataattcaaagcgtggctgtgatgataggaattatgttttttac    <  1:5025907‑M1/48‑48 (MQ=255)
                                                 tcagataattcaaagcgtggctgtgatgataggaattatgttttttac    <  1:2957071‑M1/48‑48 (MQ=255)
                                                 |                                                 
ATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAA  >  NC_000913/1197087‑1197185

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: