Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,609,336 (C)7→9 intergenic (+13/‑78) yjjG → / → prfC pyrimidine nucleotidase/peptide chain release factor RF3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,609,3291.C100.0% 148.3 / NA 43intergenic (+6/‑85)yjjG/prfCpyrimidine nucleotidase/peptide chain release factor RF3
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (28/15);  total (28/15)
*NC_0009134,609,3292.C95.3% 142.0 / 0.1 43intergenic (+6/‑85)yjjG/prfCpyrimidine nucleotidase/peptide chain release factor RF3
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (27/14);  total (28/15)
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.80e-01

TTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCT‑‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACT  >  NC_000913/4609190‑4609458
                                                                                                                                              ||                                                                                                                                  
ttCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCt                                                                                                                                      <  1:1231673/140‑1 (MQ=255)
  ctcGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCT‑‑‑cc                                                                                                                                 <  1:452119/140‑1 (MQ=255)
     ggTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCT‑‑‑ccccc                                                                                                                              <  2:764794/140‑1 (MQ=255)
          cATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCT‑‑‑cccccccccg                                                                                                                         >  1:169661/1‑137 (MQ=255)
            tCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCT‑‑‑cccccccccgtt                                                                                                                       <  2:169661/140‑6 (MQ=255)
             cAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCT‑‑‑cccccccccgttg                                                                                                                      <  1:663396/140‑7 (MQ=255)
             cAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCACGAGCAACCAGAAGGCATCGCGCCCACCTGGTCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCT‑‑‑cccccccccgttg                                                                                                                      <  1:61205/140‑7 (MQ=255)
              caCGCCGGGCTTGCGGCCTGCTGGCTGAAAGCACAACCTCGGGAGCAACCAGAAGGGATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGGCCCC‑‑‑CCCCCCCGTTGa                                                                                                                       <  2:448389/139‑1 (MQ=255)
                     ggCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGTTGGGTaa                                                                                                                >  1:448389/1‑140 (MQ=255)
                                  ctggctgAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAtt                                                                                                   >  2:1225298/1‑140 (MQ=255)
                                         aaTGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtt                                                                                            >  1:568204/1‑140 (MQ=255)
                                               cacCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                                                                      >  2:765156/1‑140 (MQ=255)
                                                   aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGt                                                                                             >  1:1037994/1‑129 (MQ=255)
                                                   aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGt                                                                                             <  2:1037994/129‑1 (MQ=255)
                                                   aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgc                                                                                  >  2:439209/1‑140 (MQ=255)
                                                   aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgc                                                                                  >  2:1237441/1‑140 (MQ=255)
                                                    tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcg                                                                                 >  1:1337922/1‑140 (MQ=255)
                                                     cgcgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAtttt                                                                                                 <  2:206025/123‑1 (MQ=255)
                                                     cgcgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAtttt                                                                                                 >  1:206025/1‑123 (MQ=255)
                                                         aGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcg                                                                            >  1:301223/1‑140 (MQ=255)
                                                          gCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcgc                                                                           >  2:1033329/1‑140 (MQ=255)
                                                            aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                                                         <  2:568204/140‑1 (MQ=255)
                                                                   aGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                  <  1:564399/140‑1 (MQ=255)
                                                                   aGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                                                                  <  2:747756/140‑1 (MQ=255)
                                                                      cATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTC‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACt                                                              >  2:1131724/1‑140 (MQ=255)
                                                                         cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCa                                                            >  1:1227612/1‑140 (MQ=255)
                                                                               cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga                                                      <  1:1225298/140‑1 (MQ=255)
                                                                                aCCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAgaagaa                                                     >  1:798584/1‑140 (MQ=255)
                                                                                     gACCGTTTCTTCGTTGCACGAAGTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAagaaagaa                                                 >  1:131023/1‑140 (MQ=255)
                                                                                          tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                           >  1:668881/1‑124 (MQ=255)
                                                                                          tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                           >  2:1226412/1‑140 (MQ=255)
                                                                                          tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                           >  2:1260456/1‑140 (MQ=255)
                                                                                          tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAAAACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                           <  2:668881/124‑1 (MQ=255)
                                                                                           ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTC‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                         <  1:1131724/140‑1 (MQ=255)
                                                                                            tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                         >  2:391158/1‑140 (MQ=255)
                                                                                            tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                         <  1:1033329/140‑1 (MQ=255)
                                                                                            tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACg                                         <  1:1381894/140‑1 (MQ=255)
                                                                                                cGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGt                                     <  2:798584/140‑1 (MQ=255)
                                                                                                    gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTcc                                 >  1:897446/1‑140 (MQ=255)
                                                                                                    gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAACAAAGAATTATGACGTTGTCTcc                                 <  1:765156/140‑1 (MQ=255)
                                                                                                     cACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCt                                >  2:100226/1‑140 (MQ=255)
                                                                                                     cACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCt                                >  1:637958/1‑140 (MQ=255)
                                                                                                      aCGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCtt                               <  2:1337922/140‑1 (MQ=255)
                                                                                                          aCTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTAttt                           >  1:36672/1‑140 (MQ=255)
                                                                                                             ggAGCAGCTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCa                        >  1:986461/1‑140 (MQ=255)
                                                                                                                    cTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg                 >  2:944766/1‑140 (MQ=255)
                                                                                                                    cTCCTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTg                 >  1:454583/1‑140 (MQ=255)
                                                                                                                       cTGTGTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCg              <  1:439209/140‑1 (MQ=255)
                                                                                                                        tgtgTAAACACTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGa             <  2:454583/140‑1 (MQ=255)
                                                                                                                                 acTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCa    >  2:66296/1‑140 (MQ=255)
                                                                                                                                 acTGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCa    >  1:156572/1‑140 (MQ=255)
                                                                                                                                   tGATTG‑‑CCTCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACt  <  1:812467/140‑1 (MQ=255)
                                                                                                                                              ||                                                                                                                                  
TTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTG‑‑CCT‑‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACT  >  NC_000913/4609190‑4609458

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