Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 986,184 C→A 25.4% intergenic (+50/‑200) U712_04810 → / → U712_04815 putative oxidoreductase yhdF/putative amino acid permease YhdG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881986,1840CA25.4% 72.3 / 2.7 39intergenic (+50/‑200)U712_04810/U712_04815putative oxidoreductase yhdF/putative amino acid permease YhdG
Reads supporting (aligned to +/- strand):  ref base C (8/18);  new base A (9/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 9.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAAT  >  CP006881/986057‑986310
                                                                                                                               |                                                                                                                               
ctgctgGCGTCTGACGAATCTTCCTATATGACAGGGCAGAACATTCATGTGAATTGCGGCCGTTTTATTTAAAAGTAATAAAAGAAAAAACACCAAAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCGAAATaaaaaa                                                                                                                      >  6:118794/1‑139 (MQ=255)
                               cAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAACGAAAAACCACCTCAAGAGCTGGTTTTTTTTGTGGTGAAAAATTTAAAAAAAAAAAATTTTTCGA‑TTTTTATAACAATTTTGttttaa                                                                                      >  8:18071/1‑137 (MQ=255)
                                  ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAAAGAAAAACCACCAAAAGAGCTGGTTTTTTTTTTGGTGTAAAATTTAAAAATAAAAAATTTTTAGT‑TTTTTGTAAAAATTTTTTTTTTAttt                                                                                   >  6:139936/1‑139 (MQ=255)
                                  ggCAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACCTAAACAACCAAAAAACACCTCAAGAGCGGGTTTTTTTTTTGGTGAAAAATTTAAAAAAAAAAAATTTTTCGA‑TTTTTATAAAAATTTTGTTTTAAttt                                                                                   >  3:236235/1‑139 (MQ=255)
                                  ggCAGACCATTCATGTGAATGGCGGCCGTGTTTTTTTAAAGTAATAAACGAAAAACCACCTCAACAGGTGGTGTTTTTTGTGGTGGAAAATTTAAAAATAAAAAAATTTTAGATTTTTTTTAAAAATTTTGttttattt                                                                                    >  1:135990/1‑135 (MQ=255)
                                  ggCAGACCATTCATGTGAATGGCGGCCGTGTTATTTAAACGTAATGAACGAAAAAACAGCTCAAGAGCGGGGTTTTTGTGTGTGGTAAAATTTGAAAATATATAATTTTTAGA‑TTTTTATATAAATTTTGTTTTTAttt                                                                                   >  8:238189/1‑139 (MQ=255)
                                    cAGACCATTCATGTGAATGGCGGCCGTTTTATTTAAACGTAATAAAAGAAAAAACAGCTCAAAAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAAATTGTCGA‑TTTTTCTAAAAAATTTTTTTTAAttttt                                                                                 >  8:85564/1‑139 (MQ=255)
                                    cAGACCATTCATGTGAATGGCGGCCGTGTTTTTTAAAACTAATAAAAGAAAAAACAACAAAAAAGCAGGGTTTTTTTGTGTTGTAAAAATTAAAAATAAAAAAATTTTAGA‑TTTTTTTATAAATTTTTTTTTAAttttt                                                                                 >  2:58463/1‑139 (MQ=255)
                                       aCCATTCATGTGAATGGCGGCCGTTTTATTTAAACCTAATCAACGAAAAACCCGCTCAAGAGCTGGTTTTTTTTGTGGTGAAAAATTTAAAAAAAAAAAATTTGTCGA‑TTTTTCTTACAATTTTGTTTTCAttttttg                                                                               >  2:35489/1‑136 (MQ=255)
                                            tCATGTGAATGGCGGCCGTTTTATTTCAACGTAAAAAACGAAAAACCAGCTCAAGAGCGCGGTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAAATTGTCGA‑TTTTTATTACAAATTTTTTTTAAATTTTGCAAtttt                                                                         >  5:244021/1‑139 (MQ=255)
                                            tCATGTGAATGGCGGCCGTTTTATTTAAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTTTAGA‑TTTTTGTAACAATTTTTTTTTAATTTTTGAAAtttt                                                                         >  6:147073/1‑139 (MQ=255)
                                            tCATGTGAATGGCGGCCGTGTTTTTTCTACGGAATCACCCAAAAACCACCGCAAGAGCTGGGTTTTTGGGTGGGGGAAATATTGAAAATAAAAAAATTGTGGG‑TTTTTGTTACAAGTTTGTTTTGAGTTTTGGAAtttt                                                                         >  2:197559/1‑139 (MQ=255)
                                             cATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAAtttt                                                                         <  3:89578/138‑1 (MQ=255)
                                              aTGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtat                                                                       <  8:6424/139‑1 (MQ=255)
                                                 tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                    <  2:135990/139‑1 (MQ=255)
                                                 tgAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                    <  4:118255/139‑1 (MQ=255)
                                                  gAATGGCGGCAGTTTTATTTCAACGTAATCAACGAAAATACAGCTCAAGAGCTGGTTTTTTGTGTGTTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTAATTTTTGCAATTTtatata                                                                    <  5:124160/138‑1 (MQ=255)
                                                   aaTGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                    <  2:12946/137‑1 (MQ=255)
                                                   aaTGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTtatata                                                                    >  1:12946/1‑137 (MQ=255)
                                                     tGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTAAAAATAAAAAATTTGTCGA‑TTTTTATAACAATTTTGTTTTAATTTTTGCAATTTTATATttaa                                                                 >  8:276194/1‑138 (MQ=255)
                                                       gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                               <  3:284062/138‑1 (MQ=255)
                                                       gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                               <  7:85564/138‑1 (MQ=255)
                                                       gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                               <  5:140939/138‑1 (MQ=255)
                                                       gcggcCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                               <  5:146309/138‑1 (MQ=255)
                                                       gaggcCGTGTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATataata                                                               <  7:245140/136‑1 (MQ=255)
                                                           ccGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                          <  5:147073/139‑1 (MQ=255)
                                                                        aCGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGaa                                             >  4:67889/1‑139 (MQ=255)
                                                                         cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                          <  7:9682/125‑1 (MQ=255)
                                                                         cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAg                                                          >  8:9682/1‑125 (MQ=255)
                                                                         cGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAAt                                            >  1:257788/1‑139 (MQ=255)
                                                                                    aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                 <  5:139936/139‑1 (MQ=255)
                                                                                    aaaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                 <  3:226310/139‑1 (MQ=255)
                                                                                     aaaaCCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGAcc                                 <  7:18071/138‑1 (MQ=255)
                                                                                                       ggTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGaa               >  3:94124/1‑138 (MQ=255)
                                                                                                       ggTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGaa               >  3:36008/1‑138 (MQ=255)
                                                                                                                  tggtgAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAAAAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGATAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaata    <  6:244021/138‑1 (MQ=255)
                                                                                                                    gttaaaaaTTTAAAAAAAAAAATTTTGTCTA‑TTTTTCTAAAATTTTTGTTTTAATTTTTGAAATTTTATATATAATATTGAGTAGGAGAATGGAATTATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataat  <  1:148250/135‑1 (MQ=255)
                                                                                                                    gtgAAAATTTTAAAAATAAAAAATTTGTCGA‑TTTTTCTAAAATTTTTTTTTTAATTTTTTAAATTTTATATATAATTTTGAGTAGGATAAAGGATTGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataat  <  4:236235/138‑1 (MQ=255)
                                                                                                                    gtgAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTAATTTTTGAAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGaataat  <  7:276194/138‑1 (MQ=255)
                                                                                                                               |                                                                                                                               
CTGCTGGCGTCTGACGAATCTTCCTATATGACAGGGCAGACCATTCATGTGAATGGCGGCCGTTTTATTTCAACGTAATCAACGAAAAACCAGCTCAAGAGCTGGTTTTTTGTGTGGTGAAAAATTTCAAAATAAAAAATTTGTCGA‑TTTTTCTAACAATTTTGTTTTCATTTTTGCAATTTTATATATAATATTGAGTAGGAGAAAGGAATGATTCGCGACCAGTACAAATATAGTAGAAAAAAAAGAATAAT  >  CP006881/986057‑986310

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