Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,598,527 T→G 9.2% intergenic (‑540/+90) yneO ← / ← lsrK pseudogene, AidA homolog/autoinducer‑2 (AI‑2) kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,598,5270TG 9.2% 162.1 / 8.6 50intergenic (‑540/+90)yneO/lsrKpseudogene, AidA homolog/autoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base T (27/18);  new base G (2/3);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTTATTAAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGC  >  NC_000913/1598381‑1598674
                                                                                                                                                  |                                                                                                                                                   
aaGTTTATTAAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTa                                                                                                                                                   <  2:422343/149‑1 (MQ=255)
 aGTTTATTAAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTaa                                                                                                                                                  <  3:195402/149‑1 (MQ=255)
   tttATTAAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACa                                                                                                                                                >  3:483350/1‑149 (MQ=255)
    ttattaAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACaa                                                                                                                                               <  4:483350/149‑1 (MQ=255)
        taAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCa                                                                                                                                           <  4:14488/149‑1 (MQ=255)
          aaCCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTGTTAAAAGCCa                                                                                                                                           >  4:330828/1‑147 (MQ=255)
          aaCCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTCCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTGTTAAAAGCCa                                                                                                                                           <  3:330828/147‑1 (MQ=255)
           aCCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGct                                                                                                                                        >  3:47590/1‑149 (MQ=255)
           aCCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTGTTAAAAGCCAGct                                                                                                                                        <  3:104223/149‑1 (MQ=255)
            ccccTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGctc                                                                                                                                       <  1:222566/149‑1 (MQ=255)
               cTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTaa                                                                                                                                    >  4:212331/1‑149 (MQ=255)
                     cATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAaagaag                                                                                                                              <  2:334169/149‑1 (MQ=255)
                        ccATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGa                                                                                                                           <  3:12541/149‑1 (MQ=255)
                          aTTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAAc                                                                                                                                                 >  3:60684/1‑125 (MQ=255)
                          aTTAGTGTCTTTAGCAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAAc                                                                                                                                                 <  4:60684/125‑1 (MQ=255)
                              gtgtCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGGTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGATGGCTTTTAACAAGCCAGCTCTAAaagaa                                                                                                                               <  2:52090/139‑1 (MQ=255)
                              gtgtCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAaagaa                                                                                                                               <  4:182732/139‑1 (MQ=255)
                              gtgtCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAaagaa                                                                                                                               >  3:182732/1‑139 (MQ=255)
                              gtgtCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAaagaa                                                                                                                               >  1:52090/1‑139 (MQ=255)
                               tgtCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGa                                                                                                                    >  3:736697/1‑149 (MQ=255)
                                gtCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGaa                                                                                                                   <  4:430353/149‑1 (MQ=255)
                                   tttAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCttt                                                                                                                                                     >  3:677387/1‑112 (MQ=255)
                                   tttAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCttt                                                                                                                                                     <  4:677387/112‑1 (MQ=255)
                                        gAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACtata                                                                                                           <  4:47590/149‑1 (MQ=255)
                                                  tataAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACt                                                                                                 >  1:314878/1‑149 (MQ=255)
                                                         tAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGcc                                                                                          >  2:798054/1‑149 (MQ=255)
                                                                       aaTGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCatgat                                                                            >  4:574241/1‑149 (MQ=255)
                                                                                aGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTc                                                                   >  2:17521/1‑149 (MQ=255)
                                                                                aGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTc                                                                   >  4:587977/1‑149 (MQ=255)
                                                                                 gAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACaa                                                                                                                                               >  4:97912/1‑72 (MQ=255)
                                                                                 gAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACaa                                                                                                                                               <  3:97912/72‑1 (MQ=255)
                                                                                    aaaTAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATg                                                               >  2:170332/1‑149 (MQ=255)
                                                                                            aTTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACtata                                                       >  3:326357/1‑149 (MQ=255)
                                                                                                  ccTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAg                                                 >  1:1890/1‑149 (MQ=255)
                                                                                                    ttGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGgc                                               >  3:200787/1‑149 (MQ=255)
                                                                                                    ttGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGgc                                               >  3:158815/1‑149 (MQ=255)
                                                                                                          atCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTc                                         >  4:771582/1‑149 (MQ=255)
                                                                                                          atCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTc                                         <  4:158815/149‑1 (MQ=255)
                                                                                                             gTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCAt                                      <  4:736697/149‑1 (MQ=255)
                                                                                                              tACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATa                                     >  1:456674/1‑149 (MQ=255)
                                                                                                                   gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAacgac                                >  3:53569/1‑149 (MQ=255)
                                                                                                                   gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAacgac                                >  3:656859/1‑149 (MQ=255)
                                                                                                                       tGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATaa                                                      >  2:808044/1‑123 (MQ=255)
                                                                                                                       tGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATaa                                                      <  1:808044/123‑1 (MQ=255)
                                                                                                                           aCGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAAc                                                     >  3:394605/1‑120 (MQ=255)
                                                                                                                           aCGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACACAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAAc                                                     <  4:394605/120‑1 (MQ=255)
                                                                                                                                gAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCAt                   <  1:170332/149‑1 (MQ=255)
                                                                                                                                         aGCTGGCTTGTTAAAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCaa                                                                 <  2:724221/94‑1 (MQ=255)
                                                                                                                                         aGCTGGCTTGTTAAAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCaa                                                                 >  1:724221/1‑94 (MQ=255)
                                                                                                                                            tGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc       >  2:427066/1‑149 (MQ=255)
                                                                                                                                            tGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc       >  3:170515/1‑149 (MQ=255)
                                                                                                                                               cTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCa    >  3:407927/1‑149 (MQ=255)
                                                                                                                                                ttttAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCa                         <  4:154402/127‑1 (MQ=255)
                                                                                                                                                ttttAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCa                         >  3:154402/1‑127 (MQ=255)
                                                                                                                                                ttttAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAg   >  2:298815/1‑149 (MQ=255)
                                                                                                                                                ttttAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAg   >  2:45006/1‑149 (MQ=255)
                                                                                                                                                 tttAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc  <  3:771582/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGTTTATTAAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGC  >  NC_000913/1598381‑1598674

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