Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,602,236 A→T 15.8% Q249H (CAA→CAT lsrA → autoinducer 2 import ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,602,2360AT15.8% 145.8 / 20.6 57Q249H (CAA→CATlsrAautoinducer 2 import ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base T (5/4);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

ATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGT  >  NC_000913/1602141‑1602334
                                                                                               |                                                                                                  
aTCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                                                                <  1:1867031/100‑1 (MQ=255)
aTCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAAt                                                                                                <  1:646735/100‑1 (MQ=255)
       tAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATAGATCTCTGCCAGCCAAAACTTATgaat                                                                                         >  1:793681/1‑97 (MQ=255)
         aGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTgg                                                                                       >  1:2285144/1‑100 (MQ=255)
          gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATATTTTTTGGCTggc                                                                                      <  1:1578170/100‑2 (MQ=255)
          gCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggc                                                                                      <  1:623802/100‑2 (MQ=255)
          gCGGCAAAACCAGCGAACTGTCTACCGACGACATCATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAACTTTGGGCTggc                                                                                      <  2:1258416/100‑2 (MQ=255)
            ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCATAATTTTTGGCTggcag                                                                                    <  2:374909/100‑4 (MQ=255)
            ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAg                                                                                     >  1:1548437/1‑99 (MQ=255)
            ggCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGt                                                                                    >  1:565044/1‑100 (MQ=255)
             gCAAAACCAGCGAACTGTCTACAGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCACCCAAAAATTATGGCTGGAg                                                                                     >  2:1057408/1‑98 (MQ=255)
              cAAAACCAGCGAACTGTCTACCGACGACATTATTCAAGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGATGAAGTTa                                                                                  >  2:412341/1‑100 (MQ=255)
                  aCCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTg                                                                              <  1:134217/100‑1 (MQ=255)
                    cAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTg                                                                              <  1:1450864/98‑1 (MQ=255)
                        gAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAAcc                                                                        <  2:855525/100‑1 (MQ=255)
                            tGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGccc                                                                    <  1:2365466/100‑1 (MQ=255)
                             gTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGc                                                                      <  1:993871/97‑1 (MQ=255)
                              tCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCcac                                                                  <  2:253640/100‑1 (MQ=255)
                                 aCCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaac                                                               >  1:582099/1‑100 (MQ=255)
                                  ccGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCacaaca                                                              <  1:2300533/100‑1 (MQ=255)
                                   cgacgaCATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTTGAGTTACCTGGTAACCGCCCACAACAt                                                             <  2:2377409/100‑1 (MQ=255)
                                        acATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGATACCTGGTAACCGCCCACAACATgccg                                                         >  2:669230/1‑99 (MQ=255)
                                        acATTATTCAGGACATCACCCCAGCGGTACGGGAAATATCGCTCTCTGCCAGCCAAAAATTATCGCTAGAGTTACCTGGTAACCGCCCACAACATgccgc                                                        >  1:1576693/1‑100 (MQ=255)
                                         cATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATgccgcc                                                       >  1:2292311/1‑100 (MQ=255)
                                             attCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGCTACCTGGTAACCGCCCACAACATGCCGCCgg                                                     >  2:729821/1‑98 (MQ=255)
                                               tCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACg                                                 <  1:1494190/100‑1 (MQ=255)
                                                   gCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCgg                                             >  2:27893/1‑100 (MQ=255)
                                                    ccATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCgg                                             <  1:2034359/99‑1 (MQ=255)
                                                      aTCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGc                                          <  1:458541/100‑1 (MQ=255)
                                                      aTCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGc                                          <  1:1693371/100‑1 (MQ=255)
                                                       tCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCt                                         <  1:2230872/100‑1 (MQ=255)
                                                        cACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTg                                        <  2:200100/100‑1 (MQ=255)
                                                         aCCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGa                                       >  1:2189372/1‑100 (MQ=255)
                                                           cccAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGaca                                     >  2:917198/1‑100 (MQ=255)
                                                             cAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACt                                   >  2:1784117/1‑100 (MQ=255)
                                                               gCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTCGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTgg                                 <  1:1379772/100‑1 (MQ=255)
                                                                 ggTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGa                                <  2:466408/99‑1 (MQ=255)
                                                                     cGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATc                           >  1:166255/1‑100 (MQ=255)
                                                                      gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCt                          >  1:586202/1‑100 (MQ=255)
                                                                      gggAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCAGGAGTTACCTGGTAACCGCCCACAAAATGCAGCCGGAACGCCGGTGCTGACACTGGAAAATCt                          >  2:1102306/1‑100 (MQ=255)
                                                                       ggAAAAATCGCTCTCTGACAGCCCAACATTCCGCCTGCAGTTACACTGTAAAAGCCAAACACCTGCCGCTGGAATGCCAGTGTTTACACTGGAAAATCTg                         <  2:2004485/100‑1 (MQ=255)
                                                                        gAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTg                         <  2:586199/99‑1 (MQ=255)
                                                                         aaaaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGAc                       >  2:1012481/1‑100 (MQ=255)
                                                                          aaaaTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGAcc                      >  1:215215/1‑100 (MQ=255)
                                                                             aTCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGc                   >  2:768150/1‑100 (MQ=255)
                                                                                gCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGaa                <  2:464823/100‑1 (MQ=255)
                                                                                gCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGaa                >  2:244355/1‑100 (MQ=255)
                                                                                    aactccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGtt            >  1:1014967/5‑100 (MQ=255)
                                                                                     actccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACAGGAAAATCTGACCGGCGAAg               >  1:615990/4‑96 (MQ=255)
                                                                                     actccagccaTAATTTTTGGATGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGttt           >  2:1513713/4‑100 (MQ=255)
                                                                                        ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACCCTGGAAAATCTGACAGGCGCAGGTTTCAg        >  1:793118/1‑100 (MQ=255)
                                                                                        ccagccaTAATTTTTGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        >  1:1842205/1‑100 (MQ=255)
                                                                                        ccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        <  2:1956868/100‑1 (MQ=255)
                                                                                        ccagccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAg        >  1:832254/1‑100 (MQ=255)
                                                                                           gccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa     >  1:1328216/1‑100 (MQ=255)
                                                                                           gccaAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGaa     <  1:1054238/100‑1 (MQ=255)
                                                                                             caAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATg   >  1:1169684/1‑100 (MQ=255)
                                                                                              aaaaaTTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGt  >  2:1215941/1‑100 (MQ=255)
                                                                                               |                                                                                                  
ATCGCCTTAAGCGGCAAAACCAGCGAACTGTCTACCGACGACATTATTCAGGCCATCACCCCAGCGGTACGGGAAAAATCGCTCTCTGCCAGCCAAAAATTATGGCTGGAGTTACCTGGTAACCGCCCACAACATGCCGCCGGAACGCCGGTGCTGACACTGGAAAATCTGACCGGCGAAGGTTTCAGGAATGT  >  NC_000913/1602141‑1602334

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: