Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,867,552 C→A 12.5% R81L (CGG→CTG)  fbpC ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,867,5520CA12.5% 101.2 / 2.3 32R81L (CGG→CTG) fbpCABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (16/12);  new base A (2/2);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

TACCAGTTCGAGCATTTCGGCAACCCGCTGGCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCCACAGCCGGATGGGCCGAGCAGGCACACCA  >  NC_002947/5867419‑5867684
                                                                                                                                     |                                                                                                                                    
tACCAGTTCGAGCATTTCGGCAACCCGCTGGCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCcgc                                                                                                                                    >  2:253091/1‑136 (MQ=255)
tACCAGTTCGAGCATTTCGGCAACCCGCTGGCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCcgc                                                                                                                                    >  2:112859/1‑136 (MQ=255)
  ccAGTTCGAGCATTTCGGCAACCCGCTGGCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCgcgc                                                                                                                                  >  2:473285/1‑136 (MQ=255)
                           cTGGCGGGCCTGCTCGCGACCACTGCCCGTCAGG‑CATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATcgc                                                                                                        <  1:486948/136‑1 (MQ=255)
                            tGGCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATcgc                                                                                                        <  1:473316/136‑1 (MQ=255)
                            tGGCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCAGCGCCTGGGGGGACAGATCGGAGATATcgc                                                                                                        >  2:258332/1‑136 (MQ=255)
                            tGGCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCAGCGCCTGGGGGGACAGATCGGAGATATcgc                                                                                                        >  2:231894/1‑136 (MQ=255)
                            gggCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATcgc                                                                                                        <  2:259270/135‑1 (MQ=255)
                                 ggCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCg                                                                                                   <  1:112859/136‑1 (MQ=255)
                                      gCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGAt                                                                                              <  2:403424/136‑1 (MQ=255)
                                                aCTGCCCGTCAGTCCATAGGCGTTATTGGCTTCGTCGGTAAGGTTGTGAAACAGAGCGTAGGACTGTAACAGGAGGCCGTAGTCACGCGCGTGGGGGGACAGTTCGGTGATTTCGCGTTCGCCGATGTACAGCGTa                                                                                    <  2:56848/136‑1 (MQ=255)
                                                    cccGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACcgc                                                                                >  2:433283/1‑136 (MQ=255)
                                                               aTAGGCGATATTTGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGgcgc                                                                     >  2:31188/1‑136 (MQ=255)
                                                                 aGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTc                                                                   <  1:447578/136‑1 (MQ=255)
                                                                 aGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTc                                                                   <  1:180061/136‑1 (MQ=255)
                                                                 aGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTc                                                                   >  2:37618/1‑136 (MQ=255)
                                                                         ttGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCgt                                                           <  1:253091/136‑2 (MQ=255)
                                                                         ttGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCgg                                                           >  1:331486/1‑136 (MQ=255)
                                                                              ttCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATg                                                      <  1:274067/136‑1 (MQ=255)
                                                                              ttCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCAGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATg                                                      <  1:258332/136‑1 (MQ=255)
                                                                              ttCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCAGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATg                                                      <  1:231894/136‑1 (MQ=255)
                                                                                    ggTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACa                                                                                                                       >  2:192136/1‑65 (MQ=255)
                                                                                    ggTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACa                                                                                                                       <  1:192136/65‑1 (MQ=255)
                                                                                    ggTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACgc                                                >  2:43557/1‑136 (MQ=255)
                                                                                    ggTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACgc                                                >  2:41239/1‑136 (MQ=255)
                                                                                                aaaCAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTc                                    >  2:246626/1‑136 (MQ=255)
                                                                                                aaaCAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTc                                    >  1:171215/1‑136 (MQ=255)
                                                                                                        cgTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCcaca                            >  2:493582/1‑136 (MQ=255)
                                                                                                                tGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCCACAGCCGGATg                    <  2:383218/136‑1 (MQ=255)
                                                                                                                     cAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCCACAGCCGGATGGGCCg               >  2:355878/1‑136 (MQ=255)
                                                                                                                     cAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCCACAGCCGGATGGGCCg               >  2:432489/1‑136 (MQ=255)
                                                                                                                     cAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCCACAGCCGGATGGGCCg               >  2:23918/1‑136 (MQ=255)
                                                                                                                          tGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCCACAGCCGGATGGGCCGAGCAg          >  2:331819/1‑136 (MQ=255)
                                                                                                                                  tCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCCACAGCCGGATGGGCCGAGCAGGcagagca  >  1:306968/1‑131 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TACCAGTTCGAGCATTTCGGCAACCCGCTGGCGGGCCTGCTCGCGACCACTGCCCGTCAGGCCATAGGCGATATTGGCTTCGACGGTAAGGTTGGGAAACAGCGCGTAGGACTGGAACAGGATGCCGTAGTCCCGCGCCTGGGGGGACAGATCGGAGATATCGCGTTCGCCGATGTACAGCGTACCGCGATCCTGGCGCTCCAGGCCGGCAATGCAACGCAGCAGCGTGGTCTTGCCACAGCCGGATGGGCCGAGCAGGCACACCA  >  NC_002947/5867419‑5867684

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: ‑