Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,366,668 C→T A36T (GCC→ACC)  dcuA ← C4‑dicarboxylate transporter DcuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,6680CT97.0% 98.6 / NA 34A36T (GCC→ACC) dcuAC4‑dicarboxylate transporter DcuA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/14);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

AGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAG  >  NC_000913/4366426‑4366915
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
aGCGATAACTGGCAGTGTAGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGCTCGGTGCGAGGATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAAcc                                                                                                                                                                                                                                                   >  2:513492/1‑250 (MQ=255)
 gCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCa                                                                                                                                                                                                                                                  >  2:426380/1‑250 (MQ=255)
                       gagagaTGTTGCCAGTAACAGCAAAGATAGTCAGGAAATAGGTCAAGATCGGTGCGAGGATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCaaa                                                                                                                                                                                                                           >  1:197236/1‑250 (MQ=255)
                                            cAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAAt                                                                                                                                                                                                     >  1:9324/1‑251 (MQ=255)
                                             aaaGATAGTCAGGAAATAGGTCACGATCGGTGCGAGAATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACAGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGCTAATGGAGATGACATCGAACGGGATGTTACCGGGTTTACCGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCactc                                                                                                                                                                                                    >  2:481449/1‑248 (MQ=255)
                                                gATAGTCAGGAAATAGGTCAAGATCGGTGAGAGGATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAAACAGATAGTCCAGACCGCCAGCAACCTGCATGGCAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGAAAGAACCAGAAC‑CCCCAATCCGCCTGCAACACCAATACATATTCCCCACAATCTgg                                                                                                                                                                                                 >  1:469042/1‑251 (MQ=255)
                                                                 tCACGATCGGTGCGAGGAACGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGcc                                                                                                                                                                                >  2:268765/1‑251 (MQ=255)
                                                                 ccacgttgggTGCGGGGATCGTGATGT‑TTTCGGTGTTGCGGCGCAGCAGCTTTTCTGCCGGATGAACAAAATACCCCACCCCGCCAACAACTTGCATGGAAAAAATAGGGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGcc                                                                                                                                                                                <  1:230547/243‑1 (MQ=255)
                                                                   aCGATCGGTGCGAGGATCGTGATGTATTTCGG‑GTTACCGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGcc                                                                                                                                                                                >  2:111026/1‑249 (MQ=255)
                                                                               aGGATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTatgat                                                                                                                                                                  >  2:3499/1‑251 (MQ=255)
                                                                                             tttCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACt                                                                                                                                                     >  1:459743/1‑250 (MQ=255)
                                                                                              ttCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACt                                                                                                                                                     >  2:196701/1‑249 (MQ=255)
                                                                                                cgggtttgcGGCGCAGCATTTTTTCTGTTTGATAAACCAGATAGCCCAGCCCGCCAGCAACCGGCATGGAAAAAATAGGGGCAAAAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAg                                                                                                                                                   <  1:196701/242‑1 (MQ=255)
                                                                                                            gcagcaGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGccttc                                                                                                                                      >  1:232397/1‑251 (MQ=255)
                                                                                                             cagcagCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGccttcc                                                                                                                                     >  1:181381/1‑251 (MQ=255)
                                                                                                                                              cAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACTGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAACACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGCTCGCCAGCAAAACTATGATGAGTTCTACCACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGtttt                                                                                                    >  2:24287/1‑249 (MQ=255)
                                                                                                                                              cAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTCGCATATTAGCCTTCCTTGTTTTTTAACAAGTTAATATTCAATTGTTAt                                                                                                    >  2:563429/1‑251 (MQ=255)
                                                                                                                                                         cAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGCGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTAc                                                                                         >  1:113634/1‑251 (MQ=255)
                                                                                                                                                           actggcagggcagaaaagggggcaaTAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTg                                                                                       <  1:193612/227‑1 (MQ=255)
                                                                                                                                                                                    tAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCCC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCACAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                              >  1:275642/1‑251 (MQ=255)
                                                                                                                                                                                    aagccGCCATGATAATGGAGATGACATCGAACGGGATGTAACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                              <  1:32914/248‑1 (MQ=255)
                                                                                                                                                                                     aaCCTCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCt                                                             <  2:253554/251‑1 (MQ=255)
                                                                                                                                                                                          ccATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGt                                                        <  1:68630/251‑1 (MQ=255)
                                                                                                                                                                                           cATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTg                                                       <  2:9324/251‑1 (MQ=255)
                                                                                                                                                                                                  aggaaattgacgacaaccGGGATGTACCCGGGTTTAACCCCAATAGCGGTAAGAACCAGCAC‑CCCAAACCCGCCTGCAAACCCAATACCTATTCCCCC‑AATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCtttttt                                               <  2:109486/234‑1 (MQ=255)
                                                                                                                                                                                                          gACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACACGAAAAAAGGCCCGTCATCTGACGTGCCTTTTTTATTTGt                                         >  2:480994/1‑250 (MQ=255)
                                                                                                                                                                                                                     ggagggtACCGGGTTTAACGCCAATAGCTGTAAGAACCAGCCC‑CCCCAATCCGCCTGCAAAACCAATACCTATGCCCCCCAATCTGGCGCCCATGAAGATCGCCAGCAAAACTATGATGAGTTCTTCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTGCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCTGTTTTTATTTGTACTACCCTGTAc                             <  2:426705/245‑1 (MQ=255)
                                                                                                                                                                                                                           tACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCGATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTa                        <  1:357509/250‑1 (MQ=255)
                                                                                                                                                                                                                                ggTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACGACCCTGTACGATTACTGTTc                  <  1:171358/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTTAACGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCg                 <  1:426380/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTTAAAGCCAAGAGCGGTAAGACCCACCACTCCCCAA‑CCTCCTGCAAAACCAATACCTATGCCCCCCAATCTTGCTCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGAGTGTTATTTTTAAGTTACTGCTCGCAAGAAAAAAGGCACGTCACCGGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGATCg                 <  2:225046/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aCGCCAATAGCGGTAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtc            <  1:469915/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tAGCGGTAAGAACCAGCAC‑ACCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGt      >  2:463080/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  gCAAGAACCTGCAC‑CCCCAATACGCCTGCAAAACCAAACCCAATTACCCCCAATCTGGCCCCCAAGCAGATCGCCACAAAAACTATGACGAGTTCTACAACTCGCCTATTAGCCTTACTTGTTTTTTCACAAGATGATATTAGTTTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCCCGACATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGCTAACTGTTCGCTTTCCTCAGTATAg  >  2:348788/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
AGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGG‑GTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAC‑CCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAG  >  NC_000913/4366426‑4366915

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