Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 378,220 C→T 11.3% G400S (GGT→AGT)  gbcA ← glycine‑betaine dioxygenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947378,2200CT11.3% 97.7 / 2.6 35G400S (GGT→AGT) gbcAglycine‑betaine dioxygenase subunit
Reads supporting (aligned to +/- strand):  ref base C (15/16);  new base T (2/2);  total (17/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 = 1.00e+00

CTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTGCGCATGTTCTCGGGGTTGTAGT  >  NC_002947/378093‑378354
                                                                                                                               |                                                                                                                                      
ctgctgCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGt                                                                                                                                >  2:347089/1‑136 (MQ=255)
    tgcCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGt                                                                                                                            >  1:121759/1‑136 (MQ=255)
      ccGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCt                                                                                                                          >  1:88970/1‑136 (MQ=255)
          tgattGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGa                                                                                                                      <  2:258262/132‑1 (MQ=255)
              tGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATAc                                                                                                                  <  2:121759/136‑1 (MQ=255)
                 cGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACggg                                                                                                               <  2:170099/136‑1 (MQ=255)
                        aTCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACTAAACTCGTAAGTCTTCGAATACGGGCCCGGcg                                                                                                        >  1:46392/1‑135 (MQ=255)
                        aTCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACTAAACTCGTAAGTCTTCGAATACGGGCCCGGCt                                                                                                        >  1:43132/1‑136 (MQ=255)
                           tgTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGt                                                                                                     <  1:483188/135‑1 (MQ=255)
                                  actGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGt                                                                                              >  1:466997/1‑136 (MQ=255)
                                  actGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGt                                                                                              >  1:272371/1‑136 (MQ=255)
                                         cTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTg                                                                                       >  1:622207/1‑136 (MQ=255)
                                         cTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTg                                                                                       >  1:339335/1‑136 (MQ=255)
                                              tttcttTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATccc                                                                                  <  1:539305/135‑1 (MQ=255)
                                              cttcttTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATccc                                                                                  >  1:187145/1‑136 (MQ=255)
                                              cttcttTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATccc                                                                                  >  2:616997/1‑136 (MQ=255)
                                                 cttTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCATGGAGTTGATCCCACg                                                                               >  2:440956/1‑136 (MQ=255)
                                                         tATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTtct                                                                       <  2:331496/136‑1 (MQ=255)
                                                            ggCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCtt                                                                    >  1:376620/1‑136 (MQ=255)
                                                              cgcCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGGACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCg                                                                  >  2:338548/1‑136 (MQ=255)
                                                               gcCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCCCGCTGGTTCTCTTCgg                                                                 >  1:563385/1‑136 (MQ=255)
                                                                                      ccccACCCGCCCGATGTACCAGTCAATTAAGTAAACCACACCAAACTCGTAAGTCTTCGAATACAGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTtg                                          <  1:93306/133‑1 (MQ=255)
                                                                                      cAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTtg                                          <  1:279179/136‑1 (MQ=255)
                                                                                      cAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTtg                                          <  2:466997/136‑1 (MQ=255)
                                                                                      cAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCATGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTtg                                          <  1:440956/136‑1 (MQ=255)
                                                                                        gCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCTGCCAGACGGCGGTCCTGGTCGTtggt                                        <  1:173509/136‑1 (MQ=255)
                                                                                                ccGCTGTACCAGTCAATGAAGTTCACCACACTAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTccc                                <  2:43132/136‑1 (MQ=255)
                                                                                                ccGCTGTACCAGTCAATGAAGTTCACCACACTAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTccc                                <  2:46392/136‑1 (MQ=255)
                                                                                                 cGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTcccc                               >  2:488609/1‑135 (MQ=255)
                                                                                                   cTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCcaca                             <  1:422207/136‑1 (MQ=255)
                                                                                                         cAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTTCGCCGTGGAGTTGATCCCACGCTGGTTCTCTGCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTgc                       <  2:512197/136‑1 (MQ=255)
                                                                                                         cAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTgc                       <  2:161201/136‑1 (MQ=255)
                                                                                                         cAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTgc                       <  2:578958/136‑1 (MQ=255)
                                                                                                              aaTGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTGCGCATg                  >  2:324974/1‑136 (MQ=255)
                                                                                                               aTGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTGCGCATGt                 <  2:106788/136‑1 (MQ=255)
                                                                                                                             accaAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTGCGCATGTTCTCGGGGTTGTAg   <  2:318063/136‑1 (MQ=255)
                                                                                                                              ccaAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTGCGCATGTTCTCGGGGTTGTAGt  <  1:587337/136‑1 (MQ=255)
                                                                                                                               |                                                                                                                                      
CTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGCCCGGCTGGTACGCCGTGGAGTTGATCCCACGCTGGTTCTCTTCGGCCAGACGGCGGTCCTGGTCGTTGGTGGCGTCCCACACCTTGCGCATGTTCTCGGGGTTGTAGT  >  NC_002947/378093‑378354

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