Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,010,898 C→T A312V (GCC→GTC)  uup → ATP‑binding protein with possible role in replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,010,8980CT96.5% 87.8 / NA 29A312V (GCC→GTC) uupATP‑binding protein with possible role in replication
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/11);  minor base A (1/0);  total (18/11)
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

GCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGCCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGA  >  NC_000913/1010775‑1011023
                                                                                                                           |                                                                                                                             
gCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGAAGGGGGAAGAGGTCAGCCGCTc                                                                                                                      >  1:292376/1‑133 (MQ=255)
 ccAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTcc                                                                                                                     >  2:376504/1‑133 (MQ=255)
 ccAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTcc                                                                                                                     >  1:524335/1‑133 (MQ=255)
                aCGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTc                                                                                                      >  2:339270/1‑133 (MQ=255)
                   ccGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCTGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGaa                                                                                                   >  1:59527/1‑133 (MQ=255)
                          cGTAATGAGGGCCGCGTAAGCGCCCTGAAATCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAg                                                                                            <  1:211917/133‑1 (MQ=255)
                            tAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGAc                                                                                          >  1:608441/1‑133 (MQ=255)
                              aTGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGt                                                                                        <  2:429414/133‑1 (MQ=255)
                                       gcgTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCa                                                                               >  2:25608/1‑133 (MQ=255)
                                             gcgcCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTaa                                                                         <  2:525268/133‑1 (MQ=255)
                                                       aGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCaa                                                               >  2:76699/1‑133 (MQ=255)
                                                       aGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCaa                                                               <  2:608441/133‑1 (MQ=255)
                                                           aTGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTGTTCGAAAGGGAAGACGTTTGCTACCAGCTTAACGGTAAGCAACTgg                                                           >  2:434071/1‑133 (MQ=255)
                                                           aTGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGCGGAAGAGGACAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTgg                                                           >  2:377585/1‑133 (MQ=255)
                                                             gCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTg                                                         >  2:271594/1‑133 (MQ=255)
                                                                            tgAACGTCGCGAAGTGATGGGTACCGCAAAGATGCTGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGcc                                          <  2:59527/133‑1 (MQ=255)
                                                                            tgAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGcc                                          <  1:25608/133‑1 (MQ=255)
                                                                                  tCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTTGTGAAAGAGTTTTCTGCCCAGGtt                                    >  2:264539/1‑133 (MQ=255)
                                                                                  tCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGtt                                    >  1:42000/1‑133 (MQ=255)
                                                                                  tCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGtt                                    >  1:89781/1‑133 (MQ=255)
                                                                                                  aCCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACaaaa                    <  1:376504/133‑1 (MQ=255)
                                                                                                  aCCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACaaaa                    <  1:377585/133‑1 (MQ=255)
                                                                                                  aCCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACaaaa                    >  1:485642/1‑133 (MQ=255)
                                                                                                    cGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAAtt                  >  1:511666/1‑133 (MQ=255)
                                                                                                    cGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAAtt                  >  2:56524/1‑133 (MQ=255)
                                                                                                           aggCAGGTGTAAGAGGTCAGCCGCACGGGTAATATCGTTTTGGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGa           <  2:292376/131‑1 (MQ=255)
                                                                                                           aTGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGa           <  2:511666/133‑1 (MQ=255)
                                                                                                                ggTGGAAGTGGTCAGACGCTCCTGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGt      <  1:189248/133‑1 (MQ=255)
                                                                                                                    gAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGa  >  2:140980/1‑133 (MQ=255)
                                                                                                                           |                                                                                                                             
GCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGCCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGA  >  NC_000913/1010775‑1011023

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