Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,958,119 G→A 11.3% P256S (CCG→TCG)  ptrA ← protease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,958,1190GA11.3% 146.9 / 9.3 53P256S (CCG→TCG) ptrAprotease III
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base A (3/3);  total (23/30)
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.58e-01

CGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCT  >  NC_000913/2957982‑2958266
                                                                                                                                         |                                                                                                                                                   
cGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCaa                                                                                                                                          >  2:492126/1‑149 (MQ=255)
  cTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAg                                                                                                                                        <  2:232896/149‑1 (MQ=255)
        tgttATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTcc                                                                                                                                  <  3:561261/149‑1 (MQ=255)
            aTCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTccgccg                                                                                                                              <  2:276881/149‑1 (MQ=255)
             tCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCcgccgc                                                                                                                             >  4:321829/1‑149 (MQ=255)
                 tGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGACACGCGACCAAAGGTGTCCGCCGCCAtt                                                                                                                         >  2:167444/1‑149 (MQ=255)
                 tGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGACACGCGACCAAAGGTGTCCGCCGCCAtt                                                                                                                         >  2:508342/1‑149 (MQ=255)
                  gCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCAttt                                                                                                                        <  2:480038/149‑1 (MQ=255)
                     aaaCTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCaaa                                                                                                                                         <  2:450826/129‑1 (MQ=255)
                     aaaCTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCaaa                                                                                                                                         >  1:450826/1‑129 (MQ=255)
                                 cAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCa                                                                                                         >  2:526359/1‑149 (MQ=255)
                                 cAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCa                                                                                                         >  2:447238/1‑149 (MQ=255)
                                      cTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGt                                                                                                    >  2:342295/1‑149 (MQ=255)
                                       tttACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGtt                                                                                                   <  3:235155/149‑1 (MQ=255)
                                         tACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTtta                                                                                                 <  2:73206/149‑1 (MQ=255)
                                            gcgGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTttatta                                                                                              >  3:438950/1‑149 (MQ=255)
                                                 aGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTa                                                                                         >  1:163932/1‑149 (MQ=255)
                                                        ggACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAAcc                                                                                  <  2:378255/149‑1 (MQ=255)
                                                             tAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCtt                                                                             >  2:511989/1‑149 (MQ=255)
                                                              aatgaatgATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTtc                                                                            <  4:62246/149‑1 (MQ=255)
                                                                tgaatgaTAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACt                                                                                            <  2:145499/131‑1 (MQ=255)
                                                                tgaatgaTAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACt                                                                                            >  1:145499/1‑131 (MQ=255)
                                                                tgaatgaTAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTtcat                                                                          >  2:317299/1‑149 (MQ=255)
                                                                 gaatgaTAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTtcatc                                                                         >  4:216907/1‑149 (MQ=255)
                                                                 gaatgaTAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTtcatc                                                                         >  2:493728/1‑149 (MQ=255)
                                                                    tgaTAATGTCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCaaa                                                                      <  1:492126/149‑1 (MQ=255)
                                                                       tAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTg                                                                   <  1:66890/149‑1 (MQ=255)
                                                                         aTGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGc                                                                 <  2:70594/149‑1 (MQ=255)
                                                                         aTGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGc                                                                 <  1:493728/149‑1 (MQ=255)
                                                                         aTGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGc                                                                 <  1:447238/149‑1 (MQ=255)
                                                                         aTGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGACACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGc                                                                 <  1:508342/149‑1 (MQ=255)
                                                                         aTGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGACACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGc                                                                 <  1:167444/149‑1 (MQ=255)
                                                                             ccTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGaa                                                             <  4:421765/149‑1 (MQ=255)
                                                                             ccTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCCGaa                                                             >  4:443253/1‑149 (MQ=255)
                                                                             ccTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTACGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAGATAACCGCCTTCATCAAATTGGCGGaa                                                             <  1:548223/149‑1 (MQ=255)
                                                                               ttttGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATa                                                           <  1:526359/149‑1 (MQ=255)
                                                                                        tCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTtctc                                                  <  2:331495/149‑1 (MQ=255)
                                                                                          ggTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAACGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGt                                                <  2:416372/149‑1 (MQ=255)
                                                                                          ggTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGt                                                <  2:547151/149‑1 (MQ=255)
                                                                                           gTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTg                                               >  1:526750/1‑149 (MQ=255)
                                                                                           gTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTg                                               >  1:565972/1‑149 (MQ=255)
                                                                                           gTGACTACCGGCACGGTGATTTCCGGTTTTTTACTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTg                                               >  3:105504/1‑149 (MQ=255)
                                                                                            tGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTtc                                                                            <  2:66946/119‑1 (MQ=255)
                                                                                            tGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTtc                                                                            >  1:66946/1‑119 (MQ=255)
                                                                                               cTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCGTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGaaa                                           <  4:246635/149‑1 (MQ=255)
                                                                                                       aCGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAg                                   >  4:356354/1‑149 (MQ=255)
                                                                                                         ggTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAgcg                                 >  2:326574/1‑149 (MQ=255)
                                                                                                            gATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCct                              <  2:163932/149‑1 (MQ=255)
                                                                                                                    gTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCAcc                      <  4:489687/149‑1 (MQ=255)
                                                                                                                            ctctctTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTAcc              <  4:535159/149‑1 (MQ=255)
                                                                                                                               tctTTGTTCGACACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGAt                  >  3:160046/1‑142 (MQ=255)
                                                                                                                               tctTTGTTCGACACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGAt                  <  4:160046/142‑1 (MQ=255)
                                                                                                                                ctTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGt          <  1:550762/149‑1 (MQ=255)
                                                                                                                                        ggCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCt  >  3:30100/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
CGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGATAATGCCCTTTTGCGCGTCGGTGACTACCGGCACGGTGATTTCCGGTTTTTTGCTCTCTTTGTTCGGCACGCGACCAAAGGTGTCCGCCGCCATTTTTGCCAACTCCGGCAGCGGTTTATTACTGTAAATAACCGCCTTCATCAAATTGGCGGAATAGTACTTCTCGTGGAAATCTTTCAGCGCCTGCTGCACCGGATTACCAGGTTTGTCGCT  >  NC_000913/2957982‑2958266

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑