Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 491,776 T→C V122A (GTT→GCT)  apt → adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913491,7760TC96.4% 91.4 / ‑6.1 28V122A (GTT→GCT) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (12/15);  total (12/16)
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

TTGGCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCG  >  NC_000913/491650‑491908
                                                                                                                              |                                                                                                                                    
ttGGCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGACCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGga                                                                                                                          <  1:570113/139‑1 (MQ=255)
             ggTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCtgg                                                                                                                                 >  1:134273/1‑119 (MQ=255)
             ggTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCtgg                                                                                                                                 <  2:134273/119‑1 (MQ=255)
              gTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCGAACCGGGCGACAAAGCTCtggtggtg                                                                                                                            >  1:464675/1‑123 (MQ=255)
              gTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCtggtggtg                                                                                                                            <  2:464675/123‑1 (MQ=255)
                     aaCCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACcggcgg                                                                                                     <  2:474179/139‑1 (MQ=255)
                       ccGGGCAAACTGCCGCGTTAAACCATCAGTTAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCa                                                                                                   <  2:550058/139‑1 (MQ=255)
                           gCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTCGTGGACGACCTGCTGGCAACCGGCGGCACTAt                                                                                               >  2:670106/1‑139 (MQ=255)
                                 tGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTCGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGc                                                                                         <  1:670106/139‑1 (MQ=255)
                                      cgTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCg                                                                                    <  2:928717/139‑1 (MQ=255)
                                                            aCGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTggg                                                              <  1:517394/139‑1 (MQ=255)
                                                            aCGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTggg                                                              >  2:504205/1‑139 (MQ=255)
                                                                           gCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGa                                               <  2:890515/139‑1 (MQ=255)
                                                                            cACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATc                                                                         >  1:229205/1‑112 (MQ=255)
                                                                            cACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATc                                                                         <  2:229205/112‑1 (MQ=255)
                                                                              ccGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGc                                            >  1:880868/1‑139 (MQ=255)
                                                                                    aGCTGGAGATCCACGTTGATGCCATCGAACCGGGCGACAAAGTTCTGGTAGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGGCCGTTAAACTGATCCGTCGTCTGGGTGCTGAAGTGGCTGACGCTGCgct                                      <  2:558688/139‑3 (MQ=255)
                                                                                      cTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCa                                    >  1:660184/1‑139 (MQ=255)
                                                                                        ggAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAATTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCAtt                                  <  2:660184/139‑1 (MQ=255)
                                                                                           gaTCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACg                                             >  2:50739/1‑125 (MQ=255)
                                                                                           gaTCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACg                                             <  1:50739/125‑1 (MQ=255)
                                                                                           gaTCCACGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATc                               >  2:361202/1‑139 (MQ=255)
                                                                                                 cGTTGATGCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACCATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTGATTATCAACCTg                         >  2:555442/1‑139 (MQ=255)
                                                                                                        gCCATCAAACCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCG‑‑TGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGAtctc                >  2:918821/1‑139 (MQ=255)
                                                                                                                aCCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGtt                                      >  1:735376/1‑111 (MQ=255)
                                                                                                                aCCGGGCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGtt                                      <  2:735376/111‑1 (MQ=255)
                                                                                                                   gggCGACAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGaa       <  2:370757/139‑1 (MQ=255)
                                                                                                                        aCAAAGCTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCg  <  2:12513/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                    
TTGGCTTTGTACCGGTCCGTAAACCGGGCAAACTGCCGCGTGAAACCATCAGTGAAACTTACGACCTGGAATACGGCACCGATCAGCTGGAGATCCACGTTGATGCCATCAAACCGGGCGACAAAGTTCTGGTGGTGGACGACCTGCTGGCAACCGGCGGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCG  >  NC_000913/491650‑491908

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