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.C90.9% 102.8 / 2.8 33intergenic (+6/‑85)yjjG/prfCdUMP phosphatase/peptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (17/13);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTT  >  NC_000913/4609192‑4609459
                                                                                                                                            |                                                                                                                                    
ctcGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑c                                                                                                                                    <  1:719859/139‑1 (MQ=255)
          tCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑cccccccccgt                                                                                                                         >  2:595810/1‑135 (MQ=255)
           cAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                        <  1:462974/139‑6 (MQ=255)
           cAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                        <  2:1417018/139‑6 (MQ=255)
                 cGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGt                                                                                                                   >  1:130093/1‑139 (MQ=255)
                 cGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                        <  2:1227842/133‑6 (MQ=255)
                 cGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                        >  1:1227842/1‑128 (MQ=255)
                      ttGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAAt                                                                                                              <  1:1588293/139‑1 (MQ=255)
                       tGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCCGTTGATGGGTAAAAt                                                                                                              <  1:595810/139‑1 (MQ=255)
                         cGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGc                                                                                                           >  1:203110/1‑139 (MQ=255)
                                                cATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAg                                                                                    <  2:130093/139‑1 (MQ=255)
                                                 aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTTGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                        >  2:356418/1‑96 (MQ=255)
                                                 aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                        <  1:356418/101‑6 (MQ=255)
                                                  tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑cccc                                                                                                                                 >  1:68517/1‑92 (MQ=255)
                                                  tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑cccc                                                                                                                                 <  2:68517/92‑1 (MQ=255)
                                                  tCGCGAGCAACCAGAAGGCATCGCGACCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCCCAATTTTTCGTTTTCAACAAgcg                                                                                  >  1:1414416/1‑139 (MQ=255)
                                                          aaCCAGAAGGCATCGCGCCCACCTGGAGCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                                                                       >  1:1121917/1‑126 (MQ=255)
                                                          aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTTCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                                                                       <  2:1121917/126‑1 (MQ=255)
                                                                                    actgTTTCTTCGTTGCGAGAGCTGGGGCTGCTCCGGTGTAAATACGGTTGGCCAC‑‑CCCCCCCGTTGATGGGTAAAATAGCTGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCATACTCAAGAAGAAAGAAtt                                                <  2:203110/136‑1 (MQ=255)
                                                                                        ttttttCGTTGCACGAACTGGAGCAGCTCCTGTTTAAACACTGATTGCCCC‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                                                                 <  2:807468/86‑1 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCC‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                                                                 >  1:807468/1‑86 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCa                                                                                                       <  2:190811/80‑1 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCa                                                                                                       >  1:190811/1‑80 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTTCTCAAGAAGAAAGAATTATGa                                            >  1:400783/1‑139 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                            <  1:463450/123‑1 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                            >  2:463450/1‑123 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                            >  2:892620/1‑139 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGa                                            >  1:965259/1‑139 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccc                                                                                                                                  <  1:281049/53‑1 (MQ=255)
                                                                                        tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccc                                                                                                                                  >  2:281049/1‑53 (MQ=255)
                                                                                          tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                        <  2:1559331/60‑6 (MQ=255)
                                                                                          tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑ccccccccgttg                                                                                                                        >  1:1559331/1‑55 (MQ=255)
                                                                                            ttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGtt                                        <  2:1414416/139‑1 (MQ=255)
                                                                                                   cACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                         <  1:259123/99‑1 (MQ=255)
                                                                                                   cACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                         >  2:259123/1‑99 (MQ=255)
                                                                                                   cACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCt                                 >  2:615788/1‑139 (MQ=255)
                                                                                                             agcagcTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAag                       <  1:892620/139‑1 (MQ=255)
                                                                                                                  cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                         >  1:137843/1‑84 (MQ=255)
                                                                                                                  cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                            >  2:1646376/1‑97 (MQ=255)
                                                                                                                  cTCCTGTGTAAACACTGATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                                                            <  1:1646376/97‑1 (MQ=255)
                                                                                                                  cTCCTGTGTAAACACTGATTGCCT‑‑CCCCCCTCGTTGATGGGTAAAATAGCCGCAACTTTTCGTTTTCAACAAGCGCGGCGCGAt                                                                         <  2:137843/84‑1 (MQ=255)
                                                                                                                                  gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTAcg                                                               >  2:1097978/1‑77 (MQ=255)
                                                                                                                                  gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACt                                                               <  1:1097978/78‑1 (MQ=255)
                                                                                                                                  gATTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACtt  >  1:547632/1‑139 (MQ=255)
                                                                                                                                   aTTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAGTTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTc                                                              >  1:1158656/1‑78 (MQ=255)
                                                                                                                                   aTTGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCCCGATGCCGCTTACTc                                                              <  2:1158656/78‑1 (MQ=255)
                                                                                                                                    ttGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGcc                                                                      >  1:1599434/1‑69 (MQ=255)
                                                                                                                                    ttGCCTCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGcc                                                                      <  2:1599434/69‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                    
CTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCT‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGACGTTGTCTCCTTATTTGCAAGAGGTGGCGAAGCGCCGCACTT  >  NC_000913/4609192‑4609459

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