Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,970,436 T→C I6T (ATA→ACA)  ygdQ → UPF0053 inner membrane protein YgdQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,970,4360TC97.8% 148.1 / NA 45I6T (ATA→ACA) ygdQUPF0053 inner membrane protein YgdQ
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (20/24);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGCACAACGTGCTCATGCGCGCCGTC  >  NC_000913/2970312‑2970570
                                                                                                                            |                                                                                                                                      
cACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGc                                                                                                                         >  1:793612/1‑140 (MQ=255)
     gTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGc                                                                                                                    >  2:492562/1‑140 (MQ=255)
        cACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTg                                                                                                                 >  1:771511/1‑140 (MQ=255)
          ctactGGCCCGTCATTTTCTGACCCAGTACCCATCGCATTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTgcg                                                                                                               <  1:412386/140‑1 (MQ=255)
            actGGCCCGTCATTTTCTGATCCAGGAGCAATCGCTTGGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCt                                                                                                             <  1:413930/140‑1 (MQ=255)
                gCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCt                                                                                                             <  1:18103/136‑1 (MQ=255)
                gCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCt                                                                                                             >  2:18103/1‑136 (MQ=255)
                gCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGt                                                                                                         >  1:63483/1‑140 (MQ=255)
                gCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAAATGTTATTCGCATGGATAACCGATCCTAACACCTGGCTTCCGCTCGGt                                                                                                         >  2:413930/1‑140 (MQ=255)
                     tCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgct                                                                                                    >  2:136017/1‑140 (MQ=255)
                     tCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGAACCTGACACCTGGCGTGCGCTCGGTacgct                                                                                                    >  1:256460/1‑140 (MQ=255)
                       aTTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctga                                                                                                  <  1:183500/140‑1 (MQ=255)
                       aTTTTCTGATCCAGTAGCCATCACTATGATCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctga                                                                                                  <  1:409617/140‑1 (MQ=255)
                         tttCTGATCCAGTAGTCATCGCTTTGACCGGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctgacg                                                                                                <  1:360884/140‑1 (MQ=255)
                         tttCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctgacg                                                                                                <  1:136017/140‑1 (MQ=255)
                         tttCTGATCCAGTAGCCATCGCCTTTACCTGCCGCTTCCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTacgctgacg                                                                                                <  1:306993/140‑1 (MQ=255)
                                tCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTgg                                                                                         >  1:21574/1‑140 (MQ=255)
                                     tAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATc                                                                                    >  2:863161/1‑140 (MQ=255)
                                            cccTTTGACCTGCCGCTCTCCGCGCATATAATTACCGCTCCATTTTTTCGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTg                                                                             <  2:256460/138‑1 (MQ=255)
                                              cTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTggg                                                                           <  1:492562/140‑1 (MQ=255)
                                                      tGCCGCTTTCCGGGCATATAATTACGGCTTCATTTTTGTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACaa                                                                   <  1:298843/140‑1 (MQ=255)
                                                       gCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGa                                                                      >  1:616189/1‑136 (MQ=255)
                                                       gCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGa                                                                      <  2:616189/136‑1 (MQ=255)
                                                        ccGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAata                                                                 <  1:863161/140‑1 (MQ=255)
                                                          gCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAAtatt                                                               <  2:793612/140‑1 (MQ=255)
                                                               ccGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTAtttt                                                          <  1:729781/140‑1 (MQ=255)
                                                                  ggCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCt                                                       <  2:771511/140‑1 (MQ=255)
                                                                       ataATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCGTGGGATCGACAATATTATTTTCCTTtctc                                                  >  1:324161/1‑140 (MQ=255)
                                                                           ttACCGCTTCATTTTTTTGGCATGGTTTTGTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCtggt                                              <  1:281590/140‑1 (MQ=255)
                                                                             aCCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGTCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCtggtggt                                           >  1:11548/1‑140 (MQ=255)
                                                                               cGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCtggtggtg                                          >  1:47791/1‑140 (MQ=255)
                                                                                 cTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGTGATCGTTCATGGGATCCACAATATTATTTTGCTTTCTCTGGTGGTGGc                                        >  2:255909/1‑140 (MQ=255)
                                                                                 cTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGCTCGTTCTTGGGATCGAAAATATTGTTTTCCGTTCTCTGGTGGTGGc                                        >  1:413972/1‑140 (MQ=255)
                                                                                 cTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGc                                        >  1:198177/1‑140 (MQ=255)
                                                                                  ttCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACATCCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCa                                       <  2:640931/140‑1 (MQ=255)
                                                                                   tCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCaa                                      >  2:829182/1‑140 (MQ=255)
                                                                                          tttGGCAGGGCTTTTCAGATGTCATATGCCTGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTcc                               <  1:255909/140‑1 (MQ=255)
                                                                                          tttGGCAGCGCTTTTTGCATGTTATTTGCAAAGTCTACCTATCCTAACACCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTcc                               <  2:413972/140‑1 (MQ=255)
                                                                                               cAGGGCTTTTTAGATGTTATTTCCATGGACATCCGATCCAAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAg                          <  2:324161/140‑1 (MQ=255)
                                                                                                aGGGCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGc                         >  1:633750/1‑140 (MQ=255)
                                                                                                 gggCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGca                        <  1:624334/140‑1 (MQ=255)
                                                                                                   gCTTTTTAGATGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGcaca                      <  1:819026/140‑1 (MQ=255)
                                                                                                             tGTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGCACAACGTGCTCAt            >  2:481008/1‑140 (MQ=255)
                                                                                                              gTTATTTGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGCACAACGTGCTCATg           >  1:537827/1‑140 (MQ=255)
                                                                                                                  tttGCATGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGCACAACGTGCTCATgcgcg       <  1:481008/140‑1 (MQ=255)
                                                                                                                       aTGGACAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGCACAACGTGCTCATGCGCGCCGTc  <  2:198177/140‑1 (MQ=255)
                                                                                                                            |                                                                                                                                      
CACCAGTGCACTACTGGCCCGTCATTTTCTGATCCAGTAGCCATCGCTTTGACCTGCCGCTTTCCGGGCATATAATTACCGCTTCATTTTTTTGGCAGGGCTTTTTAGATGTTATTTGCATGGATAACCGATCCTAACGCCTGGCTTGCGCTCGGTACGCTGACGCTGCTGGAGATCGTTCTTGGGATCGACAATATTATTTTCCTTTCTCTGGTGGTGGCAAAGCTTCCCACAGCACAACGTGCTCATGCGCGCCGTC  >  NC_000913/2970312‑2970570

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