Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,601,831 (A)7→6 coding (342/1536 nt) lsrA → Autoinducer‑2 ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,601,8250A.96.0% 93.0 / ‑1.9 25coding (336/1536 nt)lsrAAutoinducer‑2 ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (11/13);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGTATTTATCTCGTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGG  >  NC_000913/1601748‑1601901
                                                                             |                                                                            
ggTATTTATCTCGTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCaa                                                                             <  2:296729/79‑1 (MQ=255)
    tttATCTCGTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCaaaaaac                                                                        <  2:296730/80‑2 (MQ=255)
         ctcGTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGctc                                                                  >  1:296770/1‑81 (MQ=255)
            gTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGctc                                                                  <  1:296731/77‑1 (MQ=255)
            gTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTcc                                                               <  1:296732/80‑1 (MQ=255)
            gTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTcc                                                               <  2:296733/80‑1 (MQ=255)
                       aCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATg                                                   <  2:296734/81‑1 (MQ=255)
                         cGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGaaaaatga                                                 <  1:296735/81‑5 (MQ=255)
                         cGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATg                                                   >  1:296769/1‑79 (MQ=255)
                          gctgctTTTCCCAAGCCTGTCGATAAAAGAAGACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATgaag                                                <  2:296736/81‑1 (MQ=255)
                            tgctTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATgaagaa                                              >  2:296765/1‑81 (MQ=255)
                                tttCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTg                                          >  1:296752/1‑81 (MQ=255)
                                  tCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGc                                         <  1:296737/80‑1 (MQ=255)
                                     cAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTggc                                     >  2:296767/1‑81 (MQ=255)
                                      aaGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTggcg                                    <  2:296738/81‑1 (MQ=255)
                                       aGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTggcgg                                   >  1:296761/1‑81 (MQ=255)
                                             tCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTgg                             <  1:296740/81‑1 (MQ=255)
                                             tCGACAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCt                               <  2:296739/79‑1 (MQ=255)
                                                 tAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTg                         <  2:296741/81‑1 (MQ=255)
                                                   aaaGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGc                        <  1:296742/80‑1 (MQ=255)
                                                    aaGAAAACATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGc                        >  1:296776/1‑79 (MQ=255)
                                                          aCATCCTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGa                >  1:296785/1‑81 (MQ=255)
                                                               cTGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTg            <  2:296743/80‑1 (MQ=255)
                                                                tGTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCa          >  1:296779/1‑81 (MQ=255)
                                                                 gTTTGGGCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCa          >  1:296768/1‑80 (MQ=255)
                                                                     gggCTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTc     >  2:296789/1‑81 (MQ=255)
                                                                        cTGGCAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTgg  >  2:296775/1‑81 (MQ=255)
                                                                             |                                                                            
GGTATTTATCTCGTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCTGG  >  NC_000913/1601748‑1601901

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