Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,028,519 C→T 8.3% P263S (CCG→TCG)  PP_4432 → Xaa‑Pro aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,028,5190CT 8.3% 137.5 / 2.5 48P263S (CCG→TCG) PP_4432Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base T (2/2);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GGCGGGGACTACCCGGCCATCGTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCT  >  NC_002947/5028375‑5028667
                                                                                                                                                |                                                                                                                                                    
ggCGGGGACTACCCGGCCATCGTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCt                                                                                                                                                  >  3:1062/1‑149 (MQ=255)
   ggggACTACCCGGCCATCGTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACGTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTc                                                                                                                                               <  2:213461/149‑1 (MQ=255)
    gggACTACCCGGCCATCGTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTcc                                                                                                                                              <  5:165486/149‑1 (MQ=255)
    gggACTACCCGGCCATCGTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTcc                                                                                                                                              <  2:229750/149‑1 (MQ=255)
                ccATCGTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTAc                                                                                                                                  <  4:1062/149‑1 (MQ=255)
                ccATCGTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTAc                                                                                                                                  <  2:203681/149‑1 (MQ=255)
                     gTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTggg                                                                                                                             <  1:155888/149‑1 (MQ=255)
                                tCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCg                                                                                                                   >  5:91679/1‑148 (MQ=255)
                                       ggAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCAtt                                                                                                            >  6:261865/1‑148 (MQ=255)
                                           cAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATc                                                                                                       <  2:287060/149‑1 (MQ=255)
                                                gCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTccc                                                                                                                                             >  4:65657/1‑106 (MQ=255)
                                                  tAGCGCCCCACACCTGACCTGGGATGACTCTCCAATCAAGGTCGGCGCCGGCACCTACTTCGAAATCGCCGGCTGCTACAAGCGCAATCACAGCCCGCTGTCCCGAACTGTTTTCCTGGGCAAGCCACCGCAGCATTTCATCGACGGc                                                                                                 >  6:233650/1‑148 (MQ=255)
                                                         ccACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGAc                                                                                                    <  6:107069/138‑1 (MQ=255)
                                                         ccACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGAc                                                                                                    >  5:107069/1‑138 (MQ=255)
                                                           acacCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAg                                                                                       <  5:261865/149‑1 (MQ=255)
                                                               ctgacctgGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCTCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCa                                                                                        >  2:32294/1‑144 (MQ=255)
                                                               ctgacctgGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCTCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCa                                                                                        >  2:341219/1‑144 (MQ=255)
                                                               ctgacctgGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCTCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCa                                                                                        <  1:32294/144‑1 (MQ=255)
                                                               ctgacctgGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCTCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCa                                                                                        <  1:341219/144‑1 (MQ=255)
                                                               ctgacctgGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAgtcgt                                                                                   >  4:311880/1‑149 (MQ=255)
                                                                              tctcCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGc                                                                    >  8:260440/1‑149 (MQ=255)
                                                                                tcCATTCAAGGTCGGCGCCGGCACCTTATTCGAAATAGACGGCTCCTACAAGCGCTATCACTGCCCGCTGTACCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCCCTGCCg                                                                  >  2:240306/1‑149 (MQ=255)
                                                                                 ccATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCgg                                                                 >  3:95073/1‑149 (MQ=255)
                                                                                    ttCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCaa                                                                                                                          <  2:175380/89‑1 (MQ=255)
                                                                                    ttCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCaa                                                                                                                          >  1:175380/1‑89 (MQ=255)
                                                                                       aaGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGaaa                                                                                           >  8:7349/1‑117 (MQ=255)
                                                                                       aaGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGaaa                                                                                           <  7:7349/117‑1 (MQ=255)
                                                                                        aGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAg                                                          >  1:177513/1‑149 (MQ=255)
                                                                                          gTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGgc                                                        >  1:8603/1‑149 (MQ=255)
                                                                                                      aCCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTgg                                            <  4:95073/149‑1 (MQ=255)
                                                                                                       ccTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGc                                           <  3:311880/149‑1 (MQ=255)
                                                                                                        cttcttCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCa                                          <  2:177513/149‑1 (MQ=255)
                                                                                                             tcGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAacac                                      >  5:137222/1‑148 (MQ=255)
                                                                                                                aaaTCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGGGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACt                                  >  4:95843/1‑149 (MQ=255)
                                                                                                                   tCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTTCCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGt                               <  5:233650/149‑1 (MQ=255)
                                                                                                                   gtgCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGt                               <  5:201198/147‑1 (MQ=255)
                                                                                                                         gctgctACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCga                                                                                               >  4:140879/1‑79 (MQ=255)
                                                                                                                         gctgctACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCga                                                                                               <  3:140879/79‑1 (MQ=255)
                                                                                                                              tACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCAc                     >  4:172001/1‑148 (MQ=255)
                                                                                                                                cAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGATAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGcc                                               <  2:61132/120‑1 (MQ=255)
                                                                                                                                cAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGATAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGcc                                               >  1:61132/1‑120 (MQ=255)
                                                                                                                                cAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGcc                                               >  7:307568/1‑120 (MQ=255)
                                                                                                                                cAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGcc                                               <  8:307568/120‑1 (MQ=255)
                                                                                                                                cAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCAc                     <  3:172001/146‑1 (MQ=255)
                                                                                                                                        aTCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACg                    >  7:130093/1‑139 (MQ=255)
                                                                                                                                        aTCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACg                    <  8:130093/139‑1 (MQ=255)
                                                                                                                                           aCTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGTGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGc       <  3:95843/149‑1 (MQ=255)
                                                                                                                                           aCTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGc       <  7:260440/149‑1 (MQ=255)
                                                                                                                                               cccGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGc   >  2:202299/1‑149 (MQ=255)
                                                                                                                                                ccGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCt  <  6:91679/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGCGGGGACTACCCGGCCATCGTGCCGCTGCTTCCTACCGGAGCAGACGCTAGCGCCCCACACCTGACCTGGGATGACTCTCCATTCAAGGTCGGCGCCGGCACCTTCTTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCT  >  NC_002947/5028375‑5028667

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