Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 628,974 C→G 58.0% intergenic (‑41/‑68) ybdH ← / → ybdL putative oxidoreductase/methionine aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273628,9740CG58.0% 9.0 / 18.3 19intergenic (‑41/‑68)ybdH/ybdLputative oxidoreductase/methionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base G (6/5);  total (11/8)
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

GTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCACCCAGATGAGCGCGCTGGCG  >  CP009273/628834‑629122
                                                                                                                                              |                                                                                                                                                    
gttgttCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGGTAGGAGTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAACCGTCTTGCTTGcc                                                                                                                                             >  1:80629/1‑150 (MQ=255)
              aaaaaaTCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTGTTGATTGCCAAATACTAGCATCg                                                                                                                               >  1:132658/1‑150 (MQ=255)
                    tCGTGCAGGTGATCGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTtat                                                                                                                         <  1:192405/150‑1 (MQ=255)
                          aGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGcc                                                                                                                   >  2:246832/1‑150 (MQ=255)
                                    aaCTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAAATTTAGGTGTTTGCCACCCATTCGTAATGTAGACATCTCAACTTCTTGATTGCCGAATCCTAGCATCGTGTTATAGTGCCTTCAACACGc                                                                                                         >  2:262962/1‑150 (MQ=255)
                                                              gCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTCCAATAAAAATGa                                                                               >  2:24642/1‑150 (MQ=255)
                                                                    ccGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATa                                                                         <  1:261716/150‑1 (MQ=255)
                                                                       acgacCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAAcc                                                                      <  2:120542/150‑1 (MQ=255)
                                                                                      ggATTGTGAGGCATGTTGTGTTGACCTTATTGTTTTATTTAGACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACaaa                                                       >  1:218917/1‑150 (MQ=255)
                                                                                      ggATTGGGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACaaa                                                       <  2:45670/150‑1 (MQ=255)
                                                                                              aGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACtt                                               <  2:216040/150‑1 (MQ=255)
                                                                                                  aTGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCcac                                           >  1:153993/1‑150 (MQ=255)
                                                                                                        gtgttgtCCTTATTGTTTTATTTAGACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACtt                                      <  2:118324/150‑1 (MQ=255)
                                                                                                                       gTTTTATTTAGACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCACCACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTAtttt                        >  2:292525/1‑150 (MQ=255)
                                                                                                                           tatttaGACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCAcc                    <  2:111440/150‑1 (MQ=255)
                                                                                                                           tatttaGACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCAcc                    <  2:289422/150‑1 (MQ=255)
                                                                                                                                 gACATCTAAACGTGTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGCCACCATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCACCCaagaa              >  2:317749/1‑146 (MQ=255)
                                                                                                                                    aTCTAAACGTGTTGATTGCCAAATACTAGCATCGTGCTATAGTGCCTTCAACCCACAACTTCGTTAGGTACCATAAAAATGACAAATAACCGTCTGATTCCACAAAGCAAACTTCCACAACTTGGCGCCACTATTTTCACCCAGAAGAgc           >  1:299061/1‑150 (MQ=255)
                                                                                                                                             tCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCACCCAGATGAGCGCGCTGgcg  >  1:11412/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCACCCAGATGAGCGCGCTGGCG  >  CP009273/628834‑629122

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