Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,237,235 C→T Q173Q (CAG→CAA mglC ← D‑galactose/methyl‑galactoside ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,237,2350CT94.7% 48.4 / NA 20Q173Q (CAG→CAAmglCD‑galactose/methyl‑galactoside ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/11);  minor base G (0/1);  total (8/12)
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

CCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGG  >  NC_000913/2237002‑2237444
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                 
ccGGCTTTTAAAACCCCGCAAAGGGCAGAGAACACGCCAGACAACGCGTAGATCATCAGCAGTTTCAGGCCGACGTTGACACCAGATACTTTGGCCGTTTCCGGTTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGATCAAGGGTAGAGaa                                                                                                                                                                                                  <  2:217036/251‑1 (MQ=255)
   gCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGACGCTTACGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGATTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGAGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAggacaagaagcc                                                                                                                                                                                               >  1:54090/1‑241 (MQ=255)
           cATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAAACGGAAAATCCCCAGCGCGACAAAGCCATGAGCAAAGGTAGAGAAGCCACTgcca                                                                                                                                                                                        >  2:105679/1‑247 (MQ=255)
                               acacGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGcc                                                                                                                                                                   >  1:59508/1‑251 (MQ=255)
                                                                   ggCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAg                                                                                                                               <  1:57926/251‑1 (MQ=255)
                                                                                cccAGATCCTTTTGCCGCGTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCg                                                                                                                  <  2:106076/250‑1 (MQ=255)
                                                                                             gCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCGCCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGt                                                                                                     <  1:149441/251‑1 (MQ=255)
                                                                                                     cGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGAGATGTAAGCGAGAAGGAAACTCCCCAGCGCGACACAGCATTGAGCAACGGTAGAGAAGCCACTGTCAAACCAAGAAATTGGAGACGCCCCAACAAAGTCGAAATAGAGCGCGTTGATGCCATACACGATGATCATCGAGCCGAGCGTGGTAATGAACg                                                                                             >  1:37109/1‑251 (MQ=255)
                                                                                                      gggTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGATACGGAAACGCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACgg                                                                                            <  2:131471/251‑1 (MQ=255)
                                                                                                       ggTTACCACCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGc                                                                                           >  1:270443/1‑251 (MQ=255)
                                                                                                        gTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCg                                                                                          <  1:286549/251‑1 (MQ=255)
                                                                                                                 aaTGGCAAAAATTTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTTAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCAACGTGCCGAGCGTGGTAATGAACGGCGTCACGTTga                                                                                 <  2:3803/251‑1 (MQ=255)
                                                                                                                                           gggTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAACACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATc                                                       >  1:300014/1‑251 (MQ=255)
                                                                                                                                                 tGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAAcc                                                  <  1:15911/250‑1 (MQ=255)
                                                                                                                                                               cAGACGAACGCCTCCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGAGCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCaag                                   <  2:274706/251‑2 (MQ=255)
                                                                                                                                                                  aCGAACGCCACCGCAATCAACGCGTAGAAGGAGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGc                                >  1:84756/1‑251 (MQ=255)
                                                                                                                                                                          cACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAgcca                         >  2:6541/1‑247 (MQ=255)
                                                                                                                                                                             cGCAATCAACGCGTAGAAGGTGATGT‑‑‑AGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACATACAATCAGaa                  <  1:253703/251‑1 (MQ=255)
                                                                                                                                                                               cAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAGAGGTAGAGAAGCCACTGTCAAAACGAGAAACTGGCGACGCCCCGACAAAGTCGTAATAGAGCGGGTTGATGCCAAGGACGAAGATCATCGTGCCGAGCGTGGTAATTAACGGCGGCACGGTGAGATAAACGATAATAAATCCGTTGGTCAAACCGATCACCGCACCATTGTCACAGACAATCAGa                   <  2:319175/251‑1 (MQ=255)
                                                                                                                                                                                                gTGATGTAAGAGAGACGGAAACACCCCAGCGCGACAAAGCCGTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAAGTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGATGATGCCATAGACGATGATCATCGTGCCGAGCGGTGTAATAAACGGCGTCACGTTGAGATAAGCGATAATAAGTCCGTTGATCAAACCGATCACCGCACCAATGGCACAGACTATCAGAATAACCAGCGCAATCgg  <  2:100336/251‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                 
CCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGG  >  NC_000913/2237002‑2237444

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: