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,6680CT96.2% 86.2 / NA 26A36T (GCC→ACC) dcuAC4‑dicarboxylate transporter DcuA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/12);  minor base G (1/0);  total (14/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

TTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCAT  >  NC_000913/4366423‑4366907
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
taaagagaaaaggggCAGGGTCGCCAGAGAGAGGTGGCGAGAACCAGAAAAGAATGTCAGGAAAT‑GGTGCAGGAGGGGGACGGGGATGGTTAGGTATTT‑GGGTTTCGGCCCAAGCAGCTTTGGTGTCGGATGAACGAGATGGCCCCGACCGCCAGCACCCTGCAGGGAAAAAATGGGGGCAATTACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGa                                                                                                                                                                                                                                             <  2:361044/237‑1 (MQ=255)
  cAGCGATAACGGGCAGGGGCGCCAGAGAGATGTTGCCAGTACCAGAAAAGAAAGTAAGGAAATGGGG‑CACGATGGGTGGGAGGATCGTGATGTATTTGGGGTTAGGGCGCG‑GCAGTTTTTTTGTTTGATGAACCAAAAAGCCCAGCCCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAAcc                                                                                                                                                                                                                                          <  1:9891/251‑1 (MQ=255)
       aTAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAc                                                                                                                                                                                                                                     >  1:120598/1‑251 (MQ=255)
        aaccgggCAGTGTCGCAAGAGAGATGTTGCCAGTACCAGCAAAGATAGTAAGAAAATAGGT‑CACGATGGGTGCGAGGATCGTGATGTATTTGGGGTTCGGGCGCA‑GCAGTTTTTCTGTCTGAGGAACCAGATAGCCCAGCCCGCCAGCAACTTGCAGGGAAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAcc                                                                                                                                                                                                                                    <  1:147003/246‑1 (MQ=255)
          aCTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTACAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCAcccc                                                                                                                                                                                                                                  >  1:317790/1‑251 (MQ=255)
                                 gTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGCGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATAc                                                                                                                                                                                                           >  2:262182/1‑251 (MQ=255)
                                        gTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCCCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTcc                                                                                                                                                                                                    >  2:390894/1‑251 (MQ=255)
                                          aCCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTcccc                                                                                                                                                                                                  >  1:154388/1‑251 (MQ=255)
                                          aCCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTcccc                                                                                                                                                                                                  >  1:161074/1‑251 (MQ=255)
                                            cAGCAAAGATAGTCAGGAAATAGGA‑CACGAACGGTGCGAGGATCGTGAAGTATTAAGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGAGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTCACGCCAGTAGCGGGAAGAACCAGCACCACCAATCCACCTGCAAAACAAAGACCTATTccccc                                                                                                                                                                                                 >  2:403181/1‑250 (MQ=255)
                                                    aTAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGgc                                                                                                                                                                                        >  1:72574/1‑251 (MQ=255)
                                                                                    ggATCGTGATGTATTTCGGGTTTCGGAGCA‑GAAGCATTTCTGTCTGATGAACCAGATAGTCCAGACCGACAGCAACCTGCATGGCAGAAATAGCGGCGCTAAACGACCTGATAATGGAGATGACATCGAACGGGATGGTACAGGGATTAACGCCAATAGCGGTAAGCACACGCACCCCAAATCCGCCTGCAAAACCCATACATATTCCCCCCCATATGGCGCCCAAGAAGATCGACAGCAAAACTatgatg                                                                                                                                                         >  2:316783/1‑251 (MQ=255)
                                                                                         gTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTc                                                                                                                                                    >  2:387490/1‑251 (MQ=255)
                                                                                                                                                                 aacaaggaaaaaaagagagacgATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCt                                                                             <  2:44839/231‑1 (MQ=255)
                                                                                                                                                                            aaaTAGGGGCAAAAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTTCTGCTCACAAGaaaaa                                                                  <  1:352589/251‑1 (MQ=255)
                                                                                                                                                                             aaaaGCGGCAAAAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTaa                                                                                                                  <  1:357175/199‑1 (MQ=255)
                                                                                                                                                                             aaTAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTaa                                                                                                                  >  2:357175/1‑202 (MQ=255)
                                                                                                                                                                                   ggcgATAACAGCCATGATAATGGAGATGACCTCGAACGAGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAAACAGCACCCCCAATCCGCCTGCAAAACCCATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAATTTACTGCTCACAAGaacacacgtacg                                                           >  2:428860/1‑241 (MQ=255)
                                                                                                                                                                                            cGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGACTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGAc                                                  >  1:342359/1‑251 (MQ=255)
                                                                                                                                                                                              ccATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGatat                                                                                                       >  2:349653/1‑196 (MQ=255)
                                                                                                                                                                                              ccATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGatat                                                                                                       <  1:349653/196‑1 (MQ=255)
                                                                                                                                                                                                         gagaTGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTGGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                     <  1:186225/251‑1 (MQ=255)
                                                                                                                                                                                                                cATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtact                              <  1:420102/251‑1 (MQ=255)
                                                                                                                                                                                                                        gggATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTa                      <  1:374089/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCg         <  1:390894/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat  <  2:94684/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
TTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCA‑GCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCAT  >  NC_000913/4366423‑4366907

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: ‑