Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 147,213 2 bp→TT coding (246‑247/903 nt) rpnC → recombination‑promoting nuclease RpnC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913147,2130CT83.9% 68.9 / 8.1 32H82H (CAC→CATrpnCrecombination‑promoting nuclease RpnC
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (9/17);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01
*NC_000913147,2140CT84.4% 65.1 / 7.3 32Q83* (CAA→TAA) rpnCrecombination‑promoting nuclease RpnC
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (9/18);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

GGATGCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCA  >  NC_000913/146970‑147445
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                        
ggatgCACCGAGTAACACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAaagca                                                                                                                                                                                                                                    <  2:102242/251‑1 (MQ=255)
  atgCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAaagcaag                                                                                                                                                                                                                                  <  2:133713/251‑1 (MQ=255)
              ccACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGaaa                                                                                                                                                                                                                      <  1:383911/251‑1 (MQ=255)
                  acCGCTTGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACGTTCTGGAGATACATTTTCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCGTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCAGTTAACTTCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATgg                                                                                                                                                                                                                  <  2:211285/251‑1 (MQ=255)
                                                aaTGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCCAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAgccg                                                                                                                                                                                    >  1:317655/1‑251 (MQ=255)
                                                aaTGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAgccg                                                                                                                                                                                    >  1:249887/1‑251 (MQ=255)
                                                   gcatgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCa                                                                                                                                                                                 <  2:81209/251‑1 (MQ=255)
                                                     atgcGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGAAACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATg                                                                                                                                                                               >  1:46890/1‑251 (MQ=255)
                                                         ggAGACGGCTCGCGAATTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCAcc                                                                                                                                                                           <  2:231210/251‑1 (MQ=255)
                                                                      gACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGc                                                                                                                                                              <  1:296166/251‑1 (MQ=255)
                                                                          ttCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTg                                                                                                                                                            <  1:300383/249‑1 (MQ=255)
                                                                              ggAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCAc                                                                                                                                                       <  1:26563/250‑1 (MQ=255)
                                                                                                 gAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTTAAGAGAGACTGAAAGGACACAACACGGAAGTGCTCTATTCCGGGCAAAAGCCGGGCCAACCCGGTTATCTGCATGTTGTGATTGAACATTAAAGACAGCCGGATAAGAAAATGGCCTTTCGCTTGATGCGTTATTATAAAACCCCCATGCACCGGCATCTGGAGGCTgccaccgatgaggcgccgcaggtggg                                                                                                                                   >  2:437793/1‑226 (MQ=255)
                                                                                                      aCGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCg                                                                                                                               >  1:387580/1‑250 (MQ=255)
                                                                                                             cTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGGTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGTAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGtt                                                                                                                       >  1:115179/1‑251 (MQ=255)
                                                                                                                        cAACACGCTTCATTTAGAGTCTGGGAGTTTCATTGAAGAGATCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCGGCATGTTGTGATTTAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCg                                                                                                            <  1:313302/251‑1 (MQ=255)
                                                                                                                        cAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCg                                                                                                            <  1:17222/251‑1 (MQ=255)
                                                                                                                                   aTTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCAGTGCAAATGCAGGGCCATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGATGGCCACGGTAAGCTGCCGCAGGTGGTGCCGCAACTGTTTTATCAGGGCGAGGCcaccccc                                                                                                 >  2:129510/1‑247 (MQ=255)
                                                                                                                                      tAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCAGCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACCCCTTa                                                                                               <  1:268420/250‑1 (MQ=255)
                                                                                                                                                              agagCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGat                                                                     <  1:5755/251‑1 (MQ=255)
                                                                                                                                                                          gACACAGCACGAACGTGCTCTATTACGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGCACATTAAAGCAAGCCGGATAAGAAAAAGGCCTTTAGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGTCAACGCTAAACTTCCTCTGGTGGGGCCGATACTGTTTTATCAGGGCTAGTCCCCACCTTATCCGGTATCCATTGTATGGTTTGTTATGTTTTCCTCg                                                         >  2:118681/1‑251 (MQ=255)
                                                                                                                                                                                ataaggtcGTGCTCTATTACGTGCAAATTAAGGTCATTCCCGGTTTACTGCATGGTTTGATTGAACATTAAAGCACTCCGGATAGGAAAATTGCAATTCGCAGGATGAGTGAGTCAATAGCCGCCATGCACCGTCCTCGGGAGGCTGACTACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTGTATCAGGGCGAGGCCACACTTTATCCGCTATCAATTTGCTGGTTTGATATGTTTTACTCGCCGGAg                                                   <  2:96152/244‑1 (MQ=255)
                                                                                                                                                                                   cGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTg                                                <  1:129510/251‑1 (MQ=255)
                                                                                                                                                                                       cGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGTTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCTGATAAGAAAATGGCCTTTCGCATTATGCGTTATTCTATAGCCGCCATGCTCCGGCATCTGGAGTCTTACCACGATAATCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGGCACACCTTATCCGCTATCAATGTGCTGGtt                                                                         <  1:411912/222‑1 (MQ=255)
                                                                                                                                                                                       cGTGCTCTATCCCGTCCAAATGCAGGACAATCCCGGTTCTCTGCATGTTGTGATTGAACACCACAGCAAGCCGGATAAGAAAATGGCCTTTCGCCTGATGCGTTATTCTATAGCCGCAATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGCTACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGtt                                                                         >  2:411912/1‑222 (MQ=255)
                                                                                                                                                                                        gggCTCTATTCCGTGCAAATGCAGGGCAATTCCGGTTATCTGCAGCTTGGGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTGCGCATGATGGGTTATTCTATAGCCGCCATGCACCTGCAGCTTGAGGCTGACCACGATATGCTGCCGCTGGTTGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCTCCGGAGCTGGcgcg                                           <  2:374816/249‑1 (MQ=255)
                                                                                                                                                                                        gTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGcgcg                                           >  1:312744/1‑251 (MQ=255)
                                                                                                                                                                                               aTTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTc                                    <  1:266391/251‑1 (MQ=255)
                                                                                                                                                                                                 cccctacactacgctcgtccgatcGGTTATCTGCATGTTTGGATGGAACAGTAAAGCAAGCCGGATAAAAAAATGGCTTTTCTCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACGACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGTGAGGCCACGCCTTATCCGCTATCAATGTGCTGGTTGGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCt                                   <  2:220779/227‑1 (MQ=255)
                                                                                                                                                                                                                          gTTATCTGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCt                                   >  1:220779/1‑225 (MQ=255)
                                                                                                                                                                                                                                tGCATGTTGTGATTGAACATTAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGAGTTATTCTATAGCCGACATGCACCGGCATCTGGAGACTGACCACGATAAGCTGCCGCTGGTGGTGCCGATAATGTTTTATCAGGGCGAGGCAACAAATTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCAcc   >  1:49471/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCATGTTGTGATTAAACATTAAAGCAAGCGGGATAAGAAAAGGGCCTGTCGCATGATGCGTTATGCTAGAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAATCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCa  <  2:432929/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                        
GGATGCACCGAGTACCACACCGCATGACGCGGTATTTAAACAATTTTTAATGCATGCGGAGACGGCTCGCGACTTTCTGGAGATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAAGAAAATGGCCTTTCGCATGATGCGTTATTCTATAGCCGCCATGCACCGGCATCTGGAGGCTGACCACGATAAGCTGCCGCTGGTGGTGCCGATACTGTTTTATCAGGGCGAGGCCACACCTTATCCGCTATCAATGTGCTGGTTTGATATGTTTTACTCGCCGGAGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCA  >  NC_000913/146970‑147445

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