Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 431,109 A→T 11.5% intergenic (+157/+20) yajD → / ← tsx HNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913431,1090AT11.5% 161.1 / 11.4 52intergenic (+157/+20)yajD/tsxHNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K
Reads supporting (aligned to +/- strand):  ref base A (29/17);  new base T (3/3);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGAC  >  NC_000913/430967‑431257
                                                                                                                                              |                                                                                                                                                    
tGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATccc                                                                                                                                                <  3:642735/149‑1 (MQ=255)
 gCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCg                                                                                                                                               <  4:85160/149‑1 (MQ=255)
     gATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAt                                                                                                                                           <  4:187563/149‑1 (MQ=255)
       tGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAttt                                                                                                                                         <  1:322096/149‑1 (MQ=255)
         cgcTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCgg                                                                                                                                              <  2:487611/142‑1 (MQ=255)
         cgcTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCgg                                                                                                                                              >  1:487611/1‑142 (MQ=255)
                 ttATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAg                                                                                                                               >  1:391264/1‑149 (MQ=255)
                 ttATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAg                                                                                                                               >  1:411137/1‑149 (MQ=255)
                         ccTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTa                                                                                                                       <  2:391264/149‑1 (MQ=255)
                          cTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTaa                                                                                                                      >  4:131738/1‑149 (MQ=255)
                             cGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACct                                                                                                                   >  4:611949/1‑149 (MQ=255)
                                          atatatTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAac                                                                                                      >  1:400824/1‑149 (MQ=255)
                                             tatTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccac                                                                                                   >  3:267929/1‑149 (MQ=255)
                                             tatTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccac                                                                                                   >  3:424720/1‑149 (MQ=255)
                                              atTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccacc                                                                                                  >  1:661920/1‑149 (MQ=255)
                                               ttGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCg                                                                                                                                               <  1:305097/103‑1 (MQ=255)
                                               ttGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCg                                                                                                                                               >  2:305097/1‑103 (MQ=255)
                                               ttGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAccc                                                                                                 >  2:640902/1‑149 (MQ=255)
                                                  aaTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAg                                                                                              >  3:180171/1‑149 (MQ=255)
                                                  aaTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAg                                                                                              >  1:346649/1‑149 (MQ=255)
                                                         gTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTag                                                                                       >  1:611825/1‑149 (MQ=255)
                                                                  ttgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGt                                                                              >  4:387418/1‑149 (MQ=255)
                                                                   tgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGtt                                                                             <  4:424720/149‑1 (MQ=255)
                                                                   tgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGtt                                                                             >  1:613449/1‑149 (MQ=255)
                                                                   tgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGtt                                                                             <  4:267929/149‑1 (MQ=255)
                                                                    ggCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTg                                                                            <  4:649931/149‑1 (MQ=255)
                                                                    ggCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTg                                                                            <  3:131738/149‑1 (MQ=255)
                                                                         gACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgtt                                                                       <  2:39424/149‑1 (MQ=255)
                                                                            aaaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTGACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgcc                                                                    >  3:606944/1‑149 (MQ=255)
                                                                            aaaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgcc                                                                    >  4:407732/1‑149 (MQ=255)
                                                                              aGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgt                                                                  >  1:616051/1‑149 (MQ=255)
                                                                               gCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgtt                                                                 >  2:8147/1‑149 (MQ=255)
                                                                                 gTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgc                                                               >  2:622189/1‑149 (MQ=255)
                                                                                  tttACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCg                                                                                                                                               <  3:365319/68‑1 (MQ=25)
                                                                                    tACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGa                                                            >  4:12899/1‑149 (MQ=255)
                                                                                     aCGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGaa                                                           <  4:489258/149‑1 (MQ=255)
                                                                                            aTCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagt                                                    <  2:661920/149‑1 (MQ=255)
                                                                                             tCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagtt                                                   >  1:503736/1‑149 (MQ=255)
                                                                                                ggCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTg                                                >  4:265211/1‑149 (MQ=255)
                                                                                                   aTGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAt                                             >  1:469036/1‑149 (MQ=255)
                                                                                                          aaGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGtt                                      >  4:293652/1‑149 (MQ=255)
                                                                                                                            aTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccac                                                                                                   <  2:215012/70‑1 (MQ=255)
                                                                                                                            aTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccac                                                                                                   >  1:215012/1‑70 (MQ=255)
                                                                                                                             tCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGc                   >  1:339713/1‑149 (MQ=255)
                                                                                                                               gaaaaTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCa                 >  3:368889/3‑149 (MQ=255)
                                                                                                                               aGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCa                 <  1:640902/149‑1 (MQ=255)
                                                                                                                                gAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAg                <  3:12899/149‑1 (MQ=255)
                                                                                                                                gAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAg                >  1:176765/1‑149 (MQ=255)
                                                                                                                                  aaTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTa              >  1:87033/1‑149 (MQ=255)
                                                                                                                                          gAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTa      >  2:325687/1‑149 (MQ=255)
                                                                                                                                             tAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGa   <  2:346649/149‑1 (MQ=255)
                                                                                                                                             tAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGa   <  4:180171/149‑1 (MQ=255)
                                                                                                                                              aaaTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGAc  <  2:469036/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGAC  >  NC_000913/430967‑431257

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: