Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,609,336 (C)7→8 intergenic (+13/‑78) yjjG → / → prfC dUMP phosphatase/peptide chain release factor RF‑3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,609,3291.C91.8% 153.5 / 6.2 49intergenic (+6/‑85)yjjG/prfCdUMP phosphatase/peptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (23/22);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

TCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCG  >  NC_000913/4609193‑4609454
                                                                                                                                          |                                                                                                                              
tcGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTcc                                                                                                                              >  2:1092053/1‑139 (MQ=255)
         tCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgtt                                                                                                                    >  2:1697599/1‑135 (MQ=255)
         tCAACGCCGGGCTAGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgtt                                                                                                                    >  1:1379263/1‑135 (MQ=255)
                                      aaTGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtt                                                                                        >  1:755620/1‑139 (MQ=255)
                                            cacCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgttg                                                                                                                   <  1:1724404/105‑6 (MQ=255)
                                            cacCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgttg                                                                                                                   >  2:1724404/1‑100 (MQ=255)
                                                    cgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAATACTGATTGCCTccccccc                                                                                                                        <  1:1142250/92‑1 (MQ=255)
                                                    cgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccc                                                                                                                        >  2:1142250/1‑92 (MQ=255)
                                                         aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                                                     >  1:1470317/1‑139 (MQ=255)
                                                           ccAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATg                                                                   >  2:1318708/1‑139 (MQ=255)
                                                                  gCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTa                                                            >  2:494416/1‑139 (MQ=255)
                                                                            cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCTGCTCCTGTGTAAACACTGATTGCCGCCCCCCCCGTTGATGGGTAAAATAGCCGCAAGTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  <  1:1582941/139‑1 (MQ=255)
                                                                            cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  <  1:1092053/139‑1 (MQ=255)
                                                                            cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  <  2:1470317/139‑1 (MQ=255)
                                                                              ccTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGaagaaa                                                <  1:1318708/139‑1 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                        <  2:1499415/106‑1 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                        >  1:1499415/1‑106 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                        <  2:874460/106‑1 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                        >  1:874460/1‑106 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggc                                                                         <  2:524020/105‑1 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggc                                                                         >  1:524020/1‑105 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgttg                                                                                                                   >  2:2012394/1‑57 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgttg                                                                                                                   >  1:941787/1‑57 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgttg                                                                                                                   <  2:941787/62‑6 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgttg                                                                                                                   <  1:2012394/62‑6 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccc                                                                                                                        >  1:797295/1‑57 (MQ=255)
                                                                                       tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccc                                                                                                                        <  2:797295/57‑1 (MQ=255)
                                                                                        ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATg                                                                                                                 >  1:445438/1‑64 (MQ=255)
                                                                                        ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATg                                                                                                                 <  2:445438/64‑1 (MQ=255)
                                                                                        ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  >  1:1035066/1‑127 (MQ=255)
                                                                                        ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  <  2:1035066/127‑1 (MQ=255)
                                                                                        ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccc                                                                                                                          >  2:220981/1‑54 (MQ=255)
                                                                                        ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccc                                                                                                                          <  1:220981/54‑1 (MQ=255)
                                                                                            tcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTg                                  >  1:1737205/1‑139 (MQ=255)
                                                                                             cgctGCACCAAATGGATCAGCTCCTTTGGTATCCCTGGTTGCCTCCCCCCCCGTTGAAGGGTAAAATATGCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGt                                 <  1:329867/136‑1 (MQ=255)
                                                                                               ttGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                    <  2:1452394/102‑1 (MQ=255)
                                                                                               ttGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgttg                                                                                                                   >  2:465149/1‑49 (MQ=255)
                                                                                               ttGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTccccccccgttg                                                                                                                   <  1:465149/54‑6 (MQ=255)
                                                                                               ttGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCGCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                    >  1:1452394/1‑102 (MQ=255)
                                                                                                    cGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTCGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTa                          <  1:944845/139‑1 (MQ=255)
                                                                                                       aCTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGTGATGCCGCTTACTCACGAAGAAAGAATTATGACGTTGTCTCCTTAttt                       >  2:329867/1‑139 (MQ=255)
                                                                                                           gAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                     <  2:644479/121‑1 (MQ=255)
                                                                                                           gAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                     >  1:644479/1‑121 (MQ=255)
                                                                                                                 cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg             >  2:327083/1‑139 (MQ=255)
                                                                                                                 cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg             >  1:1799048/1‑139 (MQ=255)
                                                                                                                 cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg             >  1:299201/1‑139 (MQ=255)
                                                                                                                     tgtgTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGaa        >  2:489529/1‑139 (MQ=255)
                                                                                                                         tAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAgcg     <  1:327083/139‑1 (MQ=255)
                                                                                                                          aaaCACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGcgcc   <  1:489529/139‑1 (MQ=255)
                                                                                                                           aaCACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAgaagaa                                                 >  2:2022909/1‑93 (MQ=255)
                                                                                                                           aaCACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAgaagaa                                                 <  1:2022909/93‑1 (MQ=255)
                                                                                                                            acacTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGcgccg  <  2:299201/139‑1 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACt                                                          <  2:1534982/78‑1 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACt                                                          >  1:1534982/1‑78 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATg                                        <  1:46754/96‑1 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATg                                        >  2:46754/1‑96 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGTCGtt                                   <  1:855639/101‑1 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                     <  2:1264460/99‑1 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                     >  1:1264460/1‑99 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGtt                                   >  2:855639/1‑101 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTa                          >  1:848250/1‑110 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTa                          <  2:848250/110‑1 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTAt                         <  2:1878256/111‑1 (MQ=255)
                                                                                                                                 gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTAt                         >  1:1878256/1‑111 (MQ=255)
                                                                                                                                  aTTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  <  1:1658704/85‑1 (MQ=255)
                                                                                                                                  aTTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                  >  2:1658704/1‑85 (MQ=255)
                                                                                                                                          |                                                                                                                              
TCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCG  >  NC_000913/4609193‑4609454

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