Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 679,233 G→C 9.8% S256T (AGC→ACC)  U712_03320 → Putative lipoprotein yerB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881679,2330GC 9.8% 130.0 / 6.6 41S256T (AGC→ACC) U712_03320Putative lipoprotein yerB
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base C (2/2);  total (23/18)
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.29e-01

ATAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGAAACAAGTGAATGGCAGAATTGTG  >  CP006881/679095‑679357
                                                                                                                                          |                                                                                                                            
aTAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccag                                                                                                                              >  8:184406/1‑139 (MQ=255)
 tAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  8:170632/139‑1 (MQ=255)
 tAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  7:72020/139‑1 (MQ=255)
 tAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  2:125579/139‑1 (MQ=255)
 tAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  3:355607/139‑1 (MQ=255)
 tAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  6:371645/139‑1 (MQ=255)
 tAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  6:171637/139‑1 (MQ=255)
 tAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  1:171629/139‑1 (MQ=255)
     gTTACAAATCTTGTCGAATACAACTATGATAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCcatca                                                                                                                        >  5:268665/1‑139 (MQ=255)
           aaTCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCAtatta                                                                                                                   <  1:193625/139‑1 (MQ=255)
           aaTCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCAtatta                                                                                                                   <  7:254793/139‑1 (MQ=255)
               ttGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGa                                                                                                               >  6:264247/1‑139 (MQ=255)
                tGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGAt                                                                                                              >  5:135236/1‑139 (MQ=255)
                 gTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGAt                                                                                                              >  6:74109/1‑138 (MQ=255)
                  tCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATTCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCa                                                                                                            <  6:50577/139‑1 (MQ=255)
                               tGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCgg                                                                                               >  3:191064/1‑139 (MQ=255)
                                 aTAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGa                                                                                             >  5:363077/1‑139 (MQ=255)
                                    aaaaaGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACAt                                                                                          >  7:322190/1‑139 (MQ=255)
                                    aaaaaGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACAt                                                                                          >  3:349004/1‑139 (MQ=255)
                                    aaaaaGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACAt                                                                                          <  4:349004/139‑1 (MQ=255)
                                    aaaaaGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAACGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             >  6:58441/1‑104 (MQ=255)
                                    aaaaaGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAACGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  5:58441/104‑1 (MQ=255)
                                        aGCTGAATTTTATACAAGTAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGAc                                                                                      <  5:354725/139‑1 (MQ=255)
                                                  tataCAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCag                                                                            >  7:158016/1‑139 (MQ=255)
                                                         gaagCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGaa                                                                      >  4:328568/1‑138 (MQ=255)
                                                         gaagCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGaa                                                                      >  7:103428/1‑138 (MQ=255)
                                                                  aTGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGCCTTGGAATCAGGAGGAAAAGGTCtgtt                                                            >  2:118456/1‑139 (MQ=255)
                                                                     gtgtCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTgttgtt                                                         >  6:17955/1‑139 (MQ=255)
                                                                       gtCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             <  7:317912/69‑1 (MQ=255)
                                                                       gtCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAagccagc                                                                                                                             >  8:317912/1‑69 (MQ=255)
                                                                       gtCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTgttgtta                                                        >  2:263597/1‑137 (MQ=255)
                                                                                aCAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGaa                                              <  7:184406/139‑1 (MQ=255)
                                                                                aCAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTATGTTGTTTCAGCAAGGa                                               >  5:375415/1‑138 (MQ=255)
                                                                                    aCCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGt                                          >  7:136756/1‑139 (MQ=255)
                                                                                    aCCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGt                                          >  4:13798/1‑139 (MQ=255)
                                                                                    aCCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGt                                          >  2:63368/1‑139 (MQ=255)
                                                                                        ggAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCAtt                                      >  1:59194/1‑139 (MQ=255)
                                                                                         gAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCAtt                                      >  7:144131/1‑138 (MQ=255)
                                                                                         gAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCACCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCAtt                                      >  3:164998/1‑138 (MQ=255)
                                                                                         gAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCACCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCAtt                                      >  3:170671/1‑138 (MQ=255)
                                                                                            aCCGGAAAGCCGGTTGCAATTCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGTAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGaaa                                  <  3:374450/139‑1 (MQ=255)
                                                                                            aCCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGaaa                                  <  6:363077/139‑1 (MQ=255)
                                                                                            aCCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGCTTCAGCACGGAAACGTCATTGaaa                                  <  8:158016/139‑1 (MQ=255)
                                                                                                   aGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACt                            <  3:13798/138‑1 (MQ=255)
                                                                                                     ccGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGa                         <  5:17955/139‑1 (MQ=255)
                                                                                                     ccGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCACCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGa                         <  4:164998/139‑1 (MQ=255)
                                                                                                     ccGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCACCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGa                         <  4:170671/139‑1 (MQ=255)
                                                                                                             aatgcaaAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGAAACAAGTg                 <  8:144131/139‑1 (MQ=255)
                                                                                                             aatgcaaAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGAAACAAGTg                 <  1:63368/139‑1 (MQ=255)
                                                                                                               tgcaaAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGAAACAAGTGaa               <  3:328568/139‑1 (MQ=255)
                                                                                                                  aaaaTATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGAAACAAGTGAATgg            <  8:322190/139‑1 (MQ=255)
                                                                                                                         tttCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGAAACAAGTGAATGGCAGAAtt     >  6:77759/1‑139 (MQ=255)
                                                                                                                             aTTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGAAACAAGTGAATGGCAGAATtgtg  <  2:182733/138‑1 (MQ=255)
                                                                                                                                          |                                                                                                                            
ATAACGTTACAAATCTTGTCGAATACAACTATGATAAAAAAGCTGAATTTTATACAAGAAGCTCTGATGGTGTCATAACGACAGACCGGGAAACCGGAAAGCCGGTTGCAATGCAAAATATTTTCATTGTTGAAGCCAGCCATCATATTATTGATCAAGATGGAAGGCGGGACATCGACTTGGAATCAGGAGGAAAAGGTCTGTTGTTTCAGCACGGAAACGTCATTGAAACAGACTGGAAACAAGTGAATGGCAGAATTGTG  >  CP006881/679095‑679357

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