Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 431,317 T→A 7.3% N233Y (AAT→TAT)  tsx ← 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,3170TA 7.3% 173.3 / 4.2 55N233Y (AAT→TAT) tsxnucleoside channel, receptor of phage T6 and colicin K
Reads supporting (aligned to +/- strand):  ref base T (27/24);  new base A (2/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAAC  >  NC_000913/431174‑431463
                                                                                                                                               |                                                                                                                                                  
gCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTa                                                                                                                                                >  7:214576/1‑148 (MQ=255)
gCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATAAGTa                                                                                                                                                >  6:106117/1‑148 (MQ=255)
gCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATAAGTa                                                                                                                                                >  6:110937/1‑148 (MQ=255)
  gAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACg                                                                                                                                              <  3:209194/148‑1 (MQ=255)
    aCGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCGAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGCGCCAGAATATGGCTGGAAGCGCTAGAGTTAGTAGGACggg                                                                                                                                            >  6:210494/1‑148 (MQ=255)
          aaGTTGCCGTTGCCGAGGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTa                                                                                                                                       <  8:37464/147‑1 (MQ=255)
          aaGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTaa                                                                                                                                      <  8:214576/148‑1 (MQ=255)
                     gccgAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATAAGTACGGGTCTTAATACCGTTGAt                                                                                                                            <  5:110937/147‑1 (MQ=255)
                     gccgAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATAAGTACGGGTCTTAATACCGTTGAt                                                                                                                            <  5:106117/147‑1 (MQ=255)
                      ccgAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTg                                                                                                                          <  5:27865/148‑1 (MQ=255)
                        gAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCg                                                                                                                        >  1:270992/1‑148 (MQ=255)
                              cagttcTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTAcc                                                                                                                   <  8:149490/147‑1 (MQ=255)
                                         cgtcgtTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATcccc                                                                                                       >  8:116269/1‑148 (MQ=255)
                                          gtcgtTCCAGTGAACACCGTCGTGCCAGTAACGAGCTACGACCGAGTAGTGCCAGTGATCGTAGTTCGGAGCCAGAATATGGCTGGAAGCGAGTGAGTTATTAGTACGGGTCTAAATACCGTTGATTGCGTTACCGCTGTCATcccc                                                                                                       <  1:112905/147‑1 (MQ=255)
                                               tCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGAGAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTCCCGCTGTCATCCCCTAAATc                                                                                                 >  5:31809/1‑148 (MQ=255)
                                                  aCTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCgg                                                                                               >  5:167266/1‑147 (MQ=255)
                                                     gACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAAcc                                                                                           >  1:136872/1‑148 (MQ=255)
                                                      accaccGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAAcc                                                                                           <  6:52817/147‑1 (MQ=255)
                                                          ccGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGt                                                                                      <  6:167266/148‑1 (MQ=255)
                                                           cgtcgtGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAg                                                                                       >  2:253847/1‑146 (MQ=255)
                                                                      tAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTg                                                                           >  4:157347/1‑147 (MQ=255)
                                                                        aCGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAg                                                                        <  1:253847/148‑1 (MQ=255)
                                                                         cGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTaa                                                                                                    <  2:109119/119‑1 (MQ=255)
                                                                         cGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTaa                                                                                                    >  1:109119/1‑119 (MQ=255)
                                                                           aGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGcc                                                                      >  2:241202/1‑147 (MQ=255)
                                                                           aGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGcc                                                                      >  5:166737/1‑147 (MQ=255)
                                                                                cGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATg                                                                  <  7:116269/146‑1 (MQ=255)
                                                                                cGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGt                                                                 <  2:136872/147‑1 (MQ=255)
                                                                                  aCAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGCAGCGGATGTAGc                                                              >  8:311081/1‑148 (MQ=255)
                                                                                       gtagtGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGc                                                         <  1:241202/148‑1 (MQ=255)
                                                                                           tGCCAGTGATAGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGTGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGAc                                                     >  4:60955/1‑148 (MQ=255)
                                                                                            gCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGAcc                                                    >  8:227213/1‑148 (MQ=255)
                                                                                                gTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGATCCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGccc                                                >  5:114431/1‑148 (MQ=255)
                                                                                                gTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGccc                                                >  5:39984/1‑148 (MQ=255)
                                                                                                 tGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGcccc                                               <  6:114431/148‑1 (MQ=255)
                                                                                                   aTCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTg                                                       <  1:215689/138‑1 (MQ=255)
                                                                                                   aTCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTg                                                       >  2:215689/1‑138 (MQ=255)
                                                                                                     cGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCca                                              >  2:5401/1‑145 (MQ=255)
                                                                                                      gTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAg                                           <  2:270992/147‑1 (MQ=255)
                                                                                                        aGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCca                                              >  7:6128/1‑142 (MQ=255)
                                                                                                         gTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATc                                        <  3:60955/147‑1 (MQ=255)
                                                                                                           tCAGAGCCAGAATATGGCTGGAAGCTATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTc                                                            <  6:48707/125‑1 (MQ=255)
                                                                                                           tCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTc                                                            >  5:48707/1‑125 (MQ=255)
                                                                                                           tCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCgg                                      >  2:80636/1‑147 (MQ=255)
                                                                                                            cAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGt                                     <  6:166737/147‑1 (MQ=255)
                                                                                                             agagCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGAt                                                                   >  1:163356/1‑116 (MQ=255)
                                                                                                             agagCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGAt                                                                   <  2:163356/116‑1 (MQ=255)
                                                                                                             agagCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTaa                                   >  2:296711/1‑148 (MQ=255)
                                                                                                                      atatGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACaa                          <  7:227213/148‑1 (MQ=255)
                                                                                                                       tatGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACaaa                         <  5:210494/148‑1 (MQ=255)
                                                                                                                                aaGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAAt                >  3:234940/1‑148 (MQ=255)
                                                                                                                                   cGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTg                                                                           <  3:260894/86‑1 (MQ=255)
                                                                                                                                   cGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTg                                                                           >  4:260894/1‑86 (MQ=255)
                                                                                                                                     aTAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGa           <  8:6128/148‑1 (MQ=255)
                                                                                                                                          gTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACgg      >  5:16112/1‑148 (MQ=255)
                                                                                                                                            tattaGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGt     >  1:211627/1‑147 (MQ=255)
                                                                                                                                              ttaGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAAc  <  6:39984/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAGAGTAGTGCCAGTGATCGTAGTTCAGAGCCAGAATATGGCTGGAAGCGATAGAGTTATTAGTACGGGTCTTAATACCGTTGATTGCGTTACCGCTGTCATCCCCTAAATCGGAACCCCAGTCGAAGTTGGTGAAGCCGATGTAGCTCAGCTGACCGCCCCACAGATCGGTAATCGGCACAAAGTATTTAATTTTGAAACGGTAAC  >  NC_000913/431174‑431463

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