Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,158,080 G→T 12.1% R19L (CGG→CTG)  opuA → glycine betaine ABC transporter ATP‑binding protein/permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,158,0800GT12.1% 87.2 / 4.1 33R19L (CGG→CTG) opuAglycine betaine ABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base G (15/14);  new base T (2/2);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

AAAACCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGA  >  NC_002947/3157931‑3158221
                                                                                                                                                      |                                                                                                                                             
aaaaCCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGcgcg                                                                                                                                               >  1:608123/1‑150 (MQ=255)
   aCCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGAt                                                                                                                                            <  1:1137043/150‑1 (MQ=255)
     cccATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATcc                                                                                                                                          >  2:685987/1‑149 (MQ=255)
      ccATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGCCGCTCTTCACACTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCGGTTAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAg                                                                                                                                         >  1:969483/1‑150 (MQ=255)
       cATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAgc                                                                                                                                        <  2:1094554/150‑1 (MQ=255)
             gcgcTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATCCGCGCAgtg                                                                                                                                  <  2:440381/150‑3 (MQ=255)
             gcgcTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCTGATCCGCGCAgtg                                                                                                                                  <  1:685987/150‑3 (MQ=255)
                 ttGATTACGCGGGTTTTGGCCGATTGGCCGGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGatcgga                                                                                                                              >  1:37162/1‑145 (MQ=255)
                                  ggCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGAtt                                                                                                             >  1:1407050/1‑150 (MQ=255)
                                             aTGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCGTGTCCCTTTATCGATTCACTATATTgc                                                                                                  >  2:1374592/1‑150 (MQ=255)
                                                gAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCt                                                                                               <  1:1374592/150‑1 (MQ=255)
                                                          tCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCg                                                                                      <  1:620066/150‑1 (MQ=255)
                                                             cTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTTGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGc                                                                                   <  1:1139354/150‑1 (MQ=255)
                                                                  cTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCa                                                                              >  1:437714/1‑150 (MQ=255)
                                                                            ttGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTTCGGTCACCAACCGGTa                                                                    >  2:169041/1‑150 (MQ=255)
                                                                                              cGTGATCGCTGATAAATACCCACAGATAACCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGa                                                  >  1:398554/1‑150 (MQ=255)
                                                                                                   tCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTcca                                             <  2:608123/150‑1 (MQ=255)
                                                                                                        gATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAAc                                        >  2:484278/1‑150 (MQ=255)
                                                                                                                cccACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAACCGCCCTTGGCTCTGGATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGa                                <  1:169041/150‑1 (MQ=255)
                                                                                                                 ccACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAg                               <  2:1407050/150‑1 (MQ=255)
                                                                                                                  cacaGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTg                                                                                                                           >  1:1317353/1‑57 (MQ=255)
                                                                                                                  cacaGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTg                                                                                                                           <  2:1317353/57‑1 (MQ=255)
                                                                                                                           ccAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAAcc                     >  2:1049270/1‑150 (MQ=255)
                                                                                                                            cAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCgg                                                       >  1:857966/1‑115 (MQ=255)
                                                                                                                            cAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCgg                                                       <  2:857966/115‑1 (MQ=255)
                                                                                                                            cAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCt                    >  1:1239976/1‑150 (MQ=255)
                                                                                                                                          aaaCCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTg      <  1:1114533/150‑1 (MQ=255)
                                                                                                                                              gtcTGCGCTGATCCGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGa  >  1:996822/3‑150 (MQ=255)
                                                                                                                                                 tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTcgcg                <  2:429452/133‑1 (MQ=255)
                                                                                                                                                 tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTcgcg                <  2:429774/133‑1 (MQ=255)
                                                                                                                                                 tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTcgcg                >  1:429774/1‑133 (MQ=255)
                                                                                                                                                 tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTcgcg                >  1:429452/1‑133 (MQ=255)
                                                                                                                                                 tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGt              >  2:1037135/1‑135 (MQ=255)
                                                                                                                                                 tGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGt              <  1:1037135/135‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
AAAACCCCATCCCGCGCTTGATTACGCGGGTTTTGGCCGATTGGCATGGAGCCTGCTTCGCTGCCCTGGAATCCATTGAACGCCAGGGGGCGACGTGATCGCTGATAAATACCCACAGATTGCCAACGACCCTTTGGAAACCACTGCGCGGATCAGCGCAGACGCCCTTGGCTCTGTATCGATTCACGATATTGCGCTGGAATTTCGTGCGGTCACCAAACGGTACGAGGGGAGCCGGTCAGACTCCACCAACGCGGTTGAGCAAGTGAACCTCGCGGTGCGGCGTGGCGA  >  NC_002947/3157931‑3158221

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