Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,914,888 T→G Q216P (CAG→CCG)  proQ ← RNA chaperone, putative ProP translation regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,914,8880TG98.7% 287.2 / ‑4.5 76Q216P (CAG→CCG) proQRNA chaperone, putative ProP translation regulator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (45/30);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCT  >  NC_000913/1914793‑1914983
                                                                                               |                                                                                               
ggTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggac                                                                                             >  2:804041/1‑100 (MQ=255)
ggTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggac                                                                                             >  1:799728/1‑100 (MQ=255)
 gTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacg                                                                                            >  2:866066/1‑100 (MQ=255)
  tAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgc                                                                                           <  2:1160122/100‑1 (MQ=255)
    aGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcgg                                                                                         >  2:254550/1‑100 (MQ=255)
     gCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcgga                                                                                        >  1:1386482/1‑100 (MQ=255)
     gCCTAAAAAACAGGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcgga                                                                                        >  2:1120792/1‑100 (MQ=255)
      ccTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcggac                                                                                       <  2:656496/100‑1 (MQ=255)
       cTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcggacg                                                                                      >  1:679383/1‑100 (MQ=255)
       cTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcggacg                                                                                      >  2:734154/1‑100 (MQ=255)
        tAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcggacgc                                                                                     >  2:1061480/1‑100 (MQ=255)
        tAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGggacgcggacgc                                                                                     >  2:1209583/1‑100 (MQ=255)
         aaaaaaCATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGcc                                                                                    >  1:891484/1‑100 (MQ=255)
            aaaCATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTc                                                                                 >  2:1409325/1‑100 (MQ=255)
                 tGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTgg                                                                            <  2:133031/100‑1 (MQ=255)
                  gTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGt                                                                           >  2:1380127/1‑100 (MQ=255)
                    tCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGa                                                                         >  1:247530/1‑100 (MQ=255)
                    tCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGa                                                                         <  2:1686039/100‑1 (MQ=255)
                      aTGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGAtt                                                                       >  2:1183029/1‑100 (MQ=255)
                       tGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGAttt                                                                      <  1:722671/100‑1 (MQ=255)
                       tGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGAttt                                                                      >  1:1320786/1‑100 (MQ=255)
                        gCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTc                                                                     >  1:904063/1‑100 (MQ=255)
                         ccTGGCCCGGCCTCCGTTTAAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCt                                                                    <  2:76919/100‑1 (MQ=255)
                          cTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTa                                                                   >  2:1494802/1‑100 (MQ=255)
                            ggCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAAt                                                                 >  1:618514/1‑100 (MQ=255)
                                cGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACgg                                                             <  2:627660/100‑1 (MQ=255)
                                 ggCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGt                                                            >  1:693287/1‑100 (MQ=255)
                                    cTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGc                                                         >  2:684860/1‑100 (MQ=255)
                                     tCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGca                                                        >  1:383775/1‑100 (MQ=255)
                                     tCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGca                                                        <  2:1009312/100‑1 (MQ=255)
                                     tCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGca                                                        >  1:582773/1‑100 (MQ=255)
                                        gTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGcatcc                                                     >  1:788445/1‑100 (MQ=255)
                                         tttCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGcatcca                                                    >  2:1177396/1‑100 (MQ=255)
                                              gAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATcgcg                                               >  1:111530/1‑100 (MQ=255)
                                              gAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATcgcg                                               >  1:595638/1‑100 (MQ=255)
                                              gAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATcgcg                                               >  1:797944/1‑100 (MQ=255)
                                                acacCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGtt                                             <  1:1527000/100‑1 (MQ=255)
                                                 cacCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGttt                                            <  2:891486/100‑1 (MQ=255)
                                                   ccAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTg                                          <  1:1053867/100‑1 (MQ=255)
                                                    cAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGa                                         >  1:785728/1‑100 (MQ=255)
                                                     aGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGAc                                        <  1:254550/100‑1 (MQ=255)
                                                     aGGTGTTCTGCGAGCACAATCAAAGACATACCCGCATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGAc                                        <  2:490336/100‑1 (MQ=255)
                                                      ggTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGAcc                                       >  1:857945/1‑100 (MQ=255)
                                                      ggTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGAcc                                       <  2:494441/100‑1 (MQ=255)
                                                       gtgtTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGAccc                                      <  2:247530/100‑1 (MQ=255)
                                                           tCTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCtt                                  >  2:1555674/1‑100 (MQ=255)
                                                            cTGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCttt                                 <  2:1456214/100‑1 (MQ=255)
                                                             tGCGCGCACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTc                                >  2:1524421/1‑100 (MQ=255)
                                                              gcgcgcACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCa                               <  1:643849/100‑1 (MQ=255)
                                                              gcgcgcACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCa                               <  2:460172/100‑1 (MQ=255)
                                                              gcgcgcACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCa                               <  1:684858/100‑1 (MQ=255)
                                                               cgcgcACAATCAAAGACATGTCCGAAGTCAGTGGGACGCGGACGCCGGCTTTGGTGAGTTGTATTACGGTGGCATCCATCGCGTTTTGACCCGCtt                                  <  1:616377/96‑1 (MQ=255)
                                                                  gcACAATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCtt                           >  2:1681030/1‑100 (MQ=255)
                                                                   cacaATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTGATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCtta                          <  1:932709/100‑2 (MQ=255)
                                                                   cacaATCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTc                          <  1:616312/100‑1 (MQ=255)
                                                                       aTCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAggg                      >  2:542455/1‑100 (MQ=255)
                                                                       aTCAAAGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAggg                      >  2:423949/1‑100 (MQ=255)
                                                                         cAAGGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCt                    <  1:244471/100‑1 (MQ=255)
                                                                          aaaGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCtt                   >  2:961845/1‑100 (MQ=255)
                                                                          aaaGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCtt                   >  2:316423/1‑100 (MQ=255)
                                                                           aaGACATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTTATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTg                  <  2:1688816/100‑1 (MQ=255)
                                                                               cttACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCg              <  1:1611499/98‑1 (MQ=255)
                                                                               cATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCg              >  2:1437941/1‑100 (MQ=255)
                                                                               cATACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCg              <  1:1177390/100‑1 (MQ=255)
                                                                                 tACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGAc            <  2:239894/100‑1 (MQ=255)
                                                                                 tACCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGag            >  2:1611506/1‑99 (MQ=255)
                                                                                  aCCCGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACa           >  2:1588213/1‑100 (MQ=255)
                                                                                   cccGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGCCAg          >  1:771848/1‑100 (MQ=255)
                                                                                   cccGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAg          >  1:75825/1‑100 (MQ=255)
                                                                                   cccGAATTCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTGTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAg          >  2:849208/1‑100 (MQ=255)
                                                                                         ttCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCagag    <  2:189258/100‑1 (MQ=255)
                                                                                         ttCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCagag    <  2:797946/100‑1 (MQ=255)
                                                                                          tCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGc   <  2:1386489/100‑1 (MQ=255)
                                                                                          tCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGc   >  2:956430/1‑100 (MQ=255)
                                                                                          tCAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGc   <  1:866064/100‑1 (MQ=255)
                                                                                           cAGCGGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCt  >  2:1853/1‑100 (MQ=255)
                                                                                               |                                                                                               
GGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCT  >  NC_000913/1914793‑1914983

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