Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,367,695 C→T 9.7% R118H (CGT→CAT)  PP_3844 ← D‑aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,367,6950CT 9.7% 108.7 / 2.8 41R118H (CGT→CAT) PP_3844D‑aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (20/17);  new base T (2/2);  total (22/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 = 7.72e-01

GGCCCAATTACAGCATTACTCAGGGCTTCGGCCACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATGGACGCCGGCAAAGC  >  NC_002947/4367551‑4367842
                                                                                                                                                |                                                                                                                                                   
ggCCCAATTACAGCATTACTCAGGGCTTCGGCCACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTc                                                                                                                                                 <  3:197916/149‑1 (MQ=255)
    cAATTACAGCATTACTCAGGGCTTCGGCCACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCATGGTCCAg                                                                                                                                              >  5:223288/1‑148 (MQ=255)
    cAATTACAGCATTACTCAGGGCTTCGGCCACATGCTCGTGCTTCACATGAAAACCGGTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGGCATCCGGTTGTTGTTCACGGTCCAGc                                                                                                                                             <  5:267085/149‑1 (MQ=255)
         aCAGCATTACTCAGGCCTTCGGCCACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAAt                                                                                                                                        >  5:266126/1‑149 (MQ=255)
                  cTCAGGGCTTCGGCCACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATc                                                                                                                               <  8:89953/149‑1 (MQ=255)
                     aGGGCTTCGGCCTCTTGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCa                                                                                                                              <  7:206288/147‑1 (MQ=255)
                       ggCTTCGGCCACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGca                                                                                                                          >  7:151926/1‑149 (MQ=255)
                               ccACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGca                                                                                                                          >  5:7428/1‑141 (MQ=255)
                               ccACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGca                                                                                                                          <  6:7428/141‑1 (MQ=255)
                                 acaTGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAg                                                                                                                >  3:315512/1‑149 (MQ=255)
                                             ttCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGtt                                                                                                    >  1:132990/1‑149 (MQ=255)
                                             ttCACATGAAAACCGTTGATGTCGATGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGtt                                                                                                    <  6:266126/149‑1 (MQ=255)
                                                        aCCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCCCCGGCATGTTCCAGTAGATACCGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGa                                                                                         <  6:109771/149‑1 (MQ=255)
                                                         ccGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTAATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGAt                                                                                        <  6:198395/149‑1 (MQ=255)
                                                          cGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATg                                                                                       >  8:110944/1‑149 (MQ=255)
                                                                        agcagcCCGTCATAGGTCTCCAGCACCACCGGCATGTTCTAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACCCACAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCaa                                                                            >  7:87261/1‑146 (MQ=255)
                                                                        agcagcCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCaaaa                                                                          <  1:238819/148‑1 (MQ=255)
                                                                         gcagcCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGa                                                                        >  6:292333/1‑149 (MQ=255)
                                                                          cagcCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCATGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGAc                                                                       <  6:223288/149‑1 (MQ=255)
                                                                                       gTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTcgc                                                              >  1:262315/1‑145 (MQ=255)
                                                                                       gTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTcgc                                                              <  2:262315/145‑1 (MQ=255)
                                                                                        tctcCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATc                                                         >  1:239222/1‑149 (MQ=255)
                                                                                            cAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCAc                                                      <  2:239222/148‑1 (MQ=255)
                                                                                            cAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACt                                                     >  8:149822/1‑149 (MQ=255)
                                                                                             aGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCCGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTc                                                    >  5:249636/1‑149 (MQ=255)
                                                                                                  caccGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGAc                                               <  7:110944/149‑1 (MQ=255)
                                                                                                  caccGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGAc                                               >  5:270436/1‑149 (MQ=255)
                                                                                                   accGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCATGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGAcc                                              >  6:184685/1‑149 (MQ=255)
                                                                                                               cAGTAGATACGGCCGTCATCCGGTTGTTGTTCATGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGc                                   <  5:184685/148‑1 (MQ=255)
                                                                                                                        cggccggcAGCCGGTTGTTGGTCACGGTCCAGCACAATCCGTGCATCACGCACCACACCCAGGCTGGGGGTGTTGGGGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTGGCGTCGCCGTTGCCAt                          <  8:151926/148‑1 (MQ=255)
                                                                                                                            cGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCaa                      >  4:309997/1‑148 (MQ=255)
                                                                                                                             gTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAaca                    <  2:132990/149‑1 (MQ=255)
                                                                                                                                aTCCGGTTGTTGTTCACGGTCCATCACAATCAGTGCATCACGCACCACACCCCGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACAt                 <  4:315512/149‑1 (MQ=255)
                                                                                                                                aTCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACAt                 >  6:244571/1‑149 (MQ=255)
                                                                                                                                  ccGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATg                >  7:65305/1‑148 (MQ=255)
                                                                                                                                    ggTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATGGAc             >  1:177225/1‑149 (MQ=255)
                                                                                                                                         ttgttCACGGTCCAGCACAATCAGTGGATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATGGACGCCgg        >  1:182874/1‑149 (MQ=255)
                                                                                                                                         ttgttCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATGGACGCCg         <  7:149822/148‑1 (MQ=255)
                                                                                                                                            ttCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATGGACGCCGGCa      >  1:294579/1‑148 (MQ=255)
                                                                                                                                              cACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATGGACGCCGGCAAAg   >  2:125495/1‑149 (MQ=255)
                                                                                                                                               aCGGTCCAGCACAATCAGTGCATCACGCAACACACCCAGGCTGTGGGTGTTGTTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATGGACGCCGGCAAAGc  <  8:87261/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGCCCAATTACAGCATTACTCAGGGCTTCGGCCACATGCTCGTGCTTCACATGAAAACCGTTGATGTCGTTGAGCAGCCCGTCATAGGTCTCCAGCACCACCGGCATGTTCCAGTAGATACGGCCGTCATCCGGTTGTTGTTCACGGTCCAGCACAATCAGTGCATCACGCACCACACCCAGGCTGTGGGTGTTGGTGAAGGCGATGGGACTGGTCAAAAGACCCGCTTCGCGAATCCACTCCAGACCGGTTGCGTCGCCGTTGCCATTCAACACATGGACGCCGGCAAAGC  >  NC_002947/4367551‑4367842

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