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,0860CT92.6% 490.4 / 23.6 152intergenic (‑300/‑77)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (4/7);  new base T (81/57);  total (88/64)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

Only 100 of 152 total aligned reads displayed.

TTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGA  >  NC_000913/12038‑12134
                                                |                                                
 ttAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAg                                                >  1:5327496/1‑50 (MQ=255)
  tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAg                                                <  1:1873144/49‑1 (MQ=255)
  tagtagTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGt                                               >  2:807940/1‑50 (MQ=255)
   agtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGt                                               <  2:1354339/49‑1 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                              <  2:1918878/49‑1 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                              <  2:2827758/49‑1 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                              >  1:2557152/1‑49 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                             <  2:2663578/50‑1 (MQ=255)
    gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                             <  2:2896743/50‑1 (MQ=255)
     tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                             >  2:2924279/1‑49 (MQ=255)
     tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                             >  1:1520081/1‑49 (MQ=255)
     tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                             >  1:1093126/1‑49 (MQ=255)
      agtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGa                                            >  2:2724486/1‑49 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaaa                                          <  1:4795175/50‑1 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaaa                                          <  2:5461854/50‑1 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaaa                                          <  1:5025878/50‑1 (MQ=255)
       gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaa                                           <  1:217612/49‑1 (MQ=255)
        tCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                         >  2:3772764/1‑50 (MQ=255)
          aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAcc                                        <  2:5121765/49‑1 (MQ=255)
           aCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAg                                      <  1:3854947/50‑1 (MQ=255)
           aCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAg                                      <  2:1584014/50‑1 (MQ=255)
            ccGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGa                                     >  1:4382743/1‑50 (MQ=255)
             cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                    >  2:1680826/1‑50 (MQ=255)
             cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGAc                                    <  2:4373285/50‑1 (MQ=255)
              gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACg                                   >  2:2649403/1‑50 (MQ=255)
               cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGt                                  >  2:662490/1‑50 (MQ=255)
               cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGt                                  <  2:2170904/50‑1 (MQ=255)
               cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGt                                  >  1:4140607/1‑50 (MQ=255)
                agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                 >  1:3296590/1‑50 (MQ=255)
                agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                 <  2:548902/50‑1 (MQ=255)
                agtgagtgAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGtt                                 >  1:2616699/1‑50 (MQ=255)
                 gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGttt                                >  1:4237644/1‑50 (MQ=255)
                 gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGttt                                <  2:3382244/50‑1 (MQ=255)
                  tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                               <  2:2740287/50‑1 (MQ=255)
                  tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                               >  1:968791/1‑50 (MQ=255)
                  tgagtgagTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTc                               <  1:5832510/50‑1 (MQ=255)
                     gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                            >  1:84715/1‑50 (MQ=255)
                     gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                            <  1:858218/50‑1 (MQ=255)
                     gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                            <  1:2448527/50‑1 (MQ=255)
                     gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                            <  2:3547500/50‑1 (MQ=255)
                     gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                            <  1:4167282/50‑1 (MQ=255)
                     gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                            >  2:251029/1‑50 (MQ=255)
                     gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                            >  2:4678476/1‑50 (MQ=255)
                      tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                           <  2:393380/50‑1 (MQ=255)
                      tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                           >  2:2932925/1‑50 (MQ=255)
                      tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                           >  1:2901098/1‑50 (MQ=255)
                      tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                           >  1:3534762/1‑50 (MQ=255)
                      tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                           <  1:889145/50‑1 (MQ=255)
                        agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                          >  2:4205765/1‑49 (MQ=255)
                        agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCt                         <  2:523388/50‑1 (MQ=255)
                        agtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGcccc                          >  1:570146/1‑49 (MQ=255)
                         gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                        >  1:2132609/1‑50 (MQ=255)
                         gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                        >  1:3123642/1‑50 (MQ=255)
                         gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                        <  2:5056841/50‑1 (MQ=255)
                          tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                       >  1:5176082/1‑50 (MQ=255)
                          tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                       >  2:5611639/1‑50 (MQ=255)
                          tCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                       <  2:5380528/50‑1 (MQ=255)
                          tCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTAt                       >  1:5102740/1‑50 (MQ=255)
                           cTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAt                       >  1:3186832/1‑49 (MQ=255)
                            tGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     >  1:4376406/1‑50 (MQ=255)
                             gtaaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     >  2:3297380/3‑49 (MQ=39)
                             gCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     >  2:1807976/1‑49 (MQ=255)
                              cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     <  2:5186174/48‑1 (MQ=255)
                              cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     <  2:60633/48‑1 (MQ=255)
                               aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     >  1:5315042/1‑47 (MQ=255)
                               aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                     >  2:4149840/1‑47 (MQ=255)
                               aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAg                  <  2:4386399/50‑1 (MQ=255)
                                aaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGa                 >  2:1252105/1‑50 (MQ=255)
                                aaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGa                 >  2:477252/1‑50 (MQ=255)
                                 aaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc                >  2:3110689/1‑50 (MQ=255)
                                 aaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc                >  2:1126639/1‑50 (MQ=255)
                                 aaaaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc                <  2:4507383/50‑1 (MQ=255)
                                  aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               >  1:1365875/1‑50 (MQ=255)
                                  aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               >  1:2358708/1‑50 (MQ=255)
                                  aaaaTGAAATTGGGGTGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               >  2:736564/1‑50 (MQ=255)
                                   aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt               >  2:4235923/1‑49 (MQ=255)
                                   aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTc              <  1:365540/50‑1 (MQ=255)
                                    aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca             >  1:2876322/1‑50 (MQ=255)
                                    aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACCca             >  1:3710804/1‑50 (MQ=255)
                                    aaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca             >  2:5325058/1‑50 (MQ=255)
                                     aTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcac            <  2:2374771/50‑1 (MQ=255)
                                     aTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcac            <  2:3213742/50‑1 (MQ=255)
                                      tGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcaca           <  1:2285190/50‑1 (MQ=255)
                                       gAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa          >  2:5700856/1‑50 (MQ=255)
                                        aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         <  1:4482101/50‑1 (MQ=255)
                                        aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:3666018/1‑50 (MQ=255)
                                        aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         <  2:3719308/50‑1 (MQ=255)
                                        aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:4102430/1‑50 (MQ=255)
                                        aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:3892459/1‑50 (MQ=255)
                                         aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc        <  2:3864398/50‑1 (MQ=255)
                                         aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc        <  1:5269586/50‑1 (MQ=255)
                                         aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc        <  2:292424/50‑1 (MQ=255)
                                         aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:5090135/1‑49 (MQ=255)
                                         aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:1212952/1‑49 (MQ=255)
                                         aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:490451/1‑49 (MQ=255)
                                         aaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc         >  2:1380188/1‑49 (MQ=255)
                                          aTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca       <  2:1150090/50‑1 (MQ=255)
                                           ttGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcac      <  1:2250518/50‑1 (MQ=255)
                                            tGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca     >  1:3217558/1‑50 (MQ=255)
                                               gCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatga  >  2:2954838/1‑50 (MQ=255)
                                                |                                                
TTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGA  >  NC_000913/12038‑12134

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: