Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,605,953 G→A Q301Q (CAG→CAA lsrB → Autoinducer‑2 ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,605,9530GA95.5% 70.0 / NA 22Q301Q (CAG→CAAlsrBAutoinducer‑2 ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/9);  minor base C (1/0);  total (13/9)
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

GATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTATGGCAGATTTAGACGATATTAAAGATGGTAAAGATTTTCGTACCGATCAACCGCAAAAAAATATCCCTTTTACCCTGAAAGGTTGCGGTGCGCT  >  NC_000913/1605711‑1606194
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
gATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCCGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGt                                                                                                                                                                                                                                           >  1:214969/1‑251 (MQ=255)
                           caAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATg                                                                                                                                                                                                                 >  1:133339/1‑250 (MQ=255)
                                  cAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGaa                                                                                                                                                                                                         <  1:105803/251‑1 (MQ=255)
                                                         gTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACt                                                                                                                                                                                  >  1:319240/1‑251 (MQ=255)
                                                                    cGGATTCAGTACGACAAATGTGATGCGCCCGTATGTAGAGCGCGGCCCGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAgc                                                                                                                                                                       >  2:48946/1‑251 (MQ=255)
                                                                           aGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTg                                                                                                                                                                                 >  2:339420/1‑234 (MQ=255)
                                                                           aGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTg                                                                                                                                                                                 <  1:339420/234‑1 (MQ=255)
                                                                                                        agagCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCaaa                                                                                                                                       <  1:21237/247‑1 (MQ=255)
                                                                                                        agagCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCaaa                                                                                                                                       >  2:21237/1‑247 (MQ=255)
                                                                                                          agCGCGGCAAGGGGAAAGAATTTGGACTGAGGGATGAGGTTCAGCACGGCAACATTTAAGAGAATGTAGCGGCCGCCTTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCACGTTGAAGTCTCGCCAAACACCGTTCCGGGCTATGACTACGAAGCGGATGGTACTGGCATCGTACTGTTACCGGAGCGCGTGATCTTCAACAAAGAGAATCACGGCCAATACgct                                                                                                                                 >  2:162277/1‑249 (MQ=255)
                                                                                                                  aCGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTAACGGAGCGCGTGATATTCAACAAAGCGAATATCGGCAAATACGATTTCTGAtg                                                                                                                         >  2:82819/1‑251 (MQ=255)
                                                                                                                     gTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGc                                                                                                                      <  1:354267/251‑1 (MQ=255)
                                                                                                                                                                      tgtATGTCGCGGATGCATTATTGACAAAAGGATCAATGCAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTAATTACCCGGAGTAAGTTATGGCAGATTTAGACGATATTAAAGATg                                                                     >  2:242720/1‑251 (MQ=255)
                                                                                                                                                                           gTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTATGGCAGATTTAGACGATATTAAAGATGGTaaa                                                                <  2:242187/251‑1 (MQ=255)
                                                                                                                                                                                            gAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTATGGCAGATTTAGACGATATTAACGATGGTAAAGAATTTCGTACCGATc                                                >  1:262531/1‑250 (MQ=255)
                                                                                                                                                                                             aaaaaaGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGTAGCGGATGGTAATGGCATCGTACTTTTACCGGAGCGTGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTATGGCAGATTTAGACGAGATTAAAGGTGGTAAAGTTTGTCGTACCGATCaa                                              <  2:203119/251‑1 (MQ=255)
                                                                                                                                                                                                  aGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGCTTTCTGATGTGCCTTACTTAACCGGAGTAAGTTATGGAGGTTTTAGGCGATATTAAAGATGGTAAAGATTTTCGTACCGATCAACCGc                                          >  2:343962/1‑250 (MQ=255)
                                                                                                                                                                                                         aTGAAAACGTGCGTCAAGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGTGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTTCTTAACCGGAGTAAGTTATGGCATATTTAGACGATATTAAAGATGGTAAAGATTTTCGTACCGATCAACCGCAAAAAAat                                  <  1:183469/251‑1 (MQ=255)
                                                                                                                                                                                                            aaaaCGGGCGACAAGCTGGATATCAAGGGCGGAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTATGGCAGATTTAGACGATATTAAAGATGGTAAAGATTTTCGTACCGATCAACCGCAAAAAAATATc                               <  2:214969/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaGCTGGATATCAAGGGCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGCAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGCGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGAAAGTTATGGCAGATTTAGACGATATTAAAGATGGTAAAGATTTTCGTACCGATCAACCGCAAAAAAATATCACTTTTACCCTg                   >  1:232993/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggCGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCCGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGCGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTATGGCAGATTTAGACGATATTAAAGATGGTAAAGATTTTCGTACCGATCAACCGCAAAAAAATATCCCTTTTACCCTGAAAGGTTGCGGTgcg    >  2:333796/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cGTAGGTCAAGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTATGGCAGATTTAGACGATATTAAAGATGGTAAAGATTTTCGTACCGATCAACCGCAAAAAAATATCCCTTTTACCCTGAAAGGTTGCGGTGCGCt  <  2:319240/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
GATGCCAACGCCCTGCCCGCTGCCGCACAAGCCGCAGAAAACTTGAAAAATGACAAAGTAGCGATTGTCGGATTCAGTACGCCAAATGTGATGCGCCCGTATGTAGAGCGCGGCACGGTGAAAGAATTTGGCCTGTGGGATGTGGTTCAGCAAGGCAAAATTTCAGTGTATGTCGCGGATGCATTATTGAAAAAAGGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTATGGCAGATTTAGACGATATTAAAGATGGTAAAGATTTTCGTACCGATCAACCGCAAAAAAATATCCCTTTTACCCTGAAAGGTTGCGGTGCGCT  >  NC_000913/1605711‑1606194

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