Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 791,005 G→A H9Y (CAT→TAT)  galT ← galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913791,0050GA96.7% 87.3 / NA 30H9Y (CAT→TAT) galTgalactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/15);  minor base C (1/0);  total (15/15)
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

GAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTACGCGCCAGTTCAGTTC  >  NC_000913/790868‑791131
                                                                                                                                          |                                                                                                                              
gagGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCTATTGCCCCGTGAGCGGTTTGTAGCGGCGATGTGgataat                                                                                                                              >  2:319991/1‑140 (MQ=255)
      gCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACgg                                                                                                                        <  1:736573/140‑1 (MQ=255)
             ggATCGTGCGCAGGTTACACCTGTTTGGCTGGCGTTTGCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTaaa                                                                                                                 <  1:252494/140‑1 (MQ=255)
             ggATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTaaa                                                                                                                 >  2:662639/1‑140 (MQ=255)
                   tGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGt                                                                                                           >  2:621648/1‑140 (MQ=255)
                   tGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTAGTAGCGGCGATGTGGAGAATCCACGGGATTAAATTGCGt                                                                                                           >  2:229566/1‑140 (MQ=255)
                    gcgcAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTAAGCGGGTTTTAGCGGCGATGTGGATCATCAACGGGATTAAATTGCGTc                                                                                                          >  2:252494/1‑140 (MQ=255)
                               acCTGTTTGGCTGGCGTTTCCTCCGTCCCCTGCCAGGTGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTcc                                                                                               <  1:188075/140‑1 (MQ=255)
                                    tttGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGTGATGTGGAGAATTAACGGGATTAAACTGCGTCATGGTCGTTCCTTAAt                                                                                          >  2:269286/1‑140 (MQ=255)
                                    tttGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAAt                                                                                          >  2:589268/1‑140 (MQ=255)
                                    tttGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAAt                                                                                          >  1:534415/1‑140 (MQ=255)
                                    tttGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAGGGGGTTATTCAACGGGTTAAAATGGCGCAGTGGCCGTCCCTTAAt                                                                                          >  1:688636/1‑140 (MQ=255)
                                       ggctggcGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCgg                                                                                       >  1:512366/1‑140 (MQ=255)
                                         gtgggGTTTCCTGCGCCCCCTGGCCGGGGCGGTTTAGGGGGGGGGGGGAAACCAGGATTCCTTTCCCGGGTGA‑CGGGTTGTTGGGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGa                                                                                     <  2:688636/139‑1 (MQ=255)
                                         ctggcGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGa                                                                                     <  2:512366/140‑1 (MQ=255)
                                                  ccTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCtgt                                                                            >  2:577382/1‑140 (MQ=255)
                                                     gcgcCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAGCCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGa                                                                         <  1:662639/140‑1 (MQ=255)
                                                     gcgcCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGa                                                                         <  1:621648/140‑1 (MQ=255)
                                                                        cTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgccagtgcc                                                      >  1:172380/1‑140 (MQ=255)
                                                                         ttAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGGCGTTCCTTAATCGGGAGATCCCTGTTGATGGCGTGACtgccagtgcca                                                     >  2:215461/1‑140 (MQ=255)
                                                                                      gtgAAACCA‑GAATCCATTGCCCTGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGcc                                        <  1:269286/140‑1 (MQ=255)
                                                                                      gtgAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGcc                                        <  2:534415/140‑1 (MQ=255)
                                                                                                      aTTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGtgtg                         >  1:503792/1‑140 (MQ=255)
                                                                                                             ggTGAGCGGGTTGTAGCGGCGATGTGGATAATCACCGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTa                  <  1:710984/140‑1 (MQ=255)
                                                                                                             ggTGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTa                  <  2:172380/140‑1 (MQ=255)
                                                                                                               tGAGCGGGTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTAcg                >  1:376584/1‑140 (MQ=255)
                                                                                                                  gCGGGTTGTAGCGGCGATGTGGCTACTCACCGGGATTCAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTACgcgc             <  1:280387/140‑1 (MQ=255)
                                                                                                                  gCGGGTTGTAGCGGCGAGGTGGATAATCAACAGGATTAATTTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTACgcgc             <  2:270983/140‑1 (MQ=255)
                                                                                                                  gCGGGTTGTAGCGGCAATGTGGTTACTCATCGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTACgcgc             <  1:229566/140‑1 (MQ=255)
                                                                                                                    gggTTGTAGCGGCGATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGTGCCATTTCATCGAGTGTGCGCGTTACGCGCca           >  1:116957/1‑140 (MQ=255)
                                                                                                                     ggTTGTAGCGGCGATGTGGATAATCTACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTACGCGCcag          <  1:319991/140‑1 (MQ=255)
                                                                                                                             cggcgATGTGGATAATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGAGTCCTGCGGCATCTCATCGAGTGTGCGCGTAACGCGCCagttcagttc  >  2:245164/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                              
GAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCA‑GAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCCAGGTGTCCTGCGCCATTTCATCGAGTGTGCGCGTTACGCGCCAGTTCAGTTC  >  NC_000913/790868‑791131

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