Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,814,461 C→T L91L (CTG→CTA)  luxS ← S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,4610CT95.7% 74.8 / NA 24L91L (CTG→CTA) luxSS‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/13);  minor base G (0/1);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTT  >  NC_000913/2814222‑2814690
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     
aTGTTCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATaa                                                                                                                                                                                                                           <  1:283070/251‑1 (MQ=255)
  gtgCAGTTCCTGCAACTTCTCTTTCGGCAGTACCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGCATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATaaa                                                                                                                                                                                                                          >  1:360181/1‑250 (MQ=255)
          ccTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGc                                                                                                                                                                                                                 <  2:171592/251‑1 (MQ=255)
                        ttCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCAtt                                                                                                                                                                                                      <  1:253639/248‑1 (MQ=255)
                                ggCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATc                                                                                                                                                                                           <  1:273532/250‑1 (MQ=255)
                                ggCCAGTTCGTCGTTGCTGTTGATGCGTACGTCACGTTCCCGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGGAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAACCCGGGGCGGCAGCCCATTGGCGAGATATc                                                                                                                                                                                           <  2:191637/250‑1 (MQ=255)
                                                  gTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCAt                                                                                                                                                                          <  1:14910/250‑1 (MQ=255)
                                                           tACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTa                                                                                                                                                                  >  1:269040/1‑249 (MQ=255)
                                                                   gTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTAc                                                                                                                                                        <  1:251426/251‑1 (MQ=255)
                                                                     tCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGc                                                                                                                                                      >  1:151277/1‑251 (MQ=255)
                                                                          aaTGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTTCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATaa                                                                                                                                                  >  1:77308/1‑250 (MQ=255)
                                                                           aTGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGAATaaa                                                                                                                                                 >  2:348370/1‑250 (MQ=255)
                                                                            tGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAAc                                                                                                                                               >  1:170450/1‑251 (MQ=255)
                                                                                aCGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGTCTCACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGAGCCTGCACTTTCAGCACGTCTTCCAGTGCCGCGTTCTAGGCATCAGCAACATGCTGCTCAGCTGGCGTGCCA‑TGATACTCAGATAAAAACCTGGGCGTCAGCCCATTGTCGAGATATCGATAATCTCTACACCGTTGCCGGTAAGATGGTTACGCATTAATCCAg                                                                                                                                            <  2:413133/250‑1 (MQ=255)
                                                                                    gcAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAAc                                                                                                                                       >  2:312073/1‑251 (MQ=255)
                                                                                     cAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGCTTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAAc                                                                                                                                       >  2:416123/1‑250 (MQ=255)
                                                                                           ccTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGc                                                                                                                                 <  1:281740/250‑1 (MQ=255)
                                                                                                                                  actgtttgtcGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCCCGTCTTCCATTGTCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAGT‑AGACTCATATAAAATCCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTTCCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGAGCCCTCTCTCTGGCATCACTTCTTTGt                                                                                         <  2:368687/242‑1 (MQ=255)
                                                                                                                                                                                 aGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGc                                          <  1:67657/251‑1 (MQ=255)
                                                                                                                                                                                         ttCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGtt                                  >  2:393891/1‑251 (MQ=255)
                                                                                                                                                                                                   gCTTTCCAGGCATCAGCAACACGCGGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTc                        <  2:151277/251‑1 (MQ=255)
                                                                                                                                                                                                                cAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACt            <  2:120314/250‑1 (MQ=255)
                                                                                                                                                                                                                   cAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAGACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCGCCAGGGTATGGATCCCTCTGTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTGCATGGTTTTCGCGACCCGAACTGCAg        <  2:83479/251‑1 (MQ=255)
                                                                                                                                                                                                                         gctgctCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCACTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCGTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCACTGTTGTCGCCACTCGAACTCCAGGCGCtt  <  2:217913/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     
ATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTT  >  NC_000913/2814222‑2814690

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