Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 341,380 C→T 14.8% intergenic (‑271/‑107) PP_5438 ← / → gabP‑I hypothetical protein/gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947341,3800CT14.8% 66.4 / 4.5 27intergenic (‑271/‑107)PP_5438/gabP‑Ihypothetical protein/gamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base T (2/2);  total (11/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 = 8.91e-01

TTTGTAAGTGGGTATTTCAGTAATCCGCTGATTTGAGGTCCTTTCGGTCATGCAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCGGTAACAATTCCAATGACCTCGC  >  NC_002947/341252‑341515
                                                                                                                                |                                                                                                                                       
tttGTAAGTGGGTATTTCAGTAATCCGCTGATTTGAGGTCCTTTCGGGCATGCAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTAtt                                                                                                                                  >  2:316435/1‑136 (MQ=255)
   gTAAGTGGGTATTTCAGTAATCCGCTGATTTGAGGTCCTTTCGGTCATGCAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGcc                                                                                                                               >  1:150303/1‑136 (MQ=255)
   gTAAGTGGGTATTTCAGTAATCCGCTGATTTGAGGTCCTTTCGGTCATGCAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGcc                                                                                                                               >  2:486990/1‑136 (MQ=255)
                       tCCGCTGATTTGAGGTCCTTTCGGTCATGCAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGc                                                                                                           <  2:443838/136‑1 (MQ=255)
                       tCCGCTGATTTGAGGTCCTTTCGGTCATGCAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGACGCGCTTATTGCCAGTTGGAAACCGCCGATTGc                                                                                                           <  2:364431/136‑1 (MQ=255)
                       cccGCTGTTTTGTGGCCCTTGCGGTCTTGCAAGCGCCTAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGc                                                                                                           <  2:644625/135‑1 (MQ=255)
                              ttttGTGGTCCTTTCGGAATTGCATGGGCCTAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACg                                                                                                    <  1:68497/135‑1 (MQ=255)
                                      tCCTTTCGGTCATGCAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTc                                                                                            >  1:584724/1‑136 (MQ=255)
                                        ttttCGGCCTTCCAATGGCCTAACGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCAt                                                                                          <  2:646402/135‑1 (MQ=255)
                                                   ccAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGgcgcg                                                                               <  2:342113/135‑1 (MQ=255)
                                                     aaTGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGct                                                                             >  2:417056/1‑136 (MQ=255)
                                                            gAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTAcc                                                                      >  1:345557/1‑136 (MQ=255)
                                                                   cTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTc                                                               <  2:584492/136‑1 (MQ=255)
                                                                   cTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTc                                                               <  2:150303/136‑1 (MQ=255)
                                                                    tCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTTATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCg                                                              <  2:74333/136‑1 (MQ=255)
                                                                      cAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGTGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGc                                                            >  2:501651/1‑136 (MQ=255)
                                                                      cAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGTGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGc                                                            >  2:504806/1‑136 (MQ=255)
                                                                           cGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGTGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGt                                                       <  1:504806/136‑1 (MQ=255)
                                                                           cGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGTGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGt                                                       <  1:501651/136‑1 (MQ=255)
                                                                              aTCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTTGAAACCGCCGATTGCAGCGACGCCTAGGTCAt                                                                                          >  1:577055/1‑98 (MQ=255)
                                                                              aTCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTTGAAACCGCCGATTGCAGCGACGCCTAGGTCAt                                                                                          <  2:577055/98‑1 (MQ=255)
                                                                                        gCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACaa                                          <  1:417056/136‑1 (MQ=255)
                                                                                               gAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGACACAACTTACAAAAGTAAtt                                  >  2:455308/1‑136 (MQ=255)
                                                                                                          cGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCg                        >  2:65492/1‑136 (MQ=255)
                                                                                                              tGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCGGTaa                    <  2:620537/136‑1 (MQ=255)
                                                                                                                        atatGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCGGTAACAATTCCAAt          <  2:495234/136‑1 (MQ=255)
                                                                                                                         tatGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCGGTAACAATTCCAATg         <  1:316435/136‑1 (MQ=255)
                                                                                                                                cgcTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCGGTAACAATTCCAATGACCTCGc  <  2:145047/136‑1 (MQ=255)
                                                                                                                                cgcTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCGGTAACAATTCCAATGACCTCGc  <  2:584724/136‑1 (MQ=255)
                                                                                                                                cgcTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCGGTAACAATTCCAATGACCTCGc  <  1:390489/136‑1 (MQ=255)
                                                                                                                                |                                                                                                                                       
TTTGTAAGTGGGTATTTCAGTAATCCGCTGATTTGAGGTCCTTTCGGTCATGCAATGGCCGAAAGCCCTCCAGGCCGCATCCTAAAAGGCCTGCAGAATTTTCTGACGAGTGCATTCGTCATATGGCGCGCTTATTGCCAGTTGGAAACCGCCGATTGCAGCGACGCCTAGGTCATCCCCTGGCGCGCTCTGTACCCGCAGTCGGCGTGGTCCACTCCTTACAAAGGTAGTTTCAATGAGCGGTAACAATTCCAATGACCTCGC  >  NC_002947/341252‑341515

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: