Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,145 G→C 9.4% intergenic (+154/+33) 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,1450GC 9.4% 117.2 / 6.2 43intergenic (+154/+33)sibC/serAsmall RNA SibC/phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (23/16);  new base C (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 = 5.47e-01

TGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACAT  >  NC_000913/3057007‑3057278
                                                                                                                                          |                                                                                                                                     
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)
                                                                                                                                          |                                                                                                                                     
TGCCCTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACAT  >  NC_000913/3057007‑3057278

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