Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,140 A→G 9.3% intergenic (+149/+38) sibC → / ← serA small RNA SibC/phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1400AG 9.3% 100.7 / 3.5 43intergenic (+149/+38)sibC/serAsmall RNA SibC/phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (23/16);  new base G (2/2);  total (25/18)
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.81e-02

CCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACAT  >  NC_000913/3056999‑3057278
                                                                                                                                             |                                                                                                                                          
ccccATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGAATTACTTTATCGAGCACACTAGCAAAGTACAATGCGCAACGCGGCAAAGACGAGGACAAGCAATAAACACACGAGCaat                                                                                                                                           >  2:351590/1‑142 (MQ=255)
        tgccctgccATCGTAAATCCCCATTAAATAAACACAACGCAGTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAc                                                                                                                                   <  1:272418/143‑1 (MQ=255)
            cggccATCGTAGCTCCCCATTAACTAAGCACAACGCATCGATCTGACTTTGATTTATTTTCGGGAGCAGACTCGCAAAGGAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGc                                                                                                                               <  2:541154/141‑1 (MQ=255)
                cATCGTAAATCCCCATTAAATAAACACACCGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGt                                                                                                                           <  1:207465/143‑1 (MQ=255)
                  tCGTAAATCCGCATTAAATAAACACATCGCATTGATCTGACTTTGATGTATTTGCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttt                                                                                                                         <  2:128521/143‑2 (MQ=255)
                       aaTCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGGGACCTGCCCGTGGAtgtt                                                                                                                    >  1:153044/1‑143 (MQ=255)
                        agccccATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                                                   <  1:583813/141‑1 (MQ=255)
                              aTTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAGCTGCGCGTTGGTTTTCggaggg                                                                                                             >  1:89307/1‑141 (MQ=255)
                               ttAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                                                            >  1:361705/1‑143 (MQ=255)
                                   ataaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGa                                                                                                        >  2:846709/1‑143 (MQ=255)
                                    taaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGCTGATTTTCAGAGAAGGGGaa                                                                                                       >  1:289571/1‑143 (MQ=255)
                                     aaaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACTGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAACAGGGGAAt                                                                                                      <  2:152952/143‑1 (MQ=255)
                                      aaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                                                                     >  2:638808/1‑143 (MQ=255)
                                       tcctaaCCCATTGATCTGACTTTGATTTATTTTCTGGAGCGGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTa                                                                                                    <  2:52037/139‑1 (MQ=255)
                                        cacaACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGCATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAGACGGGCGAGTCGGTGACCTGGCAGTGGATTTTCAGAGGAGGGGAAGTat                                                                                                   >  2:882387/1‑142 (MQ=255)
                                                    aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGAAATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGTAATTAGTCGAGCAGACgg                                                                                       >  2:770912/1‑143 (MQ=255)
                                                      cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATCTTCAGAGAGGGGGAGTTAGTACAGCAGACGGgc                                                                                     >  2:806031/1‑143 (MQ=255)
                                                       ggACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAACAGTCCAGCAGACGGgcg                                                                                    >  2:871362/2‑143 (MQ=255)
                                                                   atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGATGATTTTCAGCGAAGGGGAATTAGTACAGCAGACGGGCGCGAACGGGAccc                                                                        >  2:551547/1‑143 (MQ=255)
                                                                                   tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAGGGTACCCGGAATAGCTTTCCTTg                                                        >  1:297938/1‑143 (MQ=255)
                                                                                    cGCAAAGTAGAATGCGCAACGCGGCGACGGTGGCGAGAAGGGATAAAAAAATGGGCAAGTCAGTGACCTGCCCGTTGATGTTCAGAGAAGGTGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                                                       <  2:36578/143‑1 (MQ=255)
                                                                                    cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                                                       <  1:666968/143‑1 (MQ=255)
                                                                                    cGCAAAGTACAAGGCGCAACGCGGCAACGGTGCGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                                                       <  2:545916/143‑1 (MQ=255)
                                                                                          gcaggaTGCGCAAGGCGGCAAGGGTGTGGAGAAGGGATAAAAAGACGGGCAAGTCAGTGACCTGCCCGGTGATTTTCAGAGAAGGGAAATTAGTACTGCAGACGGGCGCTAATGGTACCCGGAATAGCTTTCATTGCCTGCAg                                                 <  1:866717/138‑1 (MQ=255)
                                                                                                  cgcAGCGCGGCAACGGTGTGGAGAAGGGATCGAAGATGGGGCCAGTCAGTGACCTGCCCGTTGATTTTCAGAAAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt                                         <  2:216908/143‑1 (MQ=255)
                                                                                                   gcAACGCGGCAACGGTGGGGAGAAGGGATAAGAACCCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGACGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTc                                        <  1:85361/143‑1 (MQ=255)
                                                                                                   gcAACGCGGCAACGGTGAGGAGAAGATAAAAAAAAACGGGCAAGTCAGTGACCTGCCCGATGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCCGCGCGTGTTc                                        >  2:538420/1‑143 (MQ=255)
                                                                                                            cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGCATGGTACCCGCCATAGCCTTCATTGCCTGCCGCGCTTTTTCGGCAACGTc                               >  1:127613/1‑143 (MQ=255)
                                                                                                                  tgtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTtcgtc                         >  1:399399/1‑143 (MQ=255)
                                                                                                                     ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTCTTTCGGCAACGTCTCCGCCGGc                      >  1:26939/1‑143 (MQ=255)
                                                                                                                      gagaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCt                     <  1:732827/143‑1 (MQ=255)
                                                                                                                          aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACTGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaa                 <  2:425471/143‑1 (MQ=255)
                                                                                                                             gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCGACGTCTTCGTCGGCTtcaatat              >  1:698535/1‑143 (MQ=255)
                                                                                                                              aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAATAAAGCAGATGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatc             >  1:125528/1‑143 (MQ=255)
                                                                                                                                aaaaaaaCTGGTAAGTCAGTGACCTTCCCGTTGATTTTCATAGATGGTGAATTAGTACAGCAGACGTGTGCTAATGGTACCCTGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTTGGCTtcaatatcaa           <  2:266272/143‑1 (MQ=255)
                                                                                                                                  aaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCCACGTCTTCGTCGGCTtcaatatccatc         >  1:389563/1‑142 (MQ=255)
                                                                                                                                  aaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTtcaatatcaatt         >  1:383884/1‑142 (MQ=255)
                                                                                                                                    caaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTGTTCGGCAACGTCTTCGTCGGCTTCAGTATCGATAAc       >  1:248666/2‑143 (MQ=255)
                                                                                                                                    aagcGGGCAAGTCAGTTACCTGTCCGTTGATATGCAGAGAAGGGGAAGTAGGACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCGTTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAAc       <  1:70028/140‑1 (MQ=255)
                                                                                                                                    aaaCGGGCCAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGCCGGGCACGAATGGTACCCGGAATAGCTGTCATTGCCTGCAGCGCTCTTTCGGCCACGGCGTCGGCGGCGTCGATAGCAATGAc       >  2:276215/1‑143 (MQ=255)
                                                                                                                                       cGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATCTCAATTAcccc    >  1:275463/1‑141 (MQ=255)
                                                                                                                                        gggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGACTGAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCGGCGCTCTGTCGGCCACGACTTCGTCGGCTTCAATCTCAATaagccca   >  2:707045/1‑138 (MQ=255)
                                                                                                                                         ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCCACGTCTTCGTCGGCCTCAATATCGATAACcacac  >  2:576020/1‑142 (MQ=255)
                                                                                                                                         ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCAAATGGTACCCCGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACAt  <  1:777365/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
CCCCATGCTGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACAT  >  NC_000913/3056999‑3057278

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑