Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,585,681 G→T G381G (GGC→GGA ggt ← gamma‑glutamyltranspeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,6810GT97.7% 146.3 / NA 43G381G (GGC→GGAggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (30/12);  minor base A (1/0);  total (31/12)
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

CTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGAT  >  NC_000913/3585626‑3585741
                                                       |                                                            
cTTTATCCACCACAGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTATCCGGGTCGAAtt                                                  >  1:659175/1‑68 (MQ=255)
   tATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCg                                               >  2:211465/1‑68 (MQ=255)
      ccaccaCTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTg                                            >  2:1112104/1‑68 (MQ=255)
      ccaccaCTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTg                                            >  1:891411/1‑68 (MQ=255)
      ccaccaCTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTg                                            >  1:702026/1‑68 (MQ=255)
      ccaccaCTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTg                                            >  1:644023/1‑68 (MQ=255)
           acTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATgg                                       >  2:944193/1‑68 (MQ=255)
           acTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATgg                                       <  2:212817/68‑1 (MQ=255)
           acTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATgg                                       >  1:1011428/1‑68 (MQ=255)
           acTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTATCCGGGGCGAATTTCGCTGGATgg                                       >  2:815235/1‑68 (MQ=255)
           acTGAGTAATGGGAAGTTTGATTACTCTCATAAGGCGCGAGCTTTCCGGGGCGAATTTCGCTGGATgg                                       >  2:878839/1‑68 (MQ=255)
                gTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCg                                  >  2:816985/1‑68 (MQ=255)
                gTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCg                                  >  1:558717/1‑68 (MQ=255)
                gTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCg                                  >  1:904275/1‑68 (MQ=255)
                gTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCg                                  >  1:1078241/1‑68 (MQ=255)
                 tAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGc                                 <  2:644023/68‑1 (MQ=255)
                 tAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGc                                 <  1:900166/68‑1 (MQ=255)
                  aaTGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCt                                >  2:446220/1‑68 (MQ=255)
                    tGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCttt                              <  2:543044/68‑1 (MQ=255)
                       gTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTAtt                           >  1:104064/1‑68 (MQ=255)
                       gTAGTTTGATTACTCTCATAAGGCGCAAGCTTTCAGGGGCGAATTTCGCTGGATGGCTTCGCTTTAtt                           >  1:272030/1‑68 (MQ=255)
                       gTAGTTTGATTACTCTCATAAGGCGCAAGCTTAACGGGGCGAATTTAGCTGGATGGATTCGCTTTAtt                           >  1:282870/1‑68 (MQ=255)
                            ttGATTACTATCAAAAGGCGCAAGCTATCCGGGGCGAAATTCGCTGGATGGCTTCGCTTTATTGataa                      >  2:1177692/1‑67 (MQ=255)
                                  aCTCTCATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAAttt                >  2:1105044/1‑68 (MQ=255)
                                    tctcATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAAt                  >  2:1120328/1‑64 (MQ=255)
                                    tctcATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAAt                  <  1:1120330/64‑1 (MQ=255)
                                    tctcATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGa              >  2:1113837/1‑68 (MQ=255)
                                     ctcATAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGAt             >  1:170872/1‑68 (MQ=255)
                                         tAAGTCGCAAGCTTTCCGTGGCGAATTTCTCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGc         <  1:1177694/68‑1 (MQ=255)
                                         tAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGc         <  1:1113839/68‑1 (MQ=255)
                                         tAAGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGc         <  1:816985/68‑1 (MQ=255)
                                          aaGGCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCa        >  2:46195/1‑68 (MQ=255)
                                            ggCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAAt      >  1:991023/1‑68 (MQ=255)
                                            ggCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAAt      >  2:243707/1‑68 (MQ=255)
                                            ggCGCAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAAt      >  1:362903/1‑68 (MQ=255)
                                            ggCGCAAGCTTTCCGGGGAGAATTTCGCTGGATGGCTGCGCTTTAGTGATATCAATTTGATCGGCAAt      >  2:191243/1‑68 (MQ=255)
                                              cgcAGGCTTTCCGGGGCGAATGTCGCTGGATGGCTTCGCTTTATTGATATCAAGTTGATCGGCGATAg    >  2:644401/1‑68 (MQ=255)
                                              cgcAAGCTTTCCGGGGCGAATTTCGCTTGATGGCTTCGCTTTATAGATATCAATTTGATCGGCAATAg    <  1:1019465/68‑1 (MQ=255)
                                              cgcAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAg    >  2:689637/1‑68 (MQ=255)
                                              cgcAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAg    >  1:533882/1‑68 (MQ=255)
                                              cgcAAGCTTTCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAg    <  1:1120747/68‑1 (MQ=255)
                                                cAAGCTTTCCGGGGCGAATTTCGCCGGATGGCGTCGCTTTATTGATATCAATCTGATCGGCAATAGac  <  2:145365/68‑2 (MQ=255)
                                                cAAGCTTTCCGGGGCGAATTTCGCCGGACGGCTTCGCGTTATTGATATCAATTTGATCGGCAATAGac  <  2:146955/68‑2 (MQ=255)
                                                       |                                                            
CTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGAT  >  NC_000913/3585626‑3585741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: