Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 867,319 A→G Y251C (TAC→TGC)  iaaA → Isoaspartyl peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913867,3190AG90.5% 60.5 / ‑1.1 21Y251C (TAC→TGC) iaaAIsoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (13/6);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

TGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCG  >  NC_000913/867191‑867444
                                                                                                                                |                                                                                                                             
tGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTa                                                                                                                     <  2:773528/139‑1 (MQ=255)
            cTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGaa                                                                                                         >  1:194668/1‑139 (MQ=255)
                       ccGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCgg                                                                                                           <  2:435629/126‑1 (MQ=255)
                       ccGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCgg                                                                                                           >  1:435629/1‑126 (MQ=255)
                          gATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGgt                                                                                           >  2:681666/1‑139 (MQ=255)
                                      aTAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGaaaa                                                                               >  1:328061/1‑139 (MQ=255)
                                                tgtgGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTgcg                                                                     >  2:123914/1‑139 (MQ=255)
                                                  tgGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCt                                                                   >  1:272671/1‑139 (MQ=255)
                                                                    gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg                                                                                            >  1:331023/1‑96 (MQ=255)
                                                                    gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg                                                                                            <  2:331023/96‑1 (MQ=255)
                                                                    gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCtt                                                 >  1:125609/1‑139 (MQ=255)
                                                                    gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCtt                                                 >  2:237302/1‑139 (MQ=255)
                                                                                tCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGa                                     <  2:125609/139‑1 (MQ=255)
                                                                                    cATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg                                                                                            >  2:107554/1‑80 (MQ=255)
                                                                                    cATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg                                                                                            <  1:107554/80‑1 (MQ=255)
                                                                                                               cgccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAac      >  2:655072/1‑139 (MQ=255)
                                                                                                               cgccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAac      >  2:696202/1‑139 (MQ=255)
                                                                                                                gccgcGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCt                                          >  2:162244/1‑102 (MQ=255)
                                                                                                                gccgcGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCt                                          <  1:162244/102‑1 (MQ=255)
                                                                                                                 ccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAacac    >  1:438808/1‑139 (MQ=255)
                                                                                                                   gcgTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCg  <  1:237302/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                             
TGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCG  >  NC_000913/867191‑867444

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: