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% 149.7 / NA 39intergenic (‑300/‑77)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/25);  total (14/25)

GTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACA  >  NC_000913/12042‑12131
                                            |                                             
gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                           <  1:817114/49‑1 (MQ=255)
   gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaa                                        <  1:939039/49‑1 (MQ=255)
   gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaa                                        <  1:301801/49‑1 (MQ=255)
   gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGa                                         <  1:1644006/48‑1 (MQ=255)
      aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                      <  1:1322005/48‑1 (MQ=255)
      aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                      <  1:998855/48‑1 (MQ=255)
      aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                      <  1:801877/48‑1 (MQ=255)
         cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGa                                  <  1:1310809/49‑1 (MQ=255)
         cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGa                                  <  1:2038521/49‑1 (MQ=255)
          gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                 >  1:652224/1‑49 (MQ=255)
          gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                 <  1:723597/49‑1 (MQ=255)
           cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACg                                <  1:639079/49‑1 (MQ=255)
            agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGt                               >  1:19084/1‑49 (MQ=255)
             gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                              <  1:590491/49‑1 (MQ=255)
              tgagtgagTCTGCAAAAAAATGAGATTGGGTAGTTGAAACCAGACGttt                             <  1:831068/49‑1 (MQ=255)
                agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                            <  1:1956332/48‑1 (MQ=255)
                agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                            <  1:275856/48‑1 (MQ=255)
                agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                            <  1:482396/48‑1 (MQ=255)
                agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCg                           >  1:1844315/1‑49 (MQ=255)
                  tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                         <  1:2046049/49‑1 (MQ=255)
                  tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                         >  1:137722/1‑49 (MQ=255)
                    agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                       >  1:181615/1‑49 (MQ=255)
                     gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCt                      <  1:598755/49‑1 (MQ=255)
                     gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCt                      <  1:1686590/49‑1 (MQ=255)
                      tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                     >  1:1068638/1‑49 (MQ=255)
                       cTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                    >  1:1160110/1‑49 (MQ=255)
                          cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                    <  1:284670/46‑1 (MQ=255)
                          cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                  >  1:95250/1‑48 (MQ=255)
                           aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACa                <  1:650339/49‑1 (MQ=255)
                             aaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGa              <  1:534416/49‑1 (MQ=255)
                               aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt            >  1:663146/1‑49 (MQ=255)
                               aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt            >  1:728132/1‑49 (MQ=255)
                                    aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca        <  1:1351502/48‑1 (MQ=255)
                                    aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa       <  1:539651/49‑1 (MQ=255)
                                    aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa       >  1:865884/1‑49 (MQ=255)
                                      aTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc     <  1:321635/49‑1 (MQ=255)
                                        tGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca    >  1:371795/1‑48 (MQ=255)
                                         gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca  >  1:328436/1‑49 (MQ=255)
                                         gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca  >  1:1723550/1‑49 (MQ=255)
                                            |                                             
GTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACA  >  NC_000913/12042‑12131

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: