Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,371,351 C→T 8.1% A109V (GCG→GTG)  groL → Cpn60 chaperonin GroEL, large subunit of GroESL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,371,3510CT 8.1% 130.5 / 3.0 50A109V (GCG→GTG) groLCpn60 chaperonin GroEL, large subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base T (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

CGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTA  >  NC_000913/4371205‑4371499
                                                                                                                                                  |                                                                                                                                                    
cggaCTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCgg                                                                                                                                                    <  2:237864/146‑1 (MQ=255)
       gAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGGCTGAAAGCTGTTGCTGCGGGCATGaa                                                                                                                                             >  1:103844/1‑149 (MQ=255)
       gAAGACAAGTTCGAAAATATGGGTACGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGaa                                                                                                                                             >  1:596696/1‑149 (MQ=255)
         aGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAAcc                                                                                                                                           >  1:784588/1‑149 (MQ=255)
         aGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAAcc                                                                                                                                           >  1:66972/1‑149 (MQ=255)
              aGTGCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATg                                                                                                                                      <  1:575871/149‑1 (MQ=255)
                ttCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGa                                                                                                                                    <  1:657566/149‑1 (MQ=255)
                     aaaTATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGa                                                                                                                               <  1:612399/149‑1 (MQ=255)
                      aaTATGGGGGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGaa                                                                                                                              <  2:118825/149‑1 (MQ=255)
                          cgcgTGCGCAGATGGTGAAAGAAGTCGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGt                                                                                                                          <  1:154993/146‑1 (MQ=255)
                           gggTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTg                                                                                                                         >  1:492830/1‑149 (MQ=255)
                              tGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCAATGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTa                                                                                                                      >  1:503345/1‑149 (MQ=255)
                                    gATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGTGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACa                                                                                                                >  2:264883/1‑149 (MQ=255)
                                    gATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGTGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACa                                                                                                                >  2:662215/1‑149 (MQ=255)
                                           aaaGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTACCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGAACTGTAACGTGGTATCGACAAAGCGGt                                                                                                         >  1:415603/1‑149 (MQ=255)
                                            aagaagTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGtt                                                                                                        <  1:645758/149‑1 (MQ=255)
                                              gaagTGGCATCCAAAGCAAACGACGCTGCAGGCGACAGTACCACCACTGCAACGGCACTGGCTCCGACGGTCATCACCGCAAGGCTGAACGCTGTTGATGCGGGAATGAACCCGATGGACCTGAAACGTGGTAGCGCCAAAGAGGGTAc                                                                                                      >  2:258579/1‑149 (MQ=255)
                                               aagTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACACAGCGGTTAcc                                                                                                     >  1:400102/1‑149 (MQ=255)
                                               aagTTGCCTCTAAAGCAAACGACGCTGCAAGCGACGGTACCACCACTGCAACCGTACTGGCCCAGGCTATCATCCCTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACCTGGCATCGACAAAGCGGTTAcc                                                                                                     >  2:587868/1‑149 (MQ=255)
                                                    gCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGc                                                                                                >  1:374958/1‑149 (MQ=255)
                                                    gCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGAAGGACCCGAAACGTGGTCTGGACACAGCGGTTACCGCTGc                                                                                                >  2:77876/1‑149 (MQ=255)
                                                    gCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGACTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTTGTAACGACAAAGCGGTTCCCGCtcc                                                                                                >  1:645130/1‑147 (MQ=255)
                                                     ccTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGAGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCa                                                                                               <  1:424222/149‑1 (MQ=255)
                                                           aaGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAGGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCCGCAAAGCGGTTACCGCTGCAGTTgaa                                                                                         >  1:490220/1‑149 (MQ=255)
                                                              cAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTgaagaa                                                                                      >  1:227628/1‑149 (MQ=255)
                                                                  cgacgCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGa                                                                                  <  1:714436/149‑1 (MQ=255)
                                                                        tGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAgcgc                                                                            >  1:518865/1‑149 (MQ=255)
                                                                         gCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGTGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCt                                                                           <  1:264883/149‑1 (MQ=255)
                                                                         gCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGTGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCt                                                                           <  1:662215/149‑1 (MQ=255)
                                                                           gggCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGt                                                                         <  2:799280/148‑1 (MQ=255)
                                                                           aGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGt                                                                         >  1:741467/1‑149 (MQ=255)
                                                                              cGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCg                                                                      >  1:17650/1‑149 (MQ=255)
                                                                                    tACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTCTCGACAAAGCGGTTACCGCCGCAGTGGAAGAACTGAAAGCGCTGTCCGTACCag                                                                >  2:726803/1‑148 (MQ=255)
                                                                                     accaccACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATg                                                               <  2:17650/149‑1 (MQ=255)
                                                                                           acTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGa                                                         <  1:734313/149‑1 (MQ=255)
                                                                                            gtGCAACCGTACTGGGTCAGGCTATTATCACTGAAGGGCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTGGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGAc                                                        <  2:180173/148‑1 (MQ=255)
                                                                                            cTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGAc                                                        <  2:53044/149‑1 (MQ=255)
                                                                                                  ccGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGACGGACCGTAAACGTGGTATCGGCAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACGTGGCTCTGACCGTaaa                                                  >  1:154551/1‑149 (MQ=255)
                                                                                                     tACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCg                                               >  2:479130/1‑149 (MQ=255)
                                                                                                        tGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGAtt                                            <  1:724386/149‑1 (MQ=255)
                                                                                                         ggCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTg                                           <  1:173637/149‑1 (MQ=255)
                                                                                                          gCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAACGATTGc                                          >  1:221700/1‑149 (MQ=255)
                                                                                                                     catcaCTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTa                               <  2:706919/149‑1 (MQ=255)
                                                                                                                      atccCTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTAc                              <  1:710689/149‑1 (MQ=255)
                                                                                                                       tcaCTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGCTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTTAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTAcc                             <  1:461177/149‑1 (MQ=255)
                                                                                                                              aGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCg                      <  1:80608/149‑1 (MQ=255)
                                                                                                                                   tGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAAc                 >  1:470343/1‑149 (MQ=255)
                                                                                                                                        gCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCcga            >  2:388521/1‑149 (MQ=255)
                                                                                                                                         cTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCcgac           <  2:650426/149‑1 (MQ=255)
                                                                                                                                         cTGTTGCTGCCGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCcgac           <  2:721717/149‑1 (MQ=255)
                                                                                                                                             tgctgcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGaaa       <  2:503345/149‑1 (MQ=255)
                                                                                                                                               ctgcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAAcc     <  2:66972/149‑1 (MQ=255)
                                                                                                                                               ctgcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAAcc     <  2:784588/149‑1 (MQ=255)
                                                                                                                                                 gcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGt   >  1:635607/1‑149 (MQ=255)
                                                                                                                                                  cGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTa  >  1:688987/1‑149 (MQ=255)
                                                                                                                                                  cGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTa  <  2:787307/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTA  >  NC_000913/4371205‑4371499

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