Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 628,963 C→G 40.8% intergenic (‑30/‑79) 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,9630CG40.8% 10.6 / 30.1 27intergenic (‑30/‑79)ybdH/ybdLputative oxidoreductase/methionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (8/8);  new base G (6/5);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCACCC  >  CP009273/628824‑629105
                                                                                                                                           |                                                                                                                                              
gcgcgAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCt                                                                                                                                            <  2:264094/144‑1 (MQ=255)
    gagaATGTTGTTCATCAGTGAAAACATCGTGCAGGTGATTGAAACTTTCTGGATGTGAAAAGTAGTTCGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGAGATCTAAAc                                                                                                                                        <  2:354292/144‑1 (MQ=255)
    gagaAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGGTTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGACCGACGACCACGCGGATAGGAGTGTAAGGCCCGGCGTGTGGCCCATATTGATTTCTTTAGACATCTAAAc                                                                                                                                        >  2:159864/1‑144 (MQ=255)
    gagaAAGTTGTTCATAAAAAAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGAGATCTAAAc                                                                                                                                        <  2:196139/144‑1 (MQ=255)
    gagaAACTTGGTCATCAGCGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGTCATGTTGTGTTGTCCTTATTGTTTTATTTAGAGATCTAAACGTCtt                                                                                                                                   <  1:162008/149‑1 (MQ=255)
        aaGTTGTTAATAAGTGAAAAAATCGTGCAGGTGATTGAACTTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGAtt                                                                                                                               <  1:123020/149‑1 (MQ=255)
                  tcaGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGAGATCTAAACGTCTTGATTGCCAAATACt                                                                                                                     >  1:187276/1‑149 (MQ=255)
                    attGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGAGATCTAAACGTCTTGATTGCCAAATACTAg                                                                                                                   <  1:298464/147‑1 (MQ=255)
                         aaaaaTCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCg                                                                                                              >  1:6230/1‑149 (MQ=255)
                                       tGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGAGATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGccgt                                                                                                <  1:331571/149‑3 (MQ=255)
                                        gATTGAAACTTCCTGGACGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTc                                                                                               <  1:85169/149‑1 (MQ=255)
                                         aTTGAAACTTCCTGGATGAGAAAAGTAGTTAGCCGGGCGGACGACCACGCGGAAAGGATTGTGAGGCATGTTGTGGTGTCCTTATTGTTTTATTTAAAGAGCGAAACGTCTTGATTGCCAAATACTAGCATCGTGGTATTGTGCCttaa                                                                                              >  1:49896/1‑147 (MQ=255)
                                             aaaCTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAacac                                                                                          >  1:331904/1‑149 (MQ=255)
                                                                         ccGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGa                                                              <  1:275975/149‑1 (MQ=255)
                                                                            ggCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACaa                                                           >  1:228064/1‑149 (MQ=255)
                                                                                                  aTTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAAATTCGTCAGGTACAAAAAAAATGACAAATAACCCTCTGATTAca                                          >  2:289318/1‑144 (MQ=255)
                                                                                                   ttGTGAGGCATGTTGTGTTGTCCTTATTGATTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGc                                    >  1:242656/1‑149 (MQ=255)
                                                                                                        aGGCATGTTGTGTTGTCCTTATTGATTTATTTAGAGATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACt                               >  1:296626/1‑149 (MQ=255)
                                                                                                           cATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTcc                            <  1:54513/149‑1 (MQ=255)
                                                                                                                    ttgtCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGCTTCCACAAAGCAAACTTCCACaa                        >  2:349426/1‑144 (MQ=255)
                                                                                                                      gtCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCAACAAAGCAAACTTCCCCAACTTGGca                 >  1:163247/1‑149 (MQ=255)
                                                                                                                       tCCTTATTGTTTTATTTAGAGATCTAAACGTCTTGATTGCCAAATACTAGCATAGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACtt                     >  2:237733/1‑144 (MQ=255)
                                                                                                                         cTTATTGTTTGATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTgg                   <  2:337678/144‑1 (MQ=255)
                                                                                                                           tATTGTTTTATTTAGAGATCTAAACGTCTGGATTGCCAAATACTAGCATCGTGATATAGTGCCGTCAACACGCAACTTCGTCAGGTACAATAAAGACGACAAATAACACTCTGATTACACAAAGCAAACTTCCACAACTTGGca                 >  2:192799/1‑144 (MQ=255)
                                                                                                                             ttGTTTTATTTAGAGATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCGTCAACACGCAACTTCGTCAGGTACAATAAACATGACAAATAACCCTCCGATTCCACAAAGCACACTTCAACACCTTGGCACAACTAt          >  1:53532/1‑149 (MQ=255)
                                                                                                                                ttttATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTAtttt       <  1:312678/149‑1 (MQ=255)
                                                                                                                                          aCATCTAAACGTCTTGATTGCCAAATACGAGCCTCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCAccc  <  2:5411/144‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GCGCGAGAAAGTTGTTCATCAGTGAAAAAATCGTGCAGGTGATTGAAACTTCCTGGATGTGAAAAGTAGTTAGCCGGGCCGACGACCACGCGGATAGGATTGTGAGGCATGTTGTGTTGTCCTTATTGTTTTATTTAGACATCTAAACGTCTTGATTGCCAAATACTAGCATCGTGTTATAGTGCCTTCAACACGCAACTTCGTCAGGTACAATAAAAATGACAAATAACCCTCTGATTCCACAAAGCAAACTTCCACAACTTGGCACCACTATTTTCACCC  >  CP009273/628824‑629105

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