Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,033,608 T→C 16.1% L217P (CTC→CCC)  yedY → membrane‑anchored, periplasmic TMAO, DMSO reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,033,6080TC16.1% 99.8 / 2.3 53L217P (CTC→CCC) yedYmembrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base T (19/24);  new base C (8/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 7.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACC  >  CP009273/2033461‑2033756
                                                                                                                                                   |                                                                                                                                                    
ccGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTc                                                                                                                                                     <  2:83923/149‑1 (MQ=255)
        aCCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTATGGCGCAGGGGTTGGTCTCGACGAAACAAAGCATCCGCCCACACtgat                                                                                                                                             >  1:612430/1‑149 (MQ=255)
            gCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATTCTTATGTTGAAGGAGTGCGTGTCGCCGGAGCAAAGCAACCGCTCACACTGATGACCg                                                                                                                                        >  2:529236/1‑150 (MQ=255)
            gCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGCGAAAAAAACCTATGTGGAAGGATTGCGTCTCGCCGAAGCAAATTAACCGATCAAACTGATGAACg                                                                                                                                        >  2:313378/1‑150 (MQ=255)
                tGCGAGGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAgg                                                                                                                                    <  1:475668/150‑1 (MQ=255)
                   gAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTATGTCGAAGGATTGCGGCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGGAggggg                                                                                                                                 >  1:364861/1‑147 (MQ=255)
                       tATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTAt                                                                                                                             <  2:344214/150‑1 (MQ=255)
                         tGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATgg                                                                                                                           <  2:299203/150‑1 (MQ=255)
                                 tCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGc                                                                                                                     <  1:812630/148‑1 (MQ=255)
                                  cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTTTTGAAGGGTTGCGTCTCGACGGAGCAAAGCATCCCCCCCCACTGATGACCGTCGGGGGTTTTTGGAAGGCggt                                                                                                                  >  2:521971/1‑148 (MQ=255)
                                  cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGCGAAAAATTCTTTTGTTTTAGGATTACGCCTCTGCGCCGCAAAGCACCCGCCCACCCTGATGACAGTAGGGGGTTATGGCAAGGCGt                                                                                                                   >  2:453524/1‑149 (MQ=255)
                                   gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTCGAAGGATTTCGGCTCGACGAAACAAAGCATCCGCTCCCACTGATGACCGTAGGGGGTTTTGGCAAGGCGGTa                                                                                                                 >  2:424589/1‑150 (MQ=255)
                                   gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGtt                                                                                                                  <  2:594150/149‑1 (MQ=255)
                                        aaTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTTTCGAAGGATTGGGTCTCGACGAAGCAATGCATCCGCCCCCACTGATGACCGTAGGGGGTTATGGGAAGGCGGTACCGcc                                                                                                            >  2:464120/1‑150 (MQ=255)
                                         aTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCca                                                                                                           <  1:1062546/150‑1 (MQ=255)
                                             aTGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTGTCGAAGGATTTCGGCTCGGCGGAGCAATGCATCCGCTCACACTGATGGCCGTAGGGGGTTTTTGGAAGGCGTTACCGCCCCaa                                                                                                        >  1:558139/1‑149 (MQ=255)
                                                       aCAGATGCCAGGCCAGCAGGACCGTTTTATCGGCGGTGGGCTGAAATATCTTTATGTGGAGGGATTGCGTCTCGACGAAGAAATGCACCCGCTCACACTGATGACCGTGGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGcc                                                                                             <  2:721065/150‑1 (MQ=255)
                                                       aCAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATTCCTTTGTCCAAGGAGTGCGGCTCGACGCAGCAAAGCATCCCCCCACACTGATGACCGCAGGGGTTTATTGCAAAGCGGTATCCCCACAAAATTGCGCGcc                                                                                             >  1:351533/1‑150 (MQ=255)
                                                         aGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGCAAAATCCTTTTATCTAAGGAAATCGCCCCGCCGCAGCAATTCATCCGCCCACACCGCTGACCGTAGGGGGGTATTGGCAGGGGCTATCGCCCCAAAAAAGCGCGcccg                                                                                           >  2:30992/1‑148 (MQ=255)
                                                           aTGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTGTTGCAGGATTTCGTCTCGCCGAAGCAAAGCATCCCCTCACACTGATTACCGGAGGGGTTTTTTGCAAGGCGTTACCCCCACCAAAAAGCGCGCCgggg                                                                                         >  1:361091/1‑148 (MQ=255)
                                                                 ggccagCAGGACCGCTTTATCGGCGGCGGGCCGAAATATCCCTCTGTCGAAGGAATGCGGCTCCACGAAACAATGCATACGCCCACACTCATGAGCGCAGGTGGTTATGGCAAGGCGGTAACCCCCCAAAAAAGCGGGGCCGGGCGACTg                                                                                   >  1:312614/1‑150 (MQ=255)
                                                                    cagcggGACCGCTTTATCGGCGGTGGGTTAAAATATCTTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTGGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGAtt                                                                                <  1:521971/150‑1 (MQ=255)
                                                                    cagcagGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTTGAAGGATTTCGGCTCCACGAAACAATTCATCCGCCCACACTGATGACCGTAGGTGTTTTTTGCAAGGGGGTACCGCCCCCAAAAGGGGCGCCCGGGGGACTGGtt                                                                                >  2:61445/1‑150 (MQ=255)
                                                                     agcagGACCGTTTTATCGGCGGTGGGCTAAAATATCTTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACCCTGATGCCCGTGGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTACGACTGATtg                                                                               <  2:364861/150‑1 (MQ=255)
                                                                       cagGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCCCTCCCACTGATGACCGGAGGTGGTTATGGCAAGGCGTTACCGCCACAAAAAGGCGCGCCGGGGCGGCTGATtgtg                                                                             >  2:1015900/1‑150 (MQ=255)
                                                                         ggACGGCTTTATCGGCGGTGGGCTAAAATTTCCTTATGTCGAAGGATTGCTTCTCGACGAAGCATTGCACCCCCTCACTCTAATGACGGTGGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGcc                                                                           <  1:763779/150‑1 (MQ=255)
                                                                           aCCGCTTTATCGGCGGTGGGCTGAAATTTCCTTATGTCGAAGGATTGCGTCTCGACGAAGAAATGCATCCGCTCCCACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGt                                                                         <  1:636034/150‑1 (MQ=255)
                                                                             ggCTGTATCGGGGGTGGGCTGAAATATCTTTAGGTCGAAGGATTGCGTCTGGAGGAAGAAATGCATCCGCTCACACGGATGGCCGTGGGTGTTTAGGGGAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTgg                                                                       <  2:119028/149‑1 (MQ=255)
                                                                                    ccggcggTGGGCTGAATTTTCCTTATGGCGAAGGATGTCGTCCCGACAAGCCATTGCTCCCGCCCCCACTGATAACCGTTGTTTTTGTTGGCATGGAATTGCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAatat                                                                 <  2:496340/148‑1 (MQ=255)
                                                                                          gTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTa                                                          <  1:1015900/150‑1 (MQ=255)
                                                                                          gTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCACCGGTGCGACTGATTGTGCCGTGGAAATATGGCttt                                                           >  1:371351/1‑149 (MQ=255)
                                                                                           tGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTa                                                          >  2:413379/1‑149 (MQ=255)
                                                                                              gCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAgg                                                      >  2:27634/1‑150 (MQ=255)
                                                                                               cTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAgg                                                      >  1:1096963/1‑149 (MQ=255)
                                                                                                gaaaaTATCCTTATGTGGAGGGATGGCGTCTCGACGAAGAAATGAACCCGCTCACCCGGATGACCGTGGGTGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGa                                                    <  2:1079527/148‑1 (MQ=255)
                                                                                                     tttCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTaaa                                               <  1:424589/148‑1 (MQ=255)
                                                                                                       tCCTTATGTCGAGGGATTGCGTCTCGACGAAGAAATGCATCCGCTCCCAGTGATAGCGGTAGGTGTTTATGACAGGGGGTTACCGCCACAATATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAAt                                              <  1:313378/149‑1 (MQ=255)
                                                                                                          ttATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGa                                           >  1:754959/1‑149 (MQ=255)
                                                                                                            aTGTCGAAGGTTTGCGTCTCGAGAAAGAAATGCACCCGCTCACACTGATGGCGGTGGGTTTTTATGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATa                                         <  2:351533/149‑1 (MQ=255)
                                                                                                                 gAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGt                                   <  1:1058780/150‑1 (MQ=255)
                                                                                                                 gAAGGATTGCGTCGCAACGAAGCAATGCATCCGCTCACACGGGTGACCGTAGTTGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAg                                    <  2:865124/149‑1 (MQ=255)
                                                                                                                   aGGATTGCGTCTCGACGAAGCAATGCATCCGCTCCCCCTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTAt                                 <  1:839678/150‑1 (MQ=255)
                                                                                                                     gATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTAtt                                >  1:625893/1‑149 (MQ=255)
                                                                                                                         gCGTCTCGACGAAGCAATGCATCCGCCCACACTGATGACGGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGc                            <  2:381933/149‑1 (MQ=255)
                                                                                                                               cgacgaAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCcg                     >  2:300849/1‑150 (MQ=255)
                                                                                                                                gacgaAGCAATGCATCCGCTCCCACTGATGACCGTAGGTGTTTATGGCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCcgc                    <  2:139627/150‑1 (MQ=255)
                                                                                                                                 acgaAGCAATGCACCCGCTCACACTGATGACCGTGGGTGTTTAGGGAAAGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCcgcg                   <  2:361091/150‑1 (MQ=255)
                                                                                                                                       cAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTc             >  2:824978/1‑150 (MQ=255)
                                                                                                                                        aaTGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGCTTGTGCCGTGGAAATATGGCTTTACAGGGATTAAATCGACAGTCAGTATTAAGCTGACCCGCGAGCGAcc            >  2:451610/1‑150 (MQ=255)
                                                                                                                                             aTCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCaa       >  1:679896/1‑150 (MQ=255)
                                                                                                                                             aTCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCaa       >  1:150432/1‑150 (MQ=255)
                                                                                                                                              tCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAac      >  2:64399/1‑150 (MQ=255)
                                                                                                                                                  cTCCCACTAATGACCGTGGGTGTTTATGGCAAGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAaccacc  <  1:613731/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACC  >  CP009273/2033461‑2033756

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