Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 1,376,811 C→T 10.5% R260C (CGT→TGT)  ycjR → putative TIM alpha/beta barrel enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,376,8110CT10.5% 82.4 / 2.3 36R260C (CGT→TGT) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base C (15/17);  new base T (2/2);  total (17/19)
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.56e-01

GCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGA  >  NC_000913_3_bme_pgi/1376665‑1376956
                                                                                                                                                  |                                                                                                                                                 
gCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGCCTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGt                                                                                                                                                 >  1:663447/1‑149 (MQ=255)
 cTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTa                                                                                                                                                >  1:605200/1‑149 (MQ=255)
    ccAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGTGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCt                                                                                                                                             >  1:37225/1‑149 (MQ=255)
      aGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGc                                                                                                                                           >  2:258013/1‑149 (MQ=255)
           ggCTGCGGGACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGTCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTTCGTACCTGCTAaga                                                                                                                                      <  2:209316/149‑1 (MQ=255)
               gCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGACGAAAGATCCTGCAAAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTc                                                                                                                                  >  1:74512/1‑149 (MQ=255)
                  gccacccggATTTTCCCGCTGTGTTTGAACCGCCGAGCGCGGATAACTATCAGGGCTATGTCGTGTACGAAGGGCGTATCCTGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCttt                                                                                                                               <  2:136407/142‑1 (MQ=255)
                     cccTGGATTTCAACGCGCTGTTTGAACAGCTGCGCGCGGATGACTATCAGGGCTATGTAGGGTACGAAGGGCGTATCCGGGCGCAAGATCCTGCACAGGCGTACCGTGCTTCGTTGGCCGGGTTGCGAACCTGCGAAAAGGTGtttaga                                                                                                                            >  2:10487/1‑146 (MQ=255)
                         tgATTTCCACGCGATGTTTGAACAGCTGAGCGAGTATAAATATCAAGGCTATGTAGTGTACGAAGGGCGTATCCGGGAGGAAGAGCCTGACCAGGGGTACCGTGATTCGTTGGCCGTGTTGCGTACCTGCTAAGAGGTCTTTGTGaaaa                                                                                                                        <  1:390191/148‑1 (MQ=255)
                         ggATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGAGTATTCGGGAGGAAGATCCTGCACAGGCGTACAGTGATTCGTTGGCCTGGTTCCGTACCTGCTAAGAGGTATTTGTAaaaa                                                                                                                        >  1:373338/1‑149 (MQ=255)
                               ccACGCGCTGTTTGAACAGCTGCGGCCTGATAACTATCAGGGATATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCaa                                                                                                                  <  1:743431/149‑1 (MQ=255)
                                   gcgcTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGGCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGAc                                                                                                              <  1:685213/149‑1 (MQ=255)
                                     gcttgttGAACAGCTGCGCGCGGATAACTATTAGGGTTATTTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACaa                                                                                                            <  1:1037/144‑1 (MQ=255)
                                      cTGTTTGAACAGCTGCGCGCGGATAGCTATGAGGGCTGTGTAGTATATGACGGGGGAGGGCGGGGGGAAGATCCTGCCAAGGCGTACCGTGATTGGTTGGCCTGGTTGTGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAg                                                                                                           <  2:638611/149‑1 (MQ=255)
                                        tgttGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAATATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGct                                                                                                         <  2:779286/147‑1 (MQ=255)
                                           tGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATTCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCtctc                                                                                                      <  2:682323/149‑1 (MQ=255)
                                                           gATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCCCCGCTGCGGGTCGCGat                                                                                      >  1:791492/1‑149 (MQ=255)
                                                             tAACTATCAGGGCTTTGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCTCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGataa                                                                                    <  1:765081/149‑1 (MQ=255)
                                                              aaCTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTAGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGAGatcat                                                                                   >  1:239577/1‑149 (MQ=255)
                                                               aatagcAGGGCTATGTAGGGTACGGAGTGCGTATCCGGGCGGAAGATCTGGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGataata                                                                                  <  2:175794/144‑1 (MQ=255)
                                                                          tATGTAGTGTACGAAGGGCTTATCCGGTCGGAAGATCCTGCCCAGGCGTACCGGGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCa                                                                       <  1:124777/149‑1 (MQ=255)
                                                                                gtctccGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGTCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGc                                                                 <  1:464028/144‑1 (MQ=255)
                                                                                           gCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGTGGGTCGCGATACTAGGCGCAGGCCAGGTGGCGGAGAAGGTTc                                                      >  2:204214/1‑149 (MQ=255)
                                                                                           gCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGATGGCGGATAAGGTTc                                                      >  1:302826/1‑149 (MQ=255)
                                                                                                             tCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGta                                            >  1:413478/1‑141 (MQ=255)
                                                                                                             tCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGta                                            <  2:413478/141‑1 (MQ=255)
                                                                                                                             cGTGATTCGTTGGCCTGGTTGTGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCACTGCCAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGa                    >  1:576824/1‑149 (MQ=255)
                                                                                                                              gTGGTTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGaa                   <  2:37225/149‑1 (MQ=255)
                                                                                                                              gTGATTTGTTGGACTGGTTGTGTGGCGGCTGAGAGGTCGTTATGAAAGTTGGGATGACAAGCGGTGTTCTGGATGTGGCGGTAATAGTGGCGGGTGTGGTGGCGGATAGCGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGaa                   <  2:753098/149‑1 (MQ=255)
                                                                                                                               tGATTCGTTGGCCTGGTTGTGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCACTGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCGTCGTACTACTGCACCCGCAACGATCTGCAAc                  >  1:192612/1‑149 (MQ=255)
                                                                                                                                 aTTCGTTGGCCTGGTTGCGTACCTGCTAGAAGGGCCTTGTGAAACGTACAATGAAAAGGTCCGCGCTGGGGGTCGCGAAAATATGCGCAGGCAAGATGGCGGATAAGGTTCATGCTTAGTACGCCAGCAAACGCAGCAATCTGGAACtg                >  2:534722/1‑149 (MQ=255)
                                                                                                                                 aTTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTg                                      <  1:338145/127‑1 (MQ=255)
                                                                                                                                 aTTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTg                                      >  2:338145/1‑127 (MQ=255)
                                                                                                                                  ttCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCTAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACtgg               >  1:898693/1‑149 (MQ=255)
                                                                                                                                           cTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAATTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGctgtct      <  1:258013/149‑1 (MQ=255)
                                                                                                                                               ttGCGTACTGTCTGAGAGGTGTTTTTGAAAAGTGCAATGACATGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGCTGGCAGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGa  <  2:668056/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGA  >  NC_000913_3_bme_pgi/1376665‑1376956

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: ‑