Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,604,698 A→G T218A (ACG→GCG)  lsrD → autoinducer 2 import system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,6980AG98.2% 181.5 / ‑4.2 55T218A (ACG→GCG) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/28);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

ATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCA  >  NC_000913/1604551‑1604842
                                                                                                                                                   |                                                                                                                                                
aTCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCg                                                                                                                                                <  1:147946/150‑1 (MQ=255)
 tCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCg                                                                                                                                                <  2:828563/149‑1 (MQ=255)
    tCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGc                                                                                                                                            <  2:642002/150‑1 (MQ=255)
     ccTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTg                                                                                                                                          <  1:101768/151‑1 (MQ=255)
      cTGATATGTCTCCTCGGTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTgg                                                                                                                                         >  2:631354/1‑151 (MQ=255)
        gATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCg                                                                                                                                       <  2:750918/151‑1 (MQ=255)
             ttctcCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTg                                                                                                                                   >  1:721219/2‑150 (MQ=255)
              tctcCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCg                                                                                                                                 >  1:46062/1‑151 (MQ=255)
                 cctcGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCgg                                                                                                                                <  2:505954/149‑1 (MQ=255)
                   tcGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCgg                                                                                                                                >  1:368972/1‑147 (MQ=255)
                    cGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTcg                                                                                                                             <  2:418590/149‑1 (MQ=255)
                     gTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGGcgcc                                                                                                                           >  2:518929/1‑150 (MQ=255)
                              ctctggctGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCtgttgt                                                                                                                   >  2:528462/1‑146 (MQ=255)
                               tctggctGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGtgctgc                                                                                                                 >  2:508814/1‑150 (MQ=255)
                                      gCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTAt                                                                                                           <  2:791966/149‑1 (MQ=255)
                                        aTAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTa                                                                                                       <  2:286375/151‑1 (MQ=255)
                                            aaCCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGGATCGTAtttt                                                                                                   >  2:96715/1‑151 (MQ=255)
                                               ccATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGt                                                                                                <  1:839371/151‑1 (MQ=255)
                                                cATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTg                                                                                                                                          >  1:533387/1‑108 (MQ=255)
                                                cATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGc                                                                                                                                        <  2:533424/110‑1 (MQ=255)
                                                cATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGtt                                                                                               <  1:594162/151‑1 (MQ=255)
                                                      ggACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCAc                                                                                         >  1:751701/1‑151 (MQ=255)
                                                      ggACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGAGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTTTGCGGTCGCCGATGGGCTGCTGGTTTCGCATTTTTGGTCaac                                                                                           >  2:238029/1‑147 (MQ=255)
                                                          gTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACggt                                                                                      >  1:513516/1‑148 (MQ=255)
                                                            aaTGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCg                                                                                   <  1:526658/151‑1 (MQ=255)
                                                                tgtTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGAtct                                                                               >  2:252561/1‑151 (MQ=255)
                                                                  tttttGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGCATCGTATTTTGGTTCAGCACGTTCCGATCTCg                                                                             >  1:794241/1‑151 (MQ=255)
                                                                    tttGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCgg                                                                            <  1:715812/150‑1 (MQ=255)
                                                                       gATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGcg                                                                        <  2:206457/151‑1 (MQ=255)
                                                                           gggCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGcgtcg                                                                     <  2:451005/150‑1 (MQ=255)
                                                                              cAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGGGCGTCGttt                                                                  >  1:176607/1‑150 (MQ=255)
                                                                               aaaGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGGATTTTGGGTCAGCACGTTCCGATCTCGGTGCGTCGTTTCt                                                                >  2:357039/1‑151 (MQ=255)
                                                                                 aGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGa                                                              <  1:513773/151‑1 (MQ=255)
                                                                                  gCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGAt                                                             <  2:811734/151‑1 (MQ=255)
                                                                                    ccGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGc                                                           <  1:235152/151‑1 (MQ=255)
                                                                                    ccGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGc                                                           <  1:508778/151‑1 (MQ=255)
                                                                                          gTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCa                                                     <  2:21864/151‑1 (MQ=255)
                                                                                                          cgATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGc                                      >  2:807821/1‑150 (MQ=255)
                                                                                                              tCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTgg                                  >  1:242289/1‑150 (MQ=255)
                                                                                                              tCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTgg                                  >  2:22077/1‑150 (MQ=255)
                                                                                                               ccAGTTAACCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCt                                     >  2:380271/1‑146 (MQ=255)
                                                                                                                     acccGTCCCTTATGTGCGTTCTAGGCCGGGGGGGGGCTGGCGTCTGCGGTCGCCGCTGGGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGc                            <  1:283403/147‑1 (MQ=255)
                                                                                                                     aaCCGTCCCTTATGTGCGCTCTAGGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCgtggtg                                       <  1:358362/138‑1 (MQ=255)
                                                                                                                     aaCCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCgtggtg                                       >  2:358385/1‑138 (MQ=255)
                                                                                                                     aaCCGTACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGc                            <  2:405961/149‑1 (MQ=255)
                                                                                                                         gtcccTTATGTGCTCTCTATGCCGTGGGGGGGCGGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAatat                      <  1:211887/148‑1 (MQ=255)
                                                                                                                          tACCTTATGTGCGCTCTATGCCATGGCGGGGCTGGCGTCTGCGGGCGCCGCTTTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTTTTATGCCCGCCCTCACCGCCGTGGTGCTTGGCGGGGCCAATAtt                     >  2:793027/1‑151 (MQ=255)
                                                                                                                           aCCTTATGTGCGCTCTATGCCATGGCGGGGCCGGCGTCTGCGGTCGCCGCTGGGCTGCTGGGATCGTATTTTGGGTCAGCCCGTTCCGATCTCGGGGGGTCGGTTCTGATTCACGCCATCACCGACGTGGAGCTTGGCGGGGGCCATAttt                    >  1:793592/1‑151 (MQ=255)
                                                                                                                                    gcgcTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGGATCGTATTTTGGTTCAGCACGTTCCGATCTCGGGGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTAtggtgg             >  1:175214/1‑149 (MQ=255)
                                                                                                                                      gcTCTATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTAtggtgggttc         >  2:265352/1‑147 (MQ=255)
                                                                                                                                       ctctATGCATTGGCGGGGCGGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCCCGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTggttcc         <  1:807765/150‑1 (MQ=255)
                                                                                                                                         ctATGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCCCGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGGCAATATTTATGGTggttccg        <  2:698001/149‑1 (MQ=255)
                                                                                                                                            tGCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTgggtccggttcc   >  2:458792/1‑151 (MQ=255)
                                                                                                                                             gCCATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCa  <  2:183059/151‑1 (MQ=255)
                                                                                                                                              ccATGGCGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGGTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTggttccggttcc   >  2:11242/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATCTTCCTGATATGTCTCCTCGTTTTCTGGCTCTGGCTGCATAAAACCCATGCCGGACGTAATGTGTTTTTGATTGGGCAAAGCCCGCGCGTGGCGCTTTATAGCGCGATTCCAGTTAACCGTACCTTATGTGCGCTCTATGCCATGACGGGGCTGGCGTCTGCGGTCGCCGCTGTGCTGCTGGTATCGTATTTTGGTTCAGCACGTTCCGATCTCGGTGCGTCGTTTCTGATGCCCGCCATCACCGCCGTGGTGCTTGGCGGGGCCAATATTTATGGTGGTTCCGGTTCCA  >  NC_000913/1604551‑1604842

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