Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,665,311 C→T 9.5% intergenic (‑57/‑463) yhjB ← / → yhjC putative DNA‑binding transcriptional response regulator/LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,665,3110CT 9.5% 119.0 / 3.8 42intergenic (‑57/‑463)yhjB/yhjCputative DNA‑binding transcriptional response regulator/LysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (20/18);  new base T (2/2);  total (22/20)
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.82e-01

TGCAGCAACCAGTAGCAAAACTCACCATCCTGATCGCCATCGAGCATAACTAAGGCTTCAGGGTAACTTTCCAGCTTTTGCCATAACTCGTCTGCCTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCACTATTTAAAACAATTTTGAGGATTTCCTTATATTGGTGGTTAGTACGCATGCAATTAAAAATGAAATTCCGCGACCACAAGCCAAAATAACAAACGGCAAGGAGACAAAAATAAGCACAAATAGCCAACAC  >  CP009273/3665058‑3665560
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
tGCAGCAACCAGTAGCAAAACTCACCATCCTGATCGCCATCGAGCATAACTAAGGCTTCAGGGTAACTTTCCAGCTTTTGCCATAACTCGTCTGCCTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAAATGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCg                                                                                                                                                                                                                                      <  1:357109/275‑1 (MQ=255)
                                   gCCATCGAGCATAACTAAGGCTTCAGGGTAACTTTCCAGCTTTTGCCATAACTCGTCTGCCTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTTCCTGTCAAAAATGACTATTTGCCTAACTGAATCTCCACCTCAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCTGCGAGAAAAAATAAGAATTCGCCATTTGGCTGTGGCCATTCTACagagtttaatt                                                                                                                                                                                                    >  2:112280/1‑267 (MQ=255)
                                       tCGAGCATAACTAAGGCTTCAGGGTAACTTTCCAGCTTTTGCCATAACTCGTCTGCCTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACgt                                                                                                                                                                                                    >  2:48353/1‑270 (MQ=255)
                                       tCGAGCATAACTAAGGCTTCAGGGTAACTTTCCAGCTTTTGCCATAACTCGTCTGCCTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACgt                                                                                                                                                                                                    <  1:48353/270‑1 (MQ=255)
                                                         tCAGGGTAACTTTCCAGCTTTTGCCATAACTCGTCTGCCTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGATTAAGATAAAAAGACTTAAGTAGTGNAAGGGTATTACACGCGAGAAAAAATAAGAATACGCCATTTGGCGCTGGCCATTCTACAGATATGACGTGTAGAAAATAGTAACAGATaaaa                                                                                                                                                                             >  1:630033/1‑275 (MQ=255)
                                                                               gCCATAACTCGTCTGCCTGACTGGCCCCCTGAATACTCCCTCCTGGAATCCGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAATAATTCGCCATTTGGCGTTGGCCATTCTACAGAGATGCCGTGTAGAAACTAGTTACCGATATAAATAGTTACAGCTAAACGCCTGa                                                                                                                                                       >  1:196594/1‑275 (MQ=255)
                                                                                              ccTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAgaga                                                                                                                                                                                                          >  2:86981/1‑209 (MQ=255)
                                                                                               cTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCTGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAggg                                                                                                                                       <  1:112280/275‑1 (MQ=255)
                                                                                                  aCTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGAc                                                                                                                                                                                                      <  1:86981/209‑1 (MQ=255)
                                                                                                                   tcactcCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATaaa                                                                                                                                                                            >  2:591931/1‑218 (MQ=255)
                                                                                                                     actcCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCAGCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATGACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTTTAGAAAATAATTACCGATATAAATAGTTACAGATaa                                                                                                                                                               >  2:312470/1‑229 (MQ=255)
                                                                                                                       tcCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGt                                                                                                                                                                        <  1:591936/218‑1 (MQ=255)
                                                                                                                        ccTGGAATACGCTGCTGTAAACTGATTTTCATTCCATTAATAAAGAATGACTGCCTGTCAAACATGTCTATTTTCATAACTTAATCTCCACCTGGAGACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGTCAGTTGGCGGTGGCCATGCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCACTATTTAAAACAATTTTGAGGa                                                                                                               <  1:440421/274‑1 (MQ=255)
                                                                                                                         cTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGc                                                                                                                                                           <  1:312474/229‑1 (MQ=255)
                                                                                                                            gAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATa                                                                                                                                                                                          >  2:88024/1‑195 (MQ=255)
                                                                                                                                aCGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTa                                                                                                                                                                                      <  1:88024/195‑1 (MQ=255)
                                                                                                                                 cGCTGCTGTAAACTGATTTTCATTCCATGAATAAATCTTGACTGCCTGTCAAACATGACTATTTGCATAACTGCATCTCCACCTGAATACGTTAAAAAGAGTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGACGAAAGTAGTTACCGATATAAATAATTACAGCAAAACGCCTGAAATTACATGTCGAGGGCAATATTTAAAACAATTTTAAGGATATCCTtata                                                                                                     >  2:781506/1‑275 (MQ=255)
                                                                                                                                 cGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATAACTATTTGCATCACTGGATCTCCCCCTGACTACGTGAAAAAGACTTAATTCGTGGTAGGGGATTACCCGCGAGAAAAAATAAGACTTCGCCATTTGGCGGTGGCCATCTTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGGTACAGCTAAACGCATGAAATNAGCTTTCGAGGCGACTATTNAAAAAAATTTTANANATTTACTttta                                                                                                     >  2:440416/1‑275 (MQ=255)
                                                                                                                                                                         aCTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGAACTAATTAAAACAATTTTg                                                                                                                   >  2:598905/1‑221 (MQ=255)
                                                                                                                                                                             ccTGTCAAACATGACTATTTGCATAACTGAATCTCCANCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCTAGGGCACTATTTAAAACAATTTTGAGGa                                                                                                               <  1:598910/221‑1 (MQ=255)
                                                                                                                                                                             ccTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAgaga                                                                                                                                                                                                          >  1:451730/1‑130 (MQ=255)
                                                                                                                                                                              cTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCACTATTTAAAACAATTTTGAGGATTTCCTtata                                                                                                     >  1:802026/1‑230 (MQ=255)
                                                                                                                                                                                 tCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGAc                                                                                                                                                                                                      <  2:451725/130‑1 (MQ=255)
                                                                                                                                                                                  cAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCACTATTTAAAACAATTTTGAGGATTTCCTTATATtgg                                                                                                 <  2:802017/230‑1 (MQ=255)
                                                                                                                                                                                          aCTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAgaga                                                                                                                                                                                                          >  1:654737/1‑117 (MQ=255)
                                                                                                                                                                                           cTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAg                                                                                                                                         >  2:447072/1‑181 (MQ=255)
                                                                                                                                                                                              tttGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGAc                                                                                                                                                                                                      <  2:654730/117‑1 (MQ=255)
                                                                                                                                                                                               ttGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCa                                                                                                                                     <  1:447077/181‑1 (MQ=255)
                                                                                                                                                                                                      aCTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCTGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGaaa                                                                                                                                                     >  2:192048/1‑158 (MQ=255)
                                                                                                                                                                                                          aaTCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCTGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTAc                                                                                                                                                 <  1:192051/158‑1 (MQ=255)
                                                                                                                                                                                                          aaTCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCACTATTTAAAACAATTTTGAGGATTTCCTTATATTGGTGGTTAGTACGCATGCAATTAAAAATGAAATTCATCGACCACAAACCAAAATAACatcagacaaaaaa                             >  2:399181/1‑262 (MQ=255)
                                                                                                                                                                                                                 aCCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCACTATTTAAAACAATTTTGAGGATTTCCTTATATTGGTGGTTAGTACGCATGCAATTAAAAATGAAATTCCGCGACCACAAGCCAAAATAACAAACGGCAAGGAGACAAAAAt                     >  2:272489/1‑275 (MQ=255)
                                                                                                                                                                                                                      aaTACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGt                                                                                                                                                                                        >  2:306966/1‑107 (MQ=255)
                                                                                                                                                                                                                          cGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTAcc                                                                                                                                                                                    <  1:306970/107‑1 (MQ=255)
                                                                                                                                                                                                                               aaaaGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCACTATTTAAAACAATTTTGAGGATTTCCTTATATTGGTGGTTAGTACGCATGCAATTAAAAATGAAATTCCGCGACCACAAGCCAAAATAACAAACGGCAAGGAGACAAAAATAAGCACAAATAGcc       <  1:272493/275‑1 (MQ=255)
                                                                                                                                                                                                                                aaaGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTg                                                                                                                                                                                                                         >  2:820400/1‑64 (MQ=255)
                                                                                                                                                                                                                                    aCTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCa                                                                                                                                                                                                                     <  1:820409/64‑1 (MQ=255)
                                                                                                                                                                                                                                     cttcatgagtGGAAGGGTATTACGCGCGAGAAAAAATAAGAATTCGCCATTTGGCGTTGGCAATTCTACAGAGATGACGTGTAGAAAATAGTTACAGATATAAATAGTTACAGCTAAAAGCCTGAAATTACATGTCGAGGGCACTAGTTATAACAATTTAGATGATTGCCTTAGATTGGTGGTTAGTACGCATGCAATTAAAAATGAAATTCCGCGACAACAAGCCAAAATAACAAACGGCAAGGAGACAAAAATAAGCACAAATAGCCAacac  <  2:196591/267‑1 (MQ=255)
                                                                                                                                                                                                                                       tAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCAt                                                                                                                                                                                                                    >  1:444943/1‑62 (MQ=255)
                                                                                                                                                                                                                                           tagtGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTa                                                                                                                                                                                                                <  2:444938/62‑1 (MQ=255)
                                                                                                                                                                                                                                                     gTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAg                                                                                                                                         >  2:803410/1‑123 (MQ=255)
                                                                                                                                                                                                                                                         tACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACACGTCGAGGGCa                                                                                                                                     <  1:803419/123‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
TGCAGCAACCAGTAGCAAAACTCACCATCCTGATCGCCATCGAGCATAACTAAGGCTTCAGGGTAACTTTCCAGCTTTTGCCATAACTCGTCTGCCTGACTGGCCCCCTGAATACTCACTCCTGGAATACGCTGCTGTAAACTGATTTTCATTCCATGAATAAATATTGACTGCCTGTCAAACATGACTATTTGCATAACTGAATCTCCACCTGAATACGTTAAAAAGACTTAAGTAGTGGAAGGGTATTACCCGCGAGAAAAAATAAGAATTCGCCATTTGGCGGTGGCCATTCTACAGAGATGACGTGTAGAAAATAGTTACCGATATAAATAGTTACAGCTAAACGCCTGAAATTACATGTCGAGGGCACTATTTAAAACAATTTTGAGGATTTCCTTATATTGGTGGTTAGTACGCATGCAATTAAAAATGAAATTCCGCGACCACAAGCCAAAATAACAAACGGCAAGGAGACAAAAATAAGCACAAATAGCCAACAC  >  CP009273/3665058‑3665560

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