Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 12,086 C→T intergenic (‑300/‑77) yaaI ← / → dnaK DUF2541 domain‑containing protein YaaI/chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,0860CT100.0% 122.6 / NA 32intergenic (‑300/‑77)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/15);  total (17/15)

GTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCAC  >  NC_000913/12042‑12130
                                            |                                            
gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGt                                           <  1:995422/48‑1 (MQ=255)
 tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                         <  1:750807/49‑1 (MQ=255)
    tCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaaa                                      <  1:2858820/49‑1 (MQ=255)
      aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAcc                                    <  1:85450/49‑1 (MQ=255)
      aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAcc                                    <  1:261547/49‑1 (MQ=255)
       aCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                     >  1:2379431/1‑47 (MQ=255)
       aCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                     >  1:55750/1‑47 (MQ=255)
            agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGt                              <  1:1762189/49‑1 (MQ=255)
             gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                             >  1:1896539/1‑49 (MQ=255)
               gagtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                           >  1:850472/1‑49 (MQ=255)
               gagtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                           >  1:517613/1‑49 (MQ=255)
                agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                           <  1:934032/48‑1 (MQ=255)
                agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                           <  1:1377561/48‑1 (MQ=255)
                agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCg                          >  1:2146569/1‑49 (MQ=255)
                  tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                        >  1:2166946/1‑49 (MQ=255)
                  tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                        >  1:1449617/1‑49 (MQ=255)
                   gagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                       >  1:508956/1‑49 (MQ=255)
                    agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                      >  1:1100150/1‑49 (MQ=255)
                    agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                      <  1:342873/49‑1 (MQ=255)
                      tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                    >  1:2760309/1‑49 (MQ=255)
                      tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                    <  1:944332/49‑1 (MQ=255)
                       cTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                   >  1:2268162/1‑49 (MQ=255)
                       cTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                   >  1:1798832/1‑49 (MQ=255)
                         gCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                 >  1:2460134/1‑49 (MQ=255)
                           aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACa               <  1:2581534/49‑1 (MQ=255)
                                aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTc          <  1:1225849/49‑1 (MQ=255)
                                 aTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca         >  1:454784/1‑49 (MQ=255)
                                 aTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca         <  1:2756999/49‑1 (MQ=255)
                                  tGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca         >  1:305187/1‑48 (MQ=255)
                                   gAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcac        >  1:869255/1‑48 (MQ=255)
                                      aTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc    <  1:293005/49‑1 (MQ=255)
                                         gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcac  <  1:1477505/48‑1 (MQ=255)
                                            |                                            
GTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCAC  >  NC_000913/12042‑12130

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