Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 790,307 A→G 16.0% V241V (GTT→GTC galT ← galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913790,3070AG16.0% 61.3 / 5.0 25V241V (GTT→GTCgalTgalactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (12/9);  new base G (2/2);  total (14/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 = 7.84e-01

CCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTCGCGCTGAACATAATCCACCAGCATTGGTGATTTCTGTTCGGCAAAATATTCTTTTTGCAGGCGGTCTTCGCGCTCAGCTTCGTTAGGCAGGAAGCTATTTGCCCAAATCTGACCGTGCGGATGCGGGTTAGAGCAGCCCATCGCC  >  NC_000913/790090‑790562
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                               
ccaTAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGttt                                                                                                                                                                                                                                     >  1:15493/1‑246 (MQ=255)
ccaTAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCAAGCCCATAGAGTAGGGGAAGGAGCACTGGACGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATAGGTGATCCGTAAAACGTGGGCtt                                                                                                                                                                                                                                                       >  1:92247/1‑228 (MQ=255)
                                     cGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGgcagc                                                                                                                                                                                                                                               >  1:201004/1‑199 (MQ=255)
                                        aTAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAGACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGc                                                                                                                                                                                                                             >  2:252276/1‑214 (MQ=255)
                                        aTAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGt                                                                                                                                                                                                                                            <  2:201004/199‑1 (MQ=255)
                                          aaaaGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGc                                                                                                                                                                                                                             >  2:276517/1‑212 (MQ=255)
                                           aaaGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAGACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAg                                                                                                                                                                                                                          <  1:252276/214‑1 (MQ=255)
                                             aGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAg                                                                                                                                                                                                                          <  1:276517/212‑1 (MQ=255)
                                                                       ttCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGGTCGGGTTCGACa                                                                                                                                                                             >  1:115138/1‑231 (MQ=255)
                                                                                                           aTAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGc                                                                                                                                                                                                   >  2:76765/1‑173 (MQ=255)
                                                                                                             agagTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGt                                                                                                                                                                        >  1:6605/1‑198 (MQ=255)
                                                                                                              gagTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTaa                                                                                                                                                                                                <  1:76765/173‑1 (MQ=255)
                                                                                                                    gggAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGttt                                                                                                                                                                                  >  1:203210/1‑181 (MQ=255)
                                                                                                                       aaGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGa                                                                                                                                                                               <  2:203210/181‑1 (MQ=255)
                                                                                                                          gAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTcg                                                                                                                                                  >  2:202335/1‑207 (MQ=255)
                                                                                                                             cACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTCgcgc                                                                                                                                               <  1:202335/207‑1 (MQ=255)
                                                                                                                                                      tGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTCGCGCTGAACATAATCCACCAGCATTGGTGATTTCTGTTCGGCAAAATATTCTTTTTGCAGGCGGTCTTcgc                                                                            <  1:224083/249‑1 (MQ=255)
                                                                                                                                                       ggTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTaaa                                                                                                                                                                                                                                                                 >  2:131015/1‑67 (MQ=255)
                                                                                                                                                          cAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACg                                                                                                                                                                                                                                                              <  1:131015/67‑1 (MQ=255)
                                                                                                                                                                                       tgcgctgGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGc                                                                                                                                                                                                                             >  2:74319/1‑71 (MQ=255)
                                                                                                                                                                                          gctgGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAg                                                                                                                                                                                                                          <  1:74319/71‑1 (MQ=255)
                                                                                                                                                                                                    ggTCAAATCGGTGATCCGTAAGACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTCGCGCTGAACATAATCCACCAGCATTGGTGATTTCTGTTCGGCAAAATATTCTTTTTGCAGGCGGTCTTCGCGCCCAGCTTCGTTAGGCAGGAAGCTATTTGCCCAAATCTGACCGTGCGGATTCGGGTTtgacacgccactc     >  2:229712/1‑264 (MQ=255)
                                                                                                                                                                                                      tCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCt                                                                                                                                                    >  2:101406/1‑129 (MQ=255)
                                                                                                                                                                                                        aaaTCGGTGATCCGTAAGACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTCGCGCTGAACATAATCCACCAGCATTGGTGATTTCTGTTCGGCAAAATATTCTTTTTGCAGGCGGTCTTCGCGCTCAGCTTCGTTAGGCAGGAAGCTATTTGCCCAAATCTGACCGTGCGGATGCGGGTTAGAGCAGCCCATcgcc  <  1:229712/273‑1 (MQ=255)
                                                                                                                                                                                                         aaTCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTcgc                                                                                                                                                 <  1:101406/129‑1 (MQ=255)
                                                                                                                                                                                                                     gTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTCGCGCTGAACATAATCCACCAGCATTGGTGATTTCTGTTCGGCAAAATATTCTTTTTGCAGGCGGTCTTCGCGCTCAGCTTCGTTAGGCAGGAAg                                                     >  2:146200/1‑209 (MQ=255)
                                                                                                                                                                                                                        aaaCGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTCGCGCTGAACATAATCCACCAGCATTGGTGATTTCTGTTCGGCAAAATATTCTTTTTGCAGGCGGTCTTCGCGCTCAGCTTCGTTAGGCAGGAAGCTa                                                  <  1:146200/209‑1 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                               
CCATAAATTTACGTACGGTGGCGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCATAGAGTAGGGGAAGGAGCACTGGAAGAGGTTGTCATAACGACTGGTCAGCTTTTTCAACGCCAGCGCCAGATCGCTGCGCTGGGCGTCGGTCAAATCGGTGATCCGTAAAACGTGGGCTTTGGGCAGCAGTAGCGTTTCGAACGGCCAGGCAGCCCAGTAAGGCACGACGGCTAACCAGTGTTCGGTTTCGACAACGGTACGGCTACCGTCTGCCAGCTCGCGCTGAACATAATCCACCAGCATTGGTGATTTCTGTTCGGCAAAATATTCTTTTTGCAGGCGGTCTTCGCGCTCAGCTTCGTTAGGCAGGAAGCTATTTGCCCAAATCTGACCGTGCGGATGCGGGTTAGAGCAGCCCATCGCC  >  NC_000913/790090‑790562

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