Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,256,301 G→A S117S (TCC→TCT cyuA ← putative L‑cysteine desulfidase CyuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,3010GA93.3% 73.1 / ‑2.1 30S117S (TCC→TCTcyuAputative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (12/16);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCA  >  NC_000913/3256055‑3256508
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                               
cctcaacgggacgccatttagaatcTTCAGGTTCTCATTCAGCGTCCGTTGTGAAATCAACGTCAGTGGAGACTCTTGCTCGCCCTCTTCCTCACACGCTTCCTGGGTAATCACCTCCCCATCGTGCGTCTTGATATGCACAATGTTGGTATGCCCGCCGACTATGGTGACACACGCCCACTTCTCACCGTTCCATACTTTGGCGCGTGAGAAGAGGATTTCATCACAAGGTTCCTGGATCTTAACAgaga                                                                                                                                                                                                             <  2:323537/227‑1 (MQ=255)
      cGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTcc                                                                                                                                                                                                       <  1:179620/251‑1 (MQ=255)
             tCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACAGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCcgccgc                                                                                                                                                                                                >  1:349068/1‑251 (MQ=255)
               aTTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCgccgccc                                                                                                                                                                                              >  1:124025/1‑250 (MQ=255)
                           cAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTTTTGCTCGCCCTCTGCCACCCACGCCTGCTGGGTAAACACCCCCCCACCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTgg                                                                                                                                                                                  <  1:83954/251‑1 (MQ=255)
                             gtaTCTCAGCCCGCGTCGTTCTGGAAAGAACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACCCCTGCTGGGTAAACACCACACCATCGTGCGTCCCGATATGCAACATGTTGGTATGCACACCGCCGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCACAGTTACTGAATCTTAACAGAGACTTTCCACGCCGCCAGCAGGACTTTGGc                                                                                                                                                                                 >  2:218944/3‑250 (MQ=255)
                             ggATCTCAGCCCGCTTCGTTCCGGAAAGCACCGTCCGCGGAGTCTCTTGCTCGCCTCCTGCCGCACCCCCCTGGTGGGGAAACCCCACCCCAACGGGCTTCACAATTGGCACAGTGTGGGTATGCCCGCCGACGATGGTGACACACGCCCCCTTCTCACCGGTCCAGACTTTGGCGCGGGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTGCCCCGCCGCCAGCAGTGCTTTGGCa                                                                                                                                                                                <  2:383646/251‑1 (MQ=255)
                                cctcAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGATACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTGGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCg                                                                                                                                                                             <  1:309243/250‑1 (MQ=255)
                                 ctcAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCgg                                                                                                                                                                            >  1:165763/1‑251 (MQ=255)
                                   cagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGc                                                                                                                                                                           >  1:45622/1‑250 (MQ=255)
                                      ccagcGTGGTTCTGGAAAGCACCGCCAGCGGAGCCTCTTGCTCGCCCTCGGCCACACACGCCTGCTGGGTAAACCCCCCCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAAtt                                                                                                                                                                       <  1:115011/251‑1 (MQ=255)
                                      ccagcGTCGATCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACACTGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCCCCCCCAGCAGTGCTTTGGCATCGGCAAtt                                                                                                                                                                       >  2:327474/1‑251 (MQ=255)
                                                                     gataCTCTTGCGCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                                                                                                                                         <  1:384024/248‑1 (MQ=255)
                                                                       gaCTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                                                                                                                                         >  1:225288/1‑249 (MQ=255)
                                                                               cacGCCCTCGGCCCCACCCGCCTGCTGGGTAAACACCACTCCATCGTGCTTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGt                                                                                                                               <  1:100407/249‑1 (MQ=255)
                                                                                    cccctGCCACACACGCCTGCGTGTTAAACCCCACACCAGCGTGCGTCTCTATATGCAAAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGATTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCa                                                                                                                          <  1:147131/248‑1 (MQ=255)
                                                                                                   ccTGCTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCTTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGcc                                                                                                           <  1:22925/251‑1 (MQ=255)
                                                                                                       ctgGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGccccc                                                                                                        >  1:362491/1‑250 (MQ=255)
                                                                                                            aaaccaccactccatcgggctgcacGATATGCACAAGGTTGGTAGGCCCGCCGACGATGGGGACACACGCCCACTTCCCACCGTCCCAGACTTTGGCGCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGcgcc                                                                                                  <  1:318580/227‑1 (MQ=255)
                                                                                                               acacctccccaTCGTGCGTCTCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGcgccgcc                                                                                               <  1:207110/243‑1 (MQ=255)
                                                                                                                cacccccccagcGTGCGTCTCGTTATGCACAATGTGGGTAGGCCCGCCGACGATGGTGCCACACGCCCACTGCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGcgccgccg                                                                                              <  1:141892/240‑1 (MQ=255)
                                                                                                                        caTCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGc                                                                                      >  1:128450/1‑251 (MQ=255)
                                                                                                                                                gtttgttatccccGCCGACGATGGTGACACCCCCCCACTCCTCACCTTTCCAGCCTTTGGCGCGGGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCggg                                                               <  1:249391/240‑1 (MQ=255)
                                                                                                                                                       atgcccccccaggaTGGTGACCCCCCCCCACTCCCCCCCGTTCCAGTCTTTGGCGCGTGATAAGAGTATTTCATCCCAGGGTTCCTGGATTTAAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCACCGGCAATGGCCTGCGCTGTTGCGTCTTTCAGCCCTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTg                                                        <  2:165763/238‑1 (MQ=255)
                                                                                                                                                        gccccGACGACGGTGGGGACACACGCCCTCTTCTCACCGTTCCAGACGTTGGCGCGGGAGAAGAGGATTCCATCGCAAGGGTCCTGGATCTAAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGAc                                                      <  1:249406/249‑1 (MQ=255)
                                                                                                                                                            cccctagtgtggtacccctcGCCCCCTCCTCACCGTCCAGAGTTTGGGCGCGTGAGAAGAGGATTCCATAGCAAGGTTCCGGGTTCATAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTGGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCCCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGccc                                                  <  2:191776/232‑1 (MQ=255)
                                                                                                                                                                   gaTGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCATGGATCTTAACAGAGACTTTCCCCGCCGCCAGAAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTCGCATTTCCAACTAACGCCCCCAGCGCCGCCGCACTCGGCAGCCCCACCATTCCCGTGCCGGGAACGGAGCCGCCAAGACCgc                                           >  1:252185/1‑250 (MQ=255)
                                                                                                                                                                                   cACTTCTCACCGTTCCAGACTTGGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTGCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTgg                           <  1:218944/251‑1 (MQ=255)
                                                                                                                                                                                                 ccAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGACGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCtt             >  1:11429/1‑251 (MQ=255)
                                                                                                                                                                                                              gcgTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCACCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGCAACGGCGACGCCCAGACCGATCTTCCTCAGATTTGCAGAAACCCAGGCTTCTacacttaca  >  2:271787/1‑245 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                               
CGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCA  >  NC_000913/3256055‑3256508

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: