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,3190AG97.4% 132.9 / ‑4.1 38Y251C (TAC→TGC) iaaAIsoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (18/19);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

GACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACAC  >  NC_000913/867184‑867442
                                                                                                                                       |                                                                                                                           
gACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGcgg                                                                                                                          >  2:62525/1‑139 (MQ=255)
                tCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTcgcg                                                                                                          >  1:968329/1‑139 (MQ=255)
                  ccTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGa                                                                                                        <  2:968329/139‑1 (MQ=255)
                        tGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTg                                                                                                  <  2:1070944/139‑1 (MQ=255)
                          ggTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCg                                                                                                >  2:368243/1‑139 (MQ=255)
                          ggTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCg                                                                                                >  1:816276/1‑139 (MQ=255)
                            tGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAg                                                                                              >  2:865253/1‑139 (MQ=255)
                             gCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGc                                                                                             <  2:514642/139‑1 (MQ=255)
                              ccGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGGCTGCGAGCg                                                                                            >  1:653468/1‑139 (MQ=255)
                                ggATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTTTACCGGCAGGGGCGAAGGCTTCATCCGCGCGGTGGCGGCATATGACATCGTCGAGGTAATGGATTGCGGCTGATTAAGTCTCGCGGAAGCCTGCGAGCggg                                                                                          <  1:824694/139‑1 (MQ=255)
                                ggATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg                                                                                          <  2:832368/139‑1 (MQ=255)
                                      taggCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTaa                                                                                    <  1:738738/136‑1 (MQ=255)
                                             aTAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGaaaa                                                                             >  2:505437/1‑139 (MQ=255)
                                              tAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGaaaaa                                                                            >  1:518434/1‑139 (MQ=255)
                                                       tgtgGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTgcg                                                                   >  1:844279/1‑139 (MQ=255)
                                                          ggCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCtt                                                                <  1:62525/139‑1 (MQ=255)
                                                                 acTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTa                                                         >  1:692797/2‑139 (MQ=255)
                                                                           gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCtt                                                                >  1:236150/1‑122 (MQ=255)
                                                                           gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCtt                                                                <  2:236150/122‑1 (MQ=255)
                                                                           gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCtt                                               >  1:523582/1‑139 (MQ=255)
                                                                           gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCtt                                               >  1:496373/1‑139 (MQ=255)
                                                                             aCGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTaa                                             <  2:816276/139‑1 (MQ=255)
                                                                              cGGGCGAAGTCTTCATCCGCGTGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAAt                                            <  2:653468/139‑1 (MQ=255)
                                                                                 gCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGc                                         <  2:1101002/139‑1 (MQ=255)
                                                                                    aaGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTAt                                      <  2:84519/139‑1 (MQ=255)
                                                                                      gTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCg                                    >  2:214852/1‑139 (MQ=255)
                                                                                           cATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCAt                               <  2:633981/139‑1 (MQ=255)
                                                                                                   cgcTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGaa                       <  2:692797/139‑1 (MQ=255)
                                                                                                   cgcTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGaa                       >  1:150148/1‑139 (MQ=255)
                                                                                                        gcggcATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTcg                  <  1:865253/139‑1 (MQ=255)
                                                                                                             atatGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTa             <  1:368243/139‑1 (MQ=255)
                                                                                                                      cgccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGc                                                                                             >  1:857020/1‑50 (MQ=255)
                                                                                                                      cgccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGc                                                                                             <  2:857020/50‑1 (MQ=255)
                                                                                                                      cgccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAac    >  1:1034500/1‑139 (MQ=255)
                                                                                                                        ccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg                                                                                          <  2:731799/51‑1 (MQ=255)
                                                                                                                        ccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg                                                                                          >  1:731799/1‑51 (MQ=255)
                                                                                                                        ccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAacac  <  1:223340/139‑1 (MQ=255)
                                                                                                                        ccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAacac  <  2:150148/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                           
GACGAGTTGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACAC  >  NC_000913/867184‑867442

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