Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,674,403 T→C 9.5% intergenic (+28/‑383) 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,4030TC 9.5% 151.2 / 4.5 50intergenic (+28/‑383)yhjD/yhjEinner membrane putative BrbK family alternate lipid exporter/putative MFS transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base T (18/27);  new base C (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

GTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCATT  >  NC_000913/3674256‑3674550
                                                                                                                                                   |                                                                                                                                                   
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: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                                                                                                   <  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)
                                                                                                                                         gTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAAc           >  4:252938/1‑149 (MQ=255)
                                                                                                                                          tAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACa          >  1:735917/1‑149 (MQ=255)
                                                                                                                                               gCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTc     <  1:142957/149‑1 (MQ=255)
                                                                                                                                                ccTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCa    <  2:634856/149‑1 (MQ=255)
                                                                                                                                                  ttATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCAtt  >  1:441185/1‑149 (MQ=255)
                                                                                                                                                  ttATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCAtt  <  2:407088/149‑1 (MQ=255)
                                                                                                                                                  ttATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCAtt  <  1:321882/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGCTGGGGCTGATGGCATTTTTCTACTTCTTCGCCCGTTTGACGCTGTTTTGCGCGGCGTGGATTGCCACCGCCGAATATAAAGACGACCCGCGAATGCCGGGGAAAACGCAGCCTTAAAATAGGCCGGATGCGGCGTAAACGCCTTATCCGGCATACGATTTTATTCGGCCTGGCTCCCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCTATAATGCGCACATAACCTCTTGAAACTCGTCCCCAGAGCCTCTTCAGCCATCTATTTGGGAGCAAACAATTTCATT  >  NC_000913/3674256‑3674550

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