Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,674,388 G→A 10.4% intergenic (+13/‑398) yhjD → / → yhjE inner membrane putative BrbK family alternate lipid exporter/putative MFS transporter; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,674,3880GA10.4% 139.1 / 6.5 48intergenic (+13/‑398)yhjD/yhjEinner membrane putative BrbK family alternate lipid exporter/putative MFS transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (15/28);  new base A (2/3);  total (17/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GCAGCATTTGGTTCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAG  >  NC_000913/3674241‑3674536
                                                                                                                                                   |                                                                                                                                                    
gcagcaTTTGGTTCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATgc                                                                                                                                                     >  4:647767/1‑149 (MQ=255)
    caTTTGGTTCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAACTAGGCCGGATgcggcg                                                                                                                                                 >  3:113098/1‑149 (MQ=255)
        tGGTTCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTaaa                                                                                                                                             <  1:758849/149‑1 (MQ=255)
         ggTTCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAAc                                                                                                                                            <  1:703546/149‑1 (MQ=255)
           ttCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGc                                                                                                                                          <  3:298024/149‑1 (MQ=255)
               gTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                      <  4:508356/149‑1 (MQ=255)
                tGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                     <  3:379167/149‑1 (MQ=255)
                     gggCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCCCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                <  2:167531/149‑1 (MQ=255)
                             ggCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGAtt                                                                                                                        <  1:483673/149‑1 (MQ=255)
                                aTTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTa                                                                                                                     >  3:22613/1‑149 (MQ=255)
                                   tttCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTc                                                                                                                  >  1:723344/1‑149 (MQ=255)
                                      cTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGc                                                                                                               <  3:51942/149‑1 (MQ=255)
                                       tACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATAAGGCGTTTACGCCGCATCCGGcc                                                                                                                               >  4:754890/1‑131 (MQ=38)
                                       tACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATAAGGCGTTTACGCCGCATCCGGcc                                                                                                                               <  3:754890/132‑2 (MQ=38)
                                        aCTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCt                                                                                                             <  1:452198/149‑1 (MQ=255)
                                        aCTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCt                                                                                                             <  2:527386/149‑1 (MQ=255)
                                          ttcttcGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTgg                                                                                                           <  2:723344/149‑1 (MQ=255)
                                                gCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTcccc                                                                                                     >  4:364155/1‑149 (MQ=255)
                                                gCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTcccc                                                                                                     <  3:499244/149‑1 (MQ=255)
                                                  ccGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGt                                                                                                   >  2:366164/1‑149 (MQ=255)
                                                    cgtttaCGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAg                                                                                                 <  4:213794/144‑1 (MQ=255)
                                                                  tGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATgc                                                                                   <  4:3912/149‑1 (MQ=255)
                                                                    cgcgGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATgcgc                                                                                 <  4:113098/149‑1 (MQ=255)
                                                                    cgcgGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATgcgc                                                                                 <  4:316570/149‑1 (MQ=255)
                                                                         cgTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAgca                                                                            <  2:71586/149‑1 (MQ=255)
                                                                               ttGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCt                                                                                                             <  3:569528/110‑1 (MQ=255)
                                                                               ttGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCt                                                                                                             >  4:569528/1‑110 (MQ=255)
                                                                               ttGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAgcatcgcat                                                                      >  3:742283/1‑149 (MQ=255)
                                                                                  ccACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCg                                                                   >  1:569969/1‑149 (MQ=255)
                                                                                   cACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGTCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCgg                                                                  >  2:321882/1‑149 (MQ=255)
                                                                                      cgccgAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCGTCCGGCta                                                               >  2:666887/1‑149 (MQ=255)
                                                                                      cgccgAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCta                                                               >  4:498230/1‑149 (MQ=255)
                                                                                                aaaGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGAt                                                                                     <  4:447486/117‑1 (MQ=255)
                                                                                                aaaGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGAt                                                                                     >  3:447486/1‑117 (MQ=255)
                                                                                                aaaGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGcac                                                     >  4:70239/1‑149 (MQ=255)
                                                                                                  aGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACAt                                                   >  1:201280/1‑149 (MQ=255)
                                                                                                    acgacCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATaa                                                 <  2:201280/149‑1 (MQ=255)
                                                                                                      gacCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAAcc                                               <  4:218433/149‑1 (MQ=255)
                                                                                                         ccGCGAATGCCGGGGAAAACACAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACctct                                            >  3:70400/1‑149 (MQ=255)
                                                                                                                               aGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCctc                      <  2:178951/149‑1 (MQ=255)
                                                                                                                                  cTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATaa                                                 <  3:64651/119‑1 (MQ=255)
                                                                                                                                  cTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATaa                                                 >  4:64651/1‑119 (MQ=255)
                                                                                                                                  cTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATGACCTCTTGAAACTCGTCCCCAGAGCCTCTTc                   <  3:634550/149‑1 (MQ=255)
                                                                                                                                   ttAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCa                  <  1:366164/149‑1 (MQ=255)
                                                                                                                                      aaaTAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGcc               <  3:364155/149‑1 (MQ=255)
                                                                                                                                           ggCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTa          <  4:70400/149‑1 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCACATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCgg                                                                  <  1:159791/92‑1 (MQ=25)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCACATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCgg                                                                  >  2:159791/1‑92 (MQ=25)
                                                                                                                                              cGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTAttt       <  3:70239/149‑1 (MQ=255)
                                                                                                                                                gATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTgg     <  3:498230/149‑1 (MQ=255)
                                                                                                                                                   gcggcgTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAg  >  1:115345/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGCATTTGGTTCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAG  >  NC_000913/3674241‑3674536

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