Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,265,482 C→G 6.7% G169R (GGC→CGC)  U712_16795 ← putative membrane protein yvaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,265,4820CG 6.7% 200.7 / 6.0 60G169R (GGC→CGC) U712_16795putative membrane protein yvaC
Reads supporting (aligned to +/- strand):  ref base C (35/21);  new base G (2/2);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CTAAGCATGAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAAAGGCAAGGACAAAAAACATCGGCGCAGGGCCTTG  >  CP006881/3265349‑3265617
                                                                                                                                     |                                                                                                                                       
cTAAGCATGAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGAtt                                                                                                                                     <  7:271125/138‑1 (MQ=255)
   aGCATGAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAg                                                                                                                                 >  6:76663/1‑139 (MQ=255)
        gAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTa                                                                                                                                    <  4:172744/131‑1 (MQ=255)
        gAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTa                                                                                                                                    >  3:172744/1‑131 (MQ=255)
        gAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATa                                                                                                                            >  8:97185/1‑139 (MQ=255)
        gAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATa                                                                                                                            >  4:182921/1‑139 (MQ=255)
        gAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATa                                                                                                                            >  3:17583/1‑139 (MQ=255)
        gAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGGTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCGGATTAAAGCGATa                                                                                                                            >  4:197789/1‑139 (MQ=255)
        gAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGATCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCGGATTAAAGCGATa                                                                                                                            >  4:201761/1‑139 (MQ=255)
             ccGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCAt                                                                                                                       >  5:8220/1‑139 (MQ=255)
             ccGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCAt                                                                                                                       <  8:72553/139‑1 (MQ=255)
                tCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGTCTGAAAGACATGCAGCCAGCTGAGGATATGTATTTTGCAGCACCGCTCCCTCCTGTCCATGTTTGTTCCAGAACCAGCCGATTAACGCGATACCCATGGc                                                                                                                    >  8:248221/1‑139 (MQ=255)
                tCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGc                                                                                                                    >  6:279794/1‑139 (MQ=255)
                tCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGc                                                                                                                    >  4:213379/1‑139 (MQ=255)
                          tgCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATACCCAGCCGATTAAAGCGCTACCCCTGGCGAAAATgccg                                                                                                          >  7:233885/1‑139 (MQ=255)
                          tgCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGCAAATgccg                                                                                                          >  7:258042/1‑139 (MQ=255)
                          tgCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATgccg                                                                                                          >  2:36364/1‑139 (MQ=255)
                          tgCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATgccg                                                                                                          >  1:113343/1‑139 (MQ=255)
                          tgCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATgccg                                                                                                          >  3:66807/1‑139 (MQ=255)
                          tgCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATgccg                                                                                                          >  3:233646/1‑139 (MQ=255)
                          tgCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATgccg                                                                                                          >  7:34765/1‑139 (MQ=255)
                          tgCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATgccg                                                                                                          >  8:46774/1‑139 (MQ=255)
                                      aCCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTa                                                                                                                                    >  7:234587/1‑101 (MQ=255)
                                      aCCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTa                                                                                                                                    <  8:234587/101‑1 (MQ=255)
                                      aCCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGt                                                                                               >  2:4652/1‑138 (MQ=255)
                                           gcgcTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCGGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTc                                                                                          <  3:201761/138‑1 (MQ=255)
                                           gcgcTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCGGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTc                                                                                          <  3:197789/138‑1 (MQ=255)
                                           gCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTc                                                                                          <  3:182921/138‑1 (MQ=255)
                                           gCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTc                                                                                          <  7:97185/138‑1 (MQ=255)
                                              ttCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGc                                                                                      <  5:76663/139‑1 (MQ=255)
                                                ccGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGct                                                                                     <  1:36364/138‑1 (MQ=255)
                                                ccGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGct                                                                                     <  4:66807/138‑1 (MQ=255)
                                                     ggCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAAc                                                                               <  4:17583/139‑1 (MQ=255)
                                                           aaGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAgcgc                                                                          <  6:8220/138‑1 (MQ=255)
                                                                 tGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAACGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGACCTGCGCTtgtg                                                                     >  2:110228/1‑137 (MQ=255)
                                                                 tGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTtgtg                                                                     >  2:168654/1‑137 (MQ=255)
                                                                        agcTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTc                                                                                          >  4:44677/1‑109 (MQ=255)
                                                                        agcTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTc                                                                                          <  3:44677/109‑1 (MQ=255)
                                                                          cTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATcggcc                                                          >  3:103841/1‑139 (MQ=255)
                                                                          cTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATcggcc                                                          >  3:99570/1‑139 (MQ=255)
                                                                             atgatATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATcggccgg                                                        <  4:99570/138‑1 (MQ=255)
                                                                             atgatATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATcggccgg                                                        <  3:107354/138‑1 (MQ=255)
                                                                             atgatATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATcggccgg                                                        <  7:46774/138‑1 (MQ=255)
                                                                             atgatATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATcggccgg                                                        <  4:103841/138‑1 (MQ=255)
                                                                                atatGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGc                                                    >  6:245371/1‑139 (MQ=255)
                                                                                  atGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGAAACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCACCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGccc                                                  >  5:167595/1‑139 (MQ=255)
                                                                                       ttttGCAGCACCGCTCCCTCCGGTCCATGTTTTTTCCATAACCAGCCGATTAAATCGATACACATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAAc                                             >  7:70973/1‑139 (MQ=255)
                                                                                       ttttGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAg                                                                             <  4:112029/107‑1 (MQ=255)
                                                                                       ttttGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAg                                                                             >  3:112029/1‑107 (MQ=255)
                                                                                       ttttGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGaa                                              <  4:233646/138‑1 (MQ=255)
                                                                                          tGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTg                                           <  2:280408/138‑1 (MQ=255)
                                                                                                              tCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAAAGGCAAGGACaaa                       >  3:15126/1‑138 (MQ=255)
                                                                                                              tCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCAGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAATCCCAAGGACaaa                       >  7:294727/1‑138 (MQ=255)
                                                                                                                       gttCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAAAGGCAAGGACAAAAAACATCGgc             >  3:64931/1‑139 (MQ=255)
                                                                                                                        ttCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAAAGGCAAGGACAAAAAACATCGgcg            >  2:266845/1‑139 (MQ=255)
                                                                                                                          ccATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAAAGGCAAGGACAAAAAACATCGGCGCa          <  2:113343/139‑1 (MQ=255)
                                                                                                                             tAACCAGCCGATTAAAGCGATACTCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAAAGGCAAGGACAAAAAACATCGGCGCAggg       >  7:70993/1‑139 (MQ=255)
                                                                                                                             tAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTc                                                                                          <  2:276185/56‑1 (MQ=255)
                                                                                                                             tAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTc                                                                                          >  1:276185/1‑56 (MQ=255)
                                                                                                                                   gCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAAAGGCAAGGACAAACAACATCGGCGCCGGGCCTTg  >  6:170515/1‑138 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CTAAGCATGAGCACCGTCTGATGCTGTGCGTTATGAAAACCTGGCGTTCCGACGGCTGAAAGACATGCAGCCAGCTGATGATATGTTTTTTGCAGCACCGCTCCCTCCGGTCCATGTTTGTTCCATAACCAGCCGATTAAAGCGATACCCATGGCGAAAATGCCGCCGAGAAACGTCAGTCCCGCTCTGAACAGCGCTTGTGAAGGATCGGCCGGCAGCCCGGAACTGACTAAAAAGGCAAGGACAAAAAACATCGGCGCAGGGCCTTG  >  CP006881/3265349‑3265617

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