Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,150 C→T 9.6% intergenic (+159/+28) 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,1500CT 9.6% 95.8 / 3.3 42intergenic (+159/+28)sibC/serAsmall RNA SibC/phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (23/15);  new base T (2/2);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

CTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACAT  >  NC_000913/3057011‑3057278
                                                                                                                                           |                                                                                                                                
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)
                                                                                                                                           |                                                                                                                                
CTGCCATCGTAAATCCCCATTAAATAAACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTCTTCGTCGGCTTCAATATCAATAACCACAT  >  NC_000913/3057011‑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: ‑