Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,598,529 A→T 9.4% intergenic (‑542/+88) 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,5290AT 9.4% 179.0 / 8.7 53intergenic (‑542/+88)yneO/lsrKpseudogene, AidA homolog/autoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base A (30/18);  new base T (2/3);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAGTTTATTAAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGA  >  NC_000913/1598381‑1598677
                                                                                                                                                    |                                                                                                                                                    
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)
                                        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:45006/1‑149 (MQ=255)
                                                                                                                                                ttttAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAg      >  2:298815/1‑149 (MQ=255)
                                                                                                                                                 tttAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc     <  3:771582/149‑1 (MQ=255)
                                                                                                                                                    aaCAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa  >  3:426779/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGTTTATTAAACCCCTTTTGCATCCATTAGTGTCTTTAGGAATATTCGCTATAAAATAAGGGTTGTGCTTAATGCTTTAAGAAAAATAGCAATTTTCCCTTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGA  >  NC_000913/1598381‑1598677

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: ‑