Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,987,655 Δ1 bp coding (120/633 nt) yqeH → putative LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,987,6540G.96.4% 239.6 / 2.8 55coding (119/633 nt)yqeHputative LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (25/28);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AAAGTCTTACATGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTGTCGATT  >  NC_000913/2987558‑2987749
                                                                                                |                                                                                               
aaaGTCTTACATGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTAtgtgccc                                                                                              >  1:322576/1‑97 (MQ=255)
aaaGTCTTACATGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGccc                                                                                              <  1:1120635/100‑1 (MQ=255)
 atgTCTTACATGCAATAGTTGCTCAATGCATTGCCAAATAATACCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTAtgtgcccg                                                                                             >  1:882499/3‑96 (MQ=255)
  aGTCTTACATGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGcccg                                                                                             <  1:433460/99‑2 (MQ=255)
     cTTACATGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCgaaga                                                                                        >  1:50587/1‑100 (MQ=255)
     cTTACATGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCgaaga                                                                                        >  2:305042/1‑100 (MQ=255)
        aCATGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAg                                                                                     <  1:25315/100‑1 (MQ=255)
          aTGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGtt                                                                                   >  2:669156/1‑100 (MQ=255)
           tGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTc                                                                                  >  1:137891/1‑100 (MQ=255)
           tGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTc                                                                                  <  1:963468/100‑1 (MQ=255)
           tGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTc                                                                                  >  1:798182/1‑100 (MQ=255)
           tGCAATAGTTGCTAAATGCATTGCCAAATAATGCCAGACGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGtt                                                                                   >  2:804171/1‑99 (MQ=255)
            gCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCACACTAGTGTTTTTGTATGTGCCCGAAGAAAGTTCa                                                                                 >  2:314416/1‑100 (MQ=255)
            gCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTc                                                                                  >  2:264823/1‑99 (MQ=255)
            gCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCa                                                                                 <  2:274377/100‑1 (MQ=255)
             cAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCat                                                                                >  1:768525/1‑100 (MQ=255)
              aaTAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCata                                                                               <  1:527610/100‑1 (MQ=255)
                  gTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATAtttt                                                                           <  1:472687/100‑1 (MQ=255)
                   ttGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTa                                                                          <  1:1207121/100‑1 (MQ=255)
                      cTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATc                                                                       >  1:147717/1‑100 (MQ=255)
                        cAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCgt                                                                     >  2:1011238/1‑100 (MQ=255)
                          aTGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCgtgg                                                                   <  2:57448/100‑1 (MQ=255)
                          aTGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCgtgg                                                                   <  2:1176100/100‑1 (MQ=255)
                            gCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCgtggtg                                                                 <  2:399173/100‑1 (MQ=255)
                             cATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCgtggtgt                                                                >  2:1118759/1‑100 (MQ=255)
                              aTTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCgtggtgtg                                                               <  1:233104/100‑1 (MQ=255)
                                tGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCgtggtgtggt                                                             <  1:381276/100‑1 (MQ=255)
                                 gCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTa                                                            >  2:424890/1‑100 (MQ=255)
                                     aataatGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAg                                                        <  1:2075/100‑1 (MQ=255)
                                       taatGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGt                                                      >  1:819571/1‑100 (MQ=255)
                                           gCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCt                                                  >  2:862958/1‑100 (MQ=255)
                                            ccAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCtt                                                 >  1:1007911/1‑100 (MQ=255)
                                               gaagaGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAAc                                              <  2:681951/100‑1 (MQ=255)
                                               gaagaGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAAc                                              <  2:440344/100‑1 (MQ=255)
                                                aagaGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACa                                             <  2:13116/100‑1 (MQ=255)
                                                aagaGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACa                                             <  1:424894/100‑1 (MQ=255)
                                                    gTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACaaaaa                                         <  1:987984/100‑1 (MQ=255)
                                                     tCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACaaaaaa                                        <  1:1011247/100‑1 (MQ=255)
                                                      cACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAAc                                       <  2:658193/100‑1 (MQ=255)
                                                        ccGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACca                                     >  1:600919/1‑100 (MQ=255)
                                                         cgcgTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACcac                                    <  1:1046517/100‑1 (MQ=255)
                                                          gcgTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACcaca                                   >  2:579292/1‑100 (MQ=255)
                                                                 cAGTATTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCgag                            <  1:745370/100‑1 (MQ=255)
                                                                     aTTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACACAACACCACAATGCGAGATTa                        >  2:1082950/1‑100 (MQ=255)
                                                                     aTTGCGCAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTa                        >  2:361133/1‑100 (MQ=255)
                                                                        gcgcAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGc                     >  1:329496/1‑100 (MQ=255)
                                                                        gcgcAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGc                     >  2:749464/1‑100 (MQ=255)
                                                                         cgcAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTAAATTAAAAAAAAACCACAATGCGAGATTAAGCg                    <  2:824366/100‑1 (MQ=255)
                                                                          gcAAACTCGTGTTTTTGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCgg                   <  2:137889/100‑1 (MQ=255)
                                                                                     tttttGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTg        <  1:514073/100‑1 (MQ=255)
                                                                                     tttttGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTg        <  2:125654/100‑1 (MQ=255)
                                                                                       tttGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTGTc      <  1:359349/100‑1 (MQ=255)
                                                                                       tttGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTGTc      <  1:932368/100‑1 (MQ=255)
                                                                                        ttGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTGTCg     <  2:742572/100‑1 (MQ=255)
                                                                                         tGTATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTGTCGa    >  2:729651/1‑100 (MQ=255)
                                                                                           tATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTGTCGAtt  >  2:200781/1‑100 (MQ=255)
                                                                                           tATGTGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTGTCGAtt  >  2:891303/1‑100 (MQ=255)
                                                                                                |                                                                                               
AAAGTCTTACATGCAATAGTTGCTCAATGCATTGCCAAATAATGCCAGAAGAGTCACCGCGTTTGCAGTATTGCGCAAACTCGTGTTTTTGTATGTGGCCCGAAGAAAGTTCATATTTTAATCGTGGTGTGGTAGAAGGTATCTTAACAAAAAACCACAATGCGAGATTAAGCGGATATATTTTTGTCGATT  >  NC_000913/2987558‑2987749

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