Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,054,245 G→A A199V (GCG→GTG)  pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,054,2450GA97.2% 95.1 / NA 36A199V (GCG→GTG) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/18);  minor base C (0/1);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGATCATCGTTGCCGGTGCGGTGAGATTT  >  NC_000913/3054109‑3054376
                                                                                                                                        |                                                                                                                                   
gTGCAGAATGTAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCAccc                                                                                                                                  <  1:906930/140‑1 (MQ=255)
 tGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGGACTCGAACATTCCCCGGCGGCATTTTTCCATCAcccc                                                                                                                                 >  2:941884/1‑139 (MQ=255)
   cAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCgtg                                                                                                                               <  1:133269/140‑1 (MQ=255)
      aaTGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTTGTGGGGAATTTCGCCTTCCAGATGGGACTCGAACATTCCCGGACGGCATTTTTTCATCACCCgggctt                                                                                                                            >  2:632294/1‑135 (MQ=255)
      aaTGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTAAGATTCGGCCTCCCGGTGGTGCTCGAACAATCCCGGACCGCATTTTTCCATCACCCCGGTAt                                                                                                                            >  2:588235/1‑140 (MQ=255)
        tGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAAATCGTGGTGAATTTCGCCTTCCAGATGGTACTCGCACATTTCCCGACGGGATTTTTTCCTCACCCGTCTAttt                                                                                                                          >  2:831739/1‑138 (MQ=255)
            gCCGTTTTCACCGCTGCCGACACTGCTGTTATAGGACAGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCa                                                                                                                      <  1:371469/140‑1 (MQ=255)
              cGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATg                                                                                                                    <  1:998062/140‑1 (MQ=255)
                 tttCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCg                                                                                                                 <  1:878671/140‑1 (MQ=255)
                  ttCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGGGAATTTCGCCTTCCAGATGGGACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCgg                                                                                                                >  2:31152/1‑140 (MQ=255)
                                   tGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGCTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCCGACGGCATTTTTTCATTAACCGCGTATTTGCCATGGCGGCGGTTTCCTCcgccccg                                                                                               >  2:132175/1‑136 (MQ=255)
                                     gtgtTATAGGACGGATAGCGCGCACCGTGGCGGTTAATTTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCgcggcg                                                                                             <  1:415557/140‑1 (MQ=255)
                                       gtTATAGGACGGATAGCGCGCAGCGGGGCCGTTAAATTCGGGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCA‑TTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAg                                                                                          <  1:366686/140‑1 (MQ=255)
                                        ttATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAg                                                                                          <  1:831739/140‑1 (MQ=255)
                                             ggACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGGGATTTTTCCCCCGCGGCGGGGGACgg                                                                                     >  1:673467/1‑140 (MQ=255)
                                                cGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGGGATTTCTCCCCCGCGGCGGGGGACGGCaa                                                                                  >  1:224728/1‑140 (MQ=255)
                                                cGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCCTGTATGTGCCATGGCGGGGGTTTTTCCCGCCCGGCGGGGGACGGCaa                                                                                  >  1:547679/1‑140 (MQ=255)
                                                   aTAGCGCGCACCGTGGCGGTTAAAATCGTGGTGAATTTTGCCTTTCAGATGGGACTCGAAAAATTCCGGGCGGCCTTTTTTCATCACCCGCGTATGTGGCATGGCGGGGATTTTTCCCCCTCGGGGGGGTAACGCAAact                                                                               >  2:801970/1‑140 (MQ=255)
                                                         cgcACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGGACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGGGGTTTTTCCCCCGCGGCGGGGGACCGCAAAAtcttctt                                                                         >  2:630720/1‑139 (MQ=255)
                                                            aCCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGGTTTTTCCCCCGCGGCGGGGGACGGCAATCTCTTCTTGCg                                                                      >  1:10259/1‑140 (MQ=255)
                                                                 ggCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGGGATTTTTCCCCCGCGGCGGGGGACGGCAAACTCTTTTGGCGATTTg                                                                 >  1:108285/1‑140 (MQ=255)
                                                                   cGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGaa                                                               <  1:132175/140‑1 (MQ=255)
                                                                         aaTTCGTGGTGAATTTCCCCTTACAGATGGTTCTCTAACCTTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGgcg                                                         <  2:299102/140‑1 (MQ=255)
                                                                             cGTGGTGAATTTCGCCTTCCAGATGGTGCTCGAACATTCCCGGACGGCGTTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCAtt                                                     <  1:53293/140‑1 (MQ=255)
                                                                               tggtgAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTGCGGCAATTTCTACTGGCGATTGGAACAGGCGCATTTc                                                   >  1:325429/1‑140 (MQ=255)
                                                                               tggtgAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTc                                                   <  2:601190/140‑1 (MQ=255)
                                                                                ggtgAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCa                                                  >  2:871291/1‑140 (MQ=255)
                                                                                 gtgAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCAt                                                 <  1:630417/140‑1 (MQ=255)
                                                                                               ccAGATGGTACTCGAGCATTCCCGGCCGGAATTTTTCTATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGc                                   <  1:401282/140‑1 (MQ=255)
                                                                                                 aGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCa                                 <  1:239941/140‑1 (MQ=255)
                                                                                                 aGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCa                                 <  1:801970/140‑1 (MQ=255)
                                                                                                        aCTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGatc                          >  2:100529/1‑140 (MQ=255)
                                                                                                         cTCGAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGatca                         >  1:624332/1‑140 (MQ=255)
                                                                                                            gAACATTCCCGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTTATGGACAACAGGATGCCAGGCGATCATCg                      >  1:437675/1‑140 (MQ=255)
                                                                                                                 ttCCCGGATGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGatcatc                       >  2:18434/1‑134 (MQ=255)
                                                                                                                 ttCCCGGATGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGatcatc                       <  1:18434/134‑1 (MQ=255)
                                                                                                                    ccGGACGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAAGAGGCGCATCTCAAGAAGAACAGGACGCCAGCCGATCATCGTTGCcggt              >  2:286824/1‑140 (MQ=255)
                                                                                                                    ccGGACGGCATTTTTCCATCACCCGTGTATGGGCCCTGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGATCATCGTTGCcggt              <  2:370926/140‑1 (MQ=255)
                                                                                                                        aCGGCATTTTTCCATCACCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGATCATCGTTGCcggtgcgg          >  2:246905/1‑140 (MQ=255)
                                                                                                                            cATTTTTCCATCACCCGTGTATGTGCCATGGCGGCGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGGCGATCATCGGTGCCGGTGCGGTgag      >  1:419102/1‑140 (MQ=255)
                                                                                                                                tttCCATCCCCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGATCATCGTTGCCGGTGCGGTGAGAttt  <  1:961955/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
GTGCAGAATGCAGCCGTTTTCACCGCTGCCGACAATGGTGTTATAGGACGGATAGCGCGCACCGTGGCGGTTAAATTCGTGGTGAATTTCGCCTTCCAGATGGTACTCGAACATTCCCGGACGGCATTTTTCCATCGCCCGTGTATGTGCCATGGCGGTGATTTCTCCCGCGCGGCGGAGTACGGCAATCTCTTCTGGCGATTTGAACAGGCGCATTTCATGAACAACAGGACGCCAGTCGATCATCGTTGCCGGTGCGGTGAGATTT  >  NC_000913/3054109‑3054376

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