Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,814,459 Δ1 bp coding (275/516 nt) luxS ← S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,4590A.96.3% 118.6 / ‑3.2 27coding (275/516 nt)luxSS‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (12/14);  total (12/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 = 1.00e+00

GTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGG  >  NC_000913/2814229‑2814685
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
gTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCgg                                                                                                                                                                                                               >  1:55190/1‑251 (MQ=255)
 ttCCTGCAACTTCTCTTTCGGCAGTGCAAGTTCTTCGTTGCTGTTGATGAGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGGGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCGTTCCAGGCATCATCAACACGCTGCTCATCCGGCGTACCATTAGAATCATATAAAAAACggg                                                                                                                                                                                                              >  2:283607/1‑250 (MQ=255)
            gctctTCCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGccc                                                                                                                                                                                                    <  1:165029/249‑1 (MQ=255)
                                     gTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTac                                                                                                                                                                          >  1:90443/1‑251 (MQ=255)
                                                 gCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCATCCTGAAAAGAGTGCATCTGGTAAGTGCCACACTGGAAGACGTTCAGTACCGGGCTCTGATTCTGATCCTGCAATTTCAGCACGTCTTCCATTGCCGCTTTCCAGTCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTACGGCAGCCCATTGGCGAGATATCGCTAATCTCTACACCATTACAGtt                                                                                                                                                              >  1:264079/1‑251 (MQ=255)
                                                       gTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATg                                                                                                                                                        >  1:177958/1‑251 (MQ=255)
                                                            gTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCCTCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTAc                                                                                                                                                   >  1:194300/1‑251 (MQ=255)
                                                            gTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTAc                                                                                                                                                   <  2:888/251‑1 (MQ=255)
                                                                                atatCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAAGCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAgg                                                                                                                               <  2:155015/251‑1 (MQ=255)
                                                                                             ccTGCAACGAGTGCATCTGGTAAGTGCCACACTGGGAGACGTTCAGTGCCGGGATCTGAGTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTTCTCATCTGGCGTACCATTAGACGCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTa                                                                                                                   <  2:277710/250‑1 (MQ=255)
                                                                                                      aGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGTTCCctc                                                                                                         >  2:248846/1‑251 (MQ=255)
                                                                                                                  aaGTGCCACTCTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCa                                                                                             <  1:177094/251‑1 (MQ=255)
                                                                                                                  aaGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATc                                                                                              >  1:3952/1‑250 (MQ=255)
                                                                                                                                         aGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGa                                                                        >  1:153642/1‑249 (MQ=255)
                                                                                                                                            tCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAgcgc                                                                   >  2:66038/1‑251 (MQ=255)
                                                                                                                                                      gattctgctCCTACACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACTCCATTACCGTTAAGATGGTTACGCATAAATCCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTGCTTTGTTCGGCACGCAGAAGCGCAGATCg                                                             <  1:262931/247‑1 (MQ=255)
                                                                                                                                                      gattctgatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCTCTTCTTTGTTCGGCACGCAGAAGCGCAGATCg                                                             <  2:55190/247‑1 (MQ=255)
                                                                                                                                                      gattctgatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATc                                                              <  2:205837/246‑1 (MQ=255)
                                                                                                                                                       attctgatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATAGAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACTCATAAGCGCAGATCGAacac                                                        <  2:303410/251‑1 (MQ=255)
                                                                                                                                                             atCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGAt                                                  >  1:103833/1‑251 (MQ=255)
                                                                                                                                                               ccTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGAtt                                                 <  1:279673/250‑1 (MQ=255)
                                                                                                                                                                           gCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCg                                    <  1:23705/251‑1 (MQ=255)
                                                                                                                                                                             ccggtttCCATTGCCGCTTTCCAGGCAACAGCAACACGCTGCTCATCTGGCGTACCAATAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCggg                                  <  1:107076/246‑1 (MQ=255)
                                                                                                                                                                                         gccgccctCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGGGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCGCCAGGGTATGGATCCCTCTTTCTGGCATCACTGCGGTGTTCGGCACGCAGAAGCGTAGATCGAACACGGTAATTGCGCTGCCATGCGGTGAGTTCAttat                       <  2:178515/243‑3 (MQ=255)
                                                                                                                                                                                                         cAGCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTg      <  1:155336/251‑1 (MQ=255)
                                                                                                                                                                                                           gCAACACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTTTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCa    <  1:135557/251‑1 (MQ=255)
                                                                                                                                                                                                             aaCACGCTGCTCATCTGGCGTACCATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAgg  >  1:158113/1‑251 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
GTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGG  >  NC_000913/2814229‑2814685

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