Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,674,395 A→T 10.0% intergenic (+20/‑391) 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,3950AT10.0% 164.5 / 6.8 50intergenic (+20/‑391)yhjD/yhjEinner membrane putative BrbK family alternate lipid exporter/putative MFS transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base T (2/3);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTTCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACA  >  NC_000913/3674249‑3674542
                                                                                                                                                  |                                                                                                                                                   
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                                                                                                                   <  2:527386/149‑1 (MQ=255)
                                aCTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCt                                                                                                                   <  1:452198/149‑1 (MQ=255)
                                  ttcttcGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTgg                                                                                                                 <  2:723344/149‑1 (MQ=255)
                                        gCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTcccc                                                                                                           <  3:499244/149‑1 (MQ=255)
                                        gCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTcccc                                                                                                           >  4:364155/1‑149 (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:316570/149‑1 (MQ=255)
                                                            cgcgGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATgcgc                                                                                       <  4:113098/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                                                                                           >  3:447486/1‑117 (MQ=255)
                                                                                        aaaGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGAt                                                                                           <  4:447486/117‑1 (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                                                                        >  2:159791/1‑92 (MQ=25)
                                                                                                                                    gCCGGATAAGGCGTTTACGCCACATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCgg                                                                        <  1:159791/92‑1 (MQ=25)
                                                                                                                                      cGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTAttt             <  3:70239/149‑1 (MQ=255)
                                                                                                                                        gATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTgg           <  3:498230/149‑1 (MQ=255)
                                                                                                                                           gcggcgTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAg        >  1:115345/1‑149 (MQ=255)
                                                                                                                                                gTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAAc   >  4:252938/1‑149 (MQ=255)
                                                                                                                                                 tAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACa  >  1:735917/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGTTCCGTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACA  >  NC_000913/3674249‑3674542

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