Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,604,223 C→T G59G (GGC→GGT)  lsrD → Autoinducer‑2 ABC transporter membrane subunit LsrD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2230CT95.5% 48.9 / ‑3.4 22G59G (GGC→GGT) lsrDAutoinducer‑2 ABC transporter membrane subunit LsrD
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (6/15);  total (6/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 = 7.04e-01

CGCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCC  >  NC_000913/1603992‑1604457
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
tgCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCGACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGttg                                                                                                                                                                                                                          <  1:136002/249‑2 (MQ=255)
              gTTAACCCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCGCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGGGACTTTATCTGCATTGGCATGGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCggn                                                                                                                                                                                                           <  1:97663/251‑2 (MQ=255)
                                                                         tGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCACAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTg                                                                                                                                                                                                                         <  1:1375/178‑1 (MQ=255)
                                                                           ggAACNNGCNNNTGNNGCACTGCTCGTTATTGAGATTGTNGCATTTGGTGCAATTAACCCGCGAATGTTNGNTNTCAANNTGTNGCTGTTCAGCNCCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGNGNTTGTAAGTGGTGGGATCGATATTTCGTTTGGNTCGACCATCGGCCTCTGCGCCATTNCATTnnggntnnnnnnn                                                                                                                                                                             >  2:367377/1‑210 (MQ=255)
                                                                                          ngCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCg                                                                                                                                                                                                                    >  1:98599/2‑166 (MQ=255)
                                                                                                      tATTGNNATNNNCGNNTTTGGTGCAATTAACCCGCGAATNTTAGATCTCAATATGTTGCTGTTCAGCACNANTNCCTTNNTCTNCATTGGCATTNTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTNGNTTGGTTCGACCATCGGCCTCTGCGCCATTGCNTTggg                                                                                                                                                                                       >  2:379140/1‑183 (MQ=255)
                                                                                                      tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATGTCTTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTggg                                                                                                                                                                                       <  1:379140/183‑1 (MQ=255)
                                                                                                          gagaTTGTCGCATTTCGTGCAATTACCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCCGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGCTGGAGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCgg                                                                                                               >  1:158553/1‑251 (MQ=255)
                                                                                                               tGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATg                                                                                                                                                           <  1:238052/202‑1 (MQ=255)
                                                                                                                             cAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCTCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGc                                                                                                                               <  1:333960/216‑1 (MQ=255)
                                                                                                                                        cgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGc                                                                                                                                                   >  1:32825/1‑185 (MQ=255)
                                                                                                                                                                        cACCAGTGACTTTATCTGCAGTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTcn                                                                                                                                                                                                                             <  1:103392/79‑2 (MQ=255)
                                                                                                                                                                                aCTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGGTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTTCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGTATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTAt                                         <  1:487075/251‑1 (MQ=255)
                                                                                                                                                                                         gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGTGCGTACTGTTTCAAAGTGGGGTGCCGATGCCGCGGTCGATACTCCGGACCTTACTGCGCGGCGCATTGTGTGGTGTGATCAACGCCGGATTAATTATCAATACCAAAGTGATCCCGCTGGTGATTACGCTTGGCACGctgtttctgttt                                   <  1:367377/248‑1 (MQ=255)
                                                                                                                                                                                         gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCAGCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGc                                                                                                                                                   <  1:171612/136‑1 (MQ=255)
                                                                                                                                                                                                cattgTCGCCCTACCGCGAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGGTTGTTTCGAGCATCGGCCTCTTCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACt                                                                                                                                 <  1:371385/147‑1 (MQ=255)
                                                                                                                                                                                                cattgNNGCNNNACNNCTAACGATGGTGATTGTCAGTGGNGGGATCGATATTTCGTTTGGTTCGACCATNGNCNTCTGNNCCANTGCATTGGGCNTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTNANt                                                                                                                                 >  2:371385/1‑147 (MQ=255)
                                                                                                                                                                                                   tgTCGCCCTACCGCTAACGATGGTGATTGTCATTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCGAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAATTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGct                       <  1:505283/250‑1 (MQ=255)
                                                                                                                                                                                                       gCCCTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCGGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCgg                              <  1:73161/239‑1 (MQ=255)
                                                                                                                                                                                                        cccTACCGCTAACGATGGTGATTTTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCg                                                                                                                                                                                                             <  1:502348/63‑1 (MQ=255)
                                                                                                                                                                                                        cccTACCGCTAACGATGGTGATTGTCAGTGGTGGGATCGATCTTTCGTCTGGTTCGCCCATCGGACTCTGCGCAATTGCATTGGGCGTACTGTTTCCACGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGgcgc                                                                                                                        >  1:462832/1‑148 (MQ=255)
                                                                                                                                                                                                                       gggTGATTGTCAGTGGTGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAATTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCTCTCGGCTGCTTTCCGGTATGGcc  <  1:419800/250‑1 (MQ=255)
                                                                                                                                                                                                                              tNNCAGTGGNNGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGc                                                                                                                                                   <  2:32825/99‑1 (MQ=255)
                                                                                                                                                                                                                               gTCAGTGGTGGGATCGATATTTCGTTTGGTGCGACCATCGGCCCCTGCGCCATTGCATGGGGCTTACTGTTTCAAAGTGGTGTGCCTATGCCGCTGGCGATACTCCTGACCt                                                                                                                                     <  1:231755/112‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
CGCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCC  >  NC_000913/1603992‑1604457

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