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,0860CT81.9% 312.9 / 67.6 127intergenic (‑300/‑77)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (12/11);  new base T (56/48);  total (68/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-03

Only 100 of 127 total aligned reads displayed.

TTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGA  >  NC_000913/12039‑12134
                                               |                                                
 tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAg                                                <  1:2796072/49‑1 (MQ=255)
 tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGt                                               >  2:4926434/1‑50 (MQ=255)
 tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGt                                               >  1:314469/1‑50 (MQ=255)
 tagtacTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGt                                               >  1:3671247/1‑50 (MQ=255)
  agtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                              >  2:1850430/1‑50 (MQ=255)
  agtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGtt                                              >  1:4594277/1‑50 (MQ=255)
  agtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGt                                               <  2:2186212/49‑1 (MQ=255)
   gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                              <  2:1998048/49‑1 (MQ=255)
   gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGtt                                              >  1:1775952/1‑49 (MQ=255)
   gtagtCAACCGCAATGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                             <  1:2476055/50‑1 (MQ=255)
    tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                             >  1:1851285/1‑49 (MQ=255)
     agtCAACCGCAGCGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGa                                            >  1:1234783/1‑49 (MQ=255)
      gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaaa                                          >  1:4093068/1‑50 (MQ=255)
      gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaa                                           <  1:4966201/49‑1 (MQ=255)
      gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaa                                           <  2:4146166/49‑1 (MQ=255)
       tCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                         >  2:1679472/1‑50 (MQ=255)
       tCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAAc                                         >  1:2774623/1‑50 (MQ=255)
       tCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAAc                                         <  1:1707136/50‑1 (MQ=255)
        cAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                         >  2:195027/1‑49 (MQ=255)
         aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                         >  2:449241/1‑48 (MQ=255)
         aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                         >  2:500444/1‑48 (MQ=255)
         aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCa                                       <  2:3201938/50‑1 (MQ=255)
          aCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAg                                      <  2:312673/50‑1 (MQ=255)
          aCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAg                                      <  2:4171097/50‑1 (MQ=255)
           ccGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGa                                     >  2:3794206/1‑50 (MQ=255)
           ccGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGa                                     <  1:5126581/50‑1 (MQ=255)
           ccGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGa                                     <  2:1653013/50‑1 (MQ=255)
           ccGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGa                                     <  2:3547567/50‑1 (MQ=255)
            cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                    >  2:2502/1‑50 (MQ=255)
            cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                    <  2:1107919/50‑1 (MQ=255)
            cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                    <  2:3285370/50‑1 (MQ=255)
              cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGt                                  >  1:1214462/1‑50 (MQ=255)
               agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAATCCAGACGtt                                 <  2:5188045/50‑1 (MQ=39)
               agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                 >  2:2640686/1‑50 (MQ=255)
               agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                 >  2:149946/1‑50 (MQ=255)
               agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                 >  1:810537/1‑50 (MQ=255)
               agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                 >  2:2744052/1‑50 (MQ=255)
               agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                 <  1:729137/50‑1 (MQ=255)
                gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGttt                                <  2:748293/50‑1 (MQ=255)
                gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGttt                                <  1:3267261/50‑1 (MQ=255)
                 tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                               >  1:442935/1‑50 (MQ=255)
                 tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                               <  1:3453920/50‑1 (MQ=255)
                 tgagtgagTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTc                               >  1:660005/1‑50 (MQ=255)
                  gagtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAACTCAGCCGTTTCg                              >  1:2318272/1‑50 (MQ=255)
                  gagtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCg                              >  1:972224/1‑50 (MQ=255)
                   agtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCg                              <  1:4930563/49‑1 (MQ=255)
                     tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                           >  1:3398927/1‑50 (MQ=255)
                      gagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                          >  1:4282126/1‑50 (MQ=255)
                      gagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                          >  1:2298096/1‑50 (MQ=255)
                      gagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                          >  1:4567466/1‑50 (MQ=255)
                      gagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                          >  2:1459120/1‑50 (MQ=255)
                       agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCt                         <  2:2130006/50‑1 (MQ=255)
                        gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                        >  1:5175425/1‑50 (MQ=255)
                        gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                        >  2:528560/1‑50 (MQ=255)
                        gtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTa                        <  1:132931/50‑1 (MQ=255)
                         tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                       <  1:4018799/50‑1 (MQ=255)
                         tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                       <  1:3583750/50‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGTAGTTGAACCCAGACGTTTCGCCCCTATTACa                   <  2:4772781/50‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAtt                      <  1:3647998/47‑1 (MQ=255)
                             cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     >  1:274180/1‑48 (MQ=255)
                             cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACa                   <  2:1501464/50‑1 (MQ=255)
                              aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     >  2:1070263/1‑47 (MQ=255)
                              aaaaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAg                  <  2:1740060/50‑1 (MQ=255)
                              aaaaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAg                  <  1:2082731/50‑1 (MQ=255)
                               aaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGa                 >  2:3533576/1‑50 (MQ=255)
                               aaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGa                 <  1:2581958/50‑1 (MQ=255)
                               aaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGa                 >  1:2149073/1‑50 (MQ=255)
                                aaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc                >  2:5383171/1‑50 (MQ=255)
                                 aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               <  1:5160667/50‑1 (MQ=255)
                                 aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               >  1:3232989/1‑50 (MQ=255)
                                 aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               >  2:668837/1‑50 (MQ=255)
                                 aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               <  1:3122107/50‑1 (MQ=255)
                                 aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               <  1:3738733/50‑1 (MQ=255)
                                  aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTc              <  1:1269385/50‑1 (MQ=255)
                                  aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTc              <  1:5010741/50‑1 (MQ=255)
                                  aaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               >  2:3513107/1‑49 (MQ=255)
                                   aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca             >  2:289450/1‑50 (MQ=255)
                                   aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca             >  1:4129059/1‑50 (MQ=255)
                                   aaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca             >  2:5172336/1‑50 (MQ=255)
                                    aTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca             >  1:935878/1‑49 (MQ=255)
                                    aTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca             >  1:3587316/1‑49 (MQ=255)
                                      gAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa          >  1:3708210/1‑50 (MQ=255)
                                       aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa          <  1:2397407/49‑1 (MQ=255)
                                       aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:746588/1‑50 (MQ=255)
                                       aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         <  2:2670786/50‑1 (MQ=255)
                                       aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:1876801/1‑50 (MQ=255)
                                       aaaTTGGGTAGTTGAAACCAGACGATTCGCCCCTATTACAGACTCACAAc         >  1:4466204/1‑50 (MQ=255)
                                       aaaTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:1141751/1‑50 (MQ=255)
                                        aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc        <  1:3554758/50‑1 (MQ=255)
                                        aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:4333557/1‑49 (MQ=255)
                                        aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:5178389/1‑49 (MQ=255)
                                        aaTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc        <  2:4287041/50‑1 (MQ=255)
                                        aaTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc        <  2:2314683/50‑1 (MQ=255)
                                         aTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca       >  1:3946324/1‑50 (MQ=255)
                                         aTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca       >  1:3532731/1‑50 (MQ=255)
                                         aTTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca       <  2:3902823/50‑1 (MQ=255)
                                           tGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca     >  1:2775902/1‑50 (MQ=255)
                                            gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca     <  2:3914842/49‑1 (MQ=255)
                                            gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACat    >  1:3795730/1‑50 (MQ=255)
                                             ggCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatg   <  1:1526290/50‑1 (MQ=255)
                                               |                                                
TTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGA  >  NC_000913/12039‑12134

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