Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,631,531 G→A G335G (GGC→GGT guaA ← GMP synthetase (glutamine aminotransferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,631,5310GA96.4% 87.8 / ‑4.1 28G335G (GGC→GGTguaAGMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/16);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTGCCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCGATGATTTTACGTTTTGCTTCCGGATCGTTTTCGCCAGCCAGCGCTGACAGGAAG  >  NC_000913/2631395‑2631657
                                                                                                                                        |                                                                                                                              
tGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTAcc                                                                                                                              <  1:358792/139‑1 (MQ=255)
       ctctTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCgcc                                                                                                                       >  1:297172/1‑139 (MQ=255)
            tCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCgccagcca                                                                                                                  >  1:21939/1‑139 (MQ=255)
                cGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTc                                                                                                              <  2:21939/139‑1 (MQ=255)
                 ggTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCa                                                                                                             >  1:80704/1‑139 (MQ=255)
                  gTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCAc                                                                                                            <  2:970571/139‑1 (MQ=255)
                   ttCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACg                                                                                                           <  2:744379/139‑1 (MQ=255)
                    tCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGt                                                                                                          <  1:664131/139‑1 (MQ=255)
                                 tCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAg                                                                                             <  2:297172/139‑1 (MQ=255)
                                      aTCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGctt                                                                                        <  1:154635/139‑1 (MQ=255)
                                        ctctTTCGTCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGcttct                                                                                      <  1:851910/139‑1 (MQ=255)
                                            ttCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTc                                                                                               >  1:349577/1‑126 (MQ=255)
                                            ttCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTc                                                                                               <  2:349577/126‑1 (MQ=255)
                                             tCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATc                                                                                 >  2:315919/1‑139 (MQ=255)
                                              cGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCg                                                                                <  2:12138/139‑1 (MQ=255)
                                                     ccgccCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACtt                                                                         <  2:80704/139‑1 (MQ=255)
                                                                   gtgAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGa                                                           >  2:817237/1‑139 (MQ=255)
                                                                           tgatgaCGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCgatgat                                                   >  2:1030471/1‑139 (MQ=255)
                                                                                  gtgtGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTGCCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCGATGATTTTACGt                                            >  1:1064659/1‑139 (MQ=255)
                                                                                       cTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCGATGATTTTACGTTTTGc                                       >  1:787768/1‑139 (MQ=255)
                                                                                           aCCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCGATGATTTTACGTTTTGCTTcc                                   >  2:744889/1‑139 (MQ=255)
                                                                                                 tGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGcttct                                                                                      <  2:784444/82‑1 (MQ=255)
                                                                                                 tGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTACCCTACGCCAGCCACTTCACGTCTTCCAGTTTCAGCGcttct                                                                                      >  1:784444/1‑82 (MQ=255)
                                                                                                           gcAGATTCGATAACGTCAGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCGATGATTTTACGTTTTGCTTCCGGATCGTTTTCgccag                   <  1:95260/139‑1 (MQ=255)
                                                                                                             cgATTCGATAACGTCAGGGTAGATTGTACCCTGCGCCAGCCACCTCACGTCTTCCAGTTTCAGTGCTTCTTCATCGAATACTTCAACGAGAACGCGACCGATTTTTTTACGTTTTGCTTCCGGATCGTTTTCgccagcc                 <  1:1058600/138‑1 (MQ=255)
                                                                                                                            aGGGTAGATGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCGATGATTTTACGTTTTGCTTCCGGATCGTTTTCGCCAGCCAGCGCTGACAGGAAg  <  1:817237/139‑1 (MQ=255)
                                                                                                                                   aTGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCg                                                        >  2:550101/1‑78 (MQ=255)
                                                                                                                                   aTGGTACCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCg                                                        <  1:550101/78‑1 (MQ=255)
                                                                                                                                        |                                                                                                                              
TGAACAGCTCTTTCAGCGGTTCAACCAGGCCCATCTTCATCTCTTTCGGCAGGCCGCCCACGTTGTGGTGAGATTTGATGACGTGTGCTTTACCGGTTGCAGACGCCGCAGATTCGATAACGTCAGGGTAGATGGTGCCCTGCGCCAGCCACTTCACGTCTTCCAGTTTCAGCGCTTCTTCATCGAATACTTCAACGAAAACGCGACCGATGATTTTACGTTTTGCTTCCGGATCGTTTTCGCCAGCCAGCGCTGACAGGAAG  >  NC_000913/2631395‑2631657

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