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,8980CT97.6% 130.2 / ‑4.4 42A312V (GCC→GTC) uupATP‑binding protein with possible role in replication
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/22);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

GTGGATCCGCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGCCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGAATGGGTGCG  >  NC_000913/1010767‑1011032
                                                                                                                                   |                                                                                                                                      
gtgGATCCGCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCt                                                                                                                                <  2:225373/140‑1 (MQ=255)
      ccGCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTa                                                                                                                          >  2:314592/1‑140 (MQ=255)
        gCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTaaa                                                                                                                        >  2:350380/1‑140 (MQ=255)
        gCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTaaa                                                                                                                        >  2:586508/1‑140 (MQ=255)
                aTCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGtttt                                                                                                                <  1:350380/140‑1 (MQ=255)
                  cAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCg                                                                                                              <  1:642197/140‑1 (MQ=255)
                          gCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGaa                                                                                                      >  1:1357717/1‑140 (MQ=255)
                            cGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGAGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGa                                                                                                    <  2:1357717/140‑1 (MQ=255)
                             gTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGAc                                                                                                   <  1:314592/140‑1 (MQ=255)
                                    tAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAAGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCt                                                                                            >  2:686243/1‑140 (MQ=255)
                                      aTGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTAc                                                                                          <  2:359222/140‑1 (MQ=255)
                                      aTGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTAc                                                                                          <  2:1138681/140‑1 (MQ=255)
                                      aTGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTAc                                                                                          <  1:586508/140‑1 (MQ=255)
                                           ggCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGt                                                                                     >  1:1335498/1‑140 (MQ=255)
                                              cgcgTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGt                                                                                                 >  2:1373804/1‑125 (MQ=255)
                                              cgcgTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGt                                                                                                 <  1:1373804/125‑1 (MQ=255)
                                              cgcgTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTaa                                                                                  <  1:51653/140‑1 (MQ=255)
                                               gcgTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAAc                                                                                 <  1:1221622/140‑1 (MQ=255)
                                                     gcgcCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAg                                                                           >  2:46148/1‑140 (MQ=255)
                                                           tGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTg                                                                     >  1:1279591/1‑140 (MQ=255)
                                                               aGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGa                                                                 >  2:694756/1‑140 (MQ=255)
                                                                   aTGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGa                                                             >  1:1016751/1‑140 (MQ=255)
                                                                   aTGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGa                                                             >  2:1251186/1‑140 (MQ=255)
                                                                   aTGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGa                                                             >  1:1235665/1‑140 (MQ=255)
                                                                             aaCGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGcc                                                   <  1:694756/140‑1 (MQ=255)
                                                                             aaCGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGcc                                                   <  1:394241/140‑1 (MQ=255)
                                                                                gtggtgAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAg                                                <  2:1279591/140‑1 (MQ=255)
                                                                                         gTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGt                                       >  2:1051229/1‑140 (MQ=255)
                                                                                                       ggTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGc                         <  2:1335498/140‑1 (MQ=255)
                                                                                                       ggTACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGc                         >  1:163622/1‑140 (MQ=255)
                                                                                                        gTACCGCAAAGATGCAGGTGGAAGAAGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGcc                        <  1:686243/140‑1 (MQ=255)
                                                                                                         tACCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGccc                       <  2:1235665/140‑1 (MQ=255)
                                                                                                          aCCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGc                         <  2:1428462/137‑1 (MQ=255)
                                                                                                          aCCGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGc                         >  1:1428462/1‑137 (MQ=255)
                                                                                                            cGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGTAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGa                    >  2:1465051/1‑140 (MQ=255)
                                                                                                            cGCAAAGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGa                    >  1:727723/1‑140 (MQ=255)
                                                                                                               aaaGATGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTg                 <  1:580872/140‑1 (MQ=255)
                                                                                                                    tGCAGGTGGAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCg            <  2:163622/140‑1 (MQ=255)
                                                                                                                           ggAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGTAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGAATGGGt     <  1:1465051/140‑1 (MQ=255)
                                                                                                                            gAAGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGAATGGGTg    <  2:727723/140‑1 (MQ=255)
                                                                                                                             aaGAGGTCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGAATGGGTGc   >  1:784969/1‑140 (MQ=255)
                                                                                                                              agagGCCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGAATGGGTGCg  >  1:10397/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                      
GTGGATCCGCCAGGGGATCAAAGCACGCCGTACCCGTAATGAAGGCCGCGTACGCGCCCTGAAAGCGATGCGTCGCGAACGTGGTGAACGTCGCGAAGTGATGGGTACCGCAAAGATGCAGGTGGAAGAGGCCAGCCGCTCCGGTAAAATCGTTTTCGAAATGGAAGACGTTTGCTACCAGGTTAACGGTAAGCAACTGGTGAAAGATTTTTCTGCCCAGGTTCTACGTGGCGACAAAATTGCCCTGATTGGTCCGAATGGGTGCG  >  NC_000913/1010767‑1011032

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