Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 324,001 C→T 12.3% T138I (ACC→ATC)  htpG → molecular chaperone HSP90 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE324,0010CT12.3% 231.9 / 20.1 81T138I (ACC→ATC) htpGmolecular chaperone HSP90 family
Reads supporting (aligned to +/- strand):  ref base C (25/46);  new base T (7/3);  total (32/49)
Fisher's exact test for biased strand distribution p-value = 4.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

AAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAA  >  minE/323932‑324071
                                                                     |                                                                      
aaGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGAcc                                                                       <  1:1083663/71‑1 (MQ=255)
aaGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGAc                                                                        >  1:2471794/1‑70 (MQ=255)
aaGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGAc                                                                        >  1:1454971/1‑70 (MQ=255)
  gACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGt                                                                     <  1:3724489/71‑1 (MQ=255)
  gACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGt                                                                     <  1:717206/71‑1 (MQ=255)
  gACAGCCAGCTGATCGGTCAGTTTGGAGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGt                                                                     <  1:2689722/71‑1 (MQ=255)
   aCAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTg                                                                    >  1:321531/1‑71 (MQ=255)
    tagccagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGATCGTGc                                                                   >  1:879725/2‑71 (MQ=255)
    tagccagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGATCGTGc                                                                   >  1:874651/2‑71 (MQ=255)
       ccagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTa                                                                >  1:895070/1‑71 (MQ=255)
        cagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGAc                                                                        >  1:746616/1‑62 (MQ=255)
        cagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTAc                                                               >  1:1392314/1‑71 (MQ=255)
           cTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTcg                                                            >  1:1755949/1‑71 (MQ=255)
            tGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTcgc                                                           >  1:1013247/1‑71 (MQ=255)
            tGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTcgc                                                           >  1:180776/1‑71 (MQ=255)
            tGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTcgc                                                           >  1:3545207/1‑71 (MQ=255)
             gATCGGGCAGTTTGGTGTTGGGTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGAcc                                                                       <  1:1943708/58‑1 (MQ=255)
                 ggTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGATCGTGCGTACTCGCGCgg                                                       >  1:3618582/1‑70 (MQ=255)
                 ggTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGATCGTGCGTACTCGCGCgg                                                       >  1:3663737/1‑70 (MQ=255)
                  gTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCa                                                     <  1:1396134/71‑1 (MQ=255)
                  gTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCa                                                     <  1:1206334/71‑1 (MQ=255)
                  gTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGGACTCGCGCGGCa                                                     <  1:1488274/71‑1 (MQ=255)
                  gTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCCTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCa                                                     <  1:1944171/71‑1 (MQ=255)
                   tCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTcgcg                                                          <  1:1934966/65‑1 (MQ=255)
                   tCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTcgcg                                                          <  1:3680030/65‑1 (MQ=255)
                    cAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTcgcg                                                          <  1:3684728/64‑1 (MQ=255)
                      gTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGc                                                      >  1:2431946/1‑66 (MQ=255)
                       tttGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGa                                                <  1:2835781/71‑1 (MQ=255)
                       tttGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGa                                                <  1:797425/71‑1 (MQ=255)
                         tGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGaa                                               <  1:3188702/70‑1 (MQ=255)
                          ggTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGaaa                                              <  1:1994368/70‑1 (MQ=255)
                          ggTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGaaa                                              <  1:3163670/70‑1 (MQ=255)
                           gtgtTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGATCGTGCGTACTCGCGCGGCAGGCGaaaaa                                            <  1:755788/71‑1 (MQ=255)
                           gtgtTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGaaaaa                                            <  1:3740626/71‑1 (MQ=255)
                            tgtTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGaaa                                              <  1:1984206/68‑1 (MQ=255)
                              ttGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAAc                                           >  1:798107/1‑69 (MQ=255)
                                ggTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGATCGTGCGTACTCGCGCGGCAGGCGAAAAACAAg                                        >  1:3321398/1‑70 (MQ=255)
                                ggTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAg                                        >  1:2033130/1‑70 (MQ=255)
                                 gTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTa                                                                >  1:1045021/1‑45 (MQ=255)
                                 gTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTa                                                                >  1:2392519/1‑45 (MQ=255)
                                 gTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaa                                      <  1:2786219/71‑1 (MQ=255)
                                    tCTACTCTGCGTTTATCGTGGCCGAGAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGa                                       <  1:1812649/67‑1 (MQ=255)
                                    tCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaaaa                                    <  1:2126242/70‑1 (MQ=255)
                                    tCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAAt                                   >  1:3305377/1‑71 (MQ=255)
                                      tACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACt                                                              >  1:2243924/1‑42 (MQ=255)
                                        ctctGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATg                                  <  1:3061758/68‑1 (MQ=255)
                                        ctctGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATg                                  <  1:1022999/68‑1 (MQ=255)
                                        ctctGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGc                                <  1:3390496/70‑1 (MQ=255)
                                           ttcGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtct                            >  1:2248887/3‑71 (MQ=255)
                                             cGTTTATCGTGTCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtcttc                          <  1:2178232/71‑1 (MQ=255)
                                             cGTTTATCGTGGCCGACAAAGTGATCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtcttc                          <  1:2868356/71‑1 (MQ=255)
                                             cGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATg                                  <  1:575984/63‑1 (MQ=255)
                                             cGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtcttc                          <  1:1105125/71‑1 (MQ=255)
                                             cGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtcttc                          >  1:3121109/1‑71 (MQ=255)
                                               tttATCGTGGCCGACAAAGTGATCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtcttct                         <  1:2347428/70‑1 (MQ=255)
                                               tttATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaaaa                                    <  1:1656973/59‑1 (MQ=255)
                                                 tATCGTGGCCGACAAAGTGATCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTgg                       >  1:2016130/1‑70 (MQ=255)
                                                 tATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtcttct                         <  1:2360979/68‑1 (MQ=255)
                                                   tCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtcttct                         >  1:1105128/1‑66 (MQ=255)
                                                   tCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtcttc                          >  1:1961253/1‑65 (MQ=255)
                                                   tCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGa                     <  1:1135463/70‑1 (MQ=255)
                                                     gTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaa                                      >  1:2741785/1‑51 (MQ=255)
                                                     gTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaa                                      >  1:3115552/1‑51 (MQ=255)
                                                      tGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCg                 >  1:1035354/1‑71 (MQ=255)
                                                      aggCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtct                            <  1:422659/59‑1 (MQ=255)
                                                        gCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGc               <  1:3149431/71‑1 (MQ=255)
                                                        gCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGc               <  1:3144713/71‑1 (MQ=255)
                                                         ccGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCt              <  1:370310/71‑1 (MQ=255)
                                                         ccGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCt              <  1:1050401/71‑1 (MQ=255)
                                                          cGACAAAGTGACCGTGCGTGCTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTg             >  1:3235992/1‑71 (MQ=255)
                                                           gACAAAGTGATCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAAt                   >  1:1156613/1‑64 (MQ=255)
                                                           gACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCg                               <  1:35796/52‑1 (MQ=255)
                                                           gACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCg                               <  1:3262282/52‑1 (MQ=255)
                                                            aCAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaaaa                                    <  1:3573857/46‑1 (MQ=255)
                                                            aCAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaa                                      <  1:1682463/44‑1 (MQ=255)
                                                            aCAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaa                                      <  1:3464251/44‑1 (MQ=255)
                                                            aCAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGaa                                      <  1:1623541/44‑1 (MQ=255)
                                                            aCAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCg                               <  1:2969452/51‑1 (MQ=255)
                                                            aCAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGtc                             <  1:2753042/53‑1 (MQ=255)
                                                              aaaGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATg                                  >  1:3821342/1‑46 (MQ=255)
                                                              aaaGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATg                                  >  1:3116006/1‑46 (MQ=255)
                                                              aaaGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGa         <  1:2072524/71‑1 (MQ=255)
                                                              aaaGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGa         <  1:3178305/71‑1 (MQ=255)
                                                                aGTGACCGTGAGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAg       >  1:1242697/1‑71 (MQ=255)
                                                                  tGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGt     <  1:1730150/71‑1 (MQ=255)
                                                                  tGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGt     <  1:2197916/71‑1 (MQ=255)
                                                                     ccGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGaa  <  1:765725/71‑1 (MQ=255)
                                                                     |                                                                      
AAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAA  >  minE/323932‑324071

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