Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,445 (C)7→8 intergenic (+90/‑265) rnlB → / → yfjP CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,4381.C89.2% 110.1 / 5.9 37intergenic (+83/‑272)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (12/21);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 2.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGG  >  NC_000913/2767301‑2767571
                                                                                                                                           |                                                                                                                                         
aGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTaag                                                                                                                                           >  2:312355/1‑138 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTaaaac                                                                                                                                         >  1:379811/1‑136 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                       <  1:1057442/139‑1 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                       <  1:992939/139‑1 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccc                                                                                                                                    >  2:971657/1‑139 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccc                                                                                                                                    >  2:594982/1‑139 (MQ=255)
      tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccccc                                                                                                                                  >  1:599102/1‑139 (MQ=255)
                 ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                       <  1:594982/139‑1 (MQ=255)
                   gatgatGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCt                                                                                                                     <  1:971657/139‑1 (MQ=255)
                              ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                          <  1:1072611/139‑1 (MQ=255)
                                  gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccccctca                                                                                                                              <  2:310878/115‑5 (MQ=255)
                                  gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccccctca                                                                                                                              >  1:310878/1‑111 (MQ=255)
                                     cTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCACAGGGAACACGtttt                                                                                                   >  1:204102/1‑139 (MQ=255)
                                        tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAACCCCCCCTCACAACATCCTGCATACGGAACACGTATTccc                                                                                                >  2:605122/1‑139 (MQ=255)
                                            tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                            <  1:250497/107‑1 (MQ=255)
                                            tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                            >  2:250497/1‑107 (MQ=255)
                                            tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCtgt                                                                                            <  2:599102/139‑1 (MQ=255)
                                             ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCtgtg                                                                                           <  2:204102/139‑1 (MQ=255)
                                             ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                     <  2:20497/97‑1 (MQ=255)
                                             ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                     >  1:20497/1‑97 (MQ=255)
                                                  tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                     >  1:795967/1‑92 (MQ=255)
                                                  tttGAAACTCTTCTTCCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                     <  2:795967/92‑1 (MQ=255)
                                                      aaaCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                     <  1:723326/88‑1 (MQ=255)
                                                      aaaCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                     >  2:723326/1‑88 (MQ=255)
                                                                 tctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                      <  1:195202/139‑1 (MQ=255)
                                                                    aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                       >  1:504698/1‑88 (MQ=255)
                                                                    aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                       <  2:504698/88‑1 (MQ=255)
                                                                    aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCACCCCTCACAACGTCCTGCATAGTGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcccccc                                                                   >  1:881695/1‑139 (MQ=255)
                                                                     aaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCATCAAAACGACTTACATACGGAAAACGTTTTCCCATGTGCACCCTCAATTAAATTTccccccc                                                                  >  2:195202/1‑139 (MQ=255)
                                                                         aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa                                                              <  1:312355/139‑1 (MQ=255)
                                                                          aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                                <  1:863062/105‑1 (MQ=255)
                                                                          aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                                >  2:863062/1‑105 (MQ=255)
                                                                            cAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTaaa                                                            <  2:213571/139‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAACCCCACCTCAAACCGTCCTGCATAGGGAACAAGATTTACCCTGTGCACCAACGACTAAATTTCCACACCTATAAACTTCTCaaaaacat                                              >  1:997307/1‑136 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                            >  2:1084734/1‑62 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                            <  1:1084734/62‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                          <  2:410317/80‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                          >  1:410317/1‑80 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCtgtg                                                                                           >  1:196647/1‑95 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCtgtg                                                                                           <  2:196647/95‑1 (MQ=255)
                                                                                                       atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                <  2:997307/139‑1 (MQ=255)
                                                                                                       atatTAACTTGTCCAATATGATCAAATAGCATTAACCGCCCCTTAACAACGGCATACATAGGCAACACGTTTTCTCAACTACAACCACGACTAAATTTACCCCGCTTTAAACTTCCCAATAATCACACCaactaactca                                >  2:1005315/1‑131 (MQ=255)
                                                                                                          ttAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                             <  1:605122/139‑1 (MQ=255)
                                                                                                            aaCTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATc                           <  2:881695/139‑1 (MQ=255)
                                                                                                                 gTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                            <  2:468944/38‑1 (MQ=255)
                                                                                                                 gTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                            >  1:468944/1‑38 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                                <  1:863099/65‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                                >  2:863099/1‑65 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACAAATTTTCTGGAAATCCTatga                     >  1:1063822/1‑139 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCAACCCCTAACAACGTTCTGCATAGGGAAAAAGTTTTACCCCGTGTACCCAAGACTACATTTCCCCCCCTTTAAAATTTTCAATAATATCACCCATTTTCTGCAAATCCTatga                     >  1:515907/1‑139 (MQ=255)
                                                                                                                       tatGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCt                <  2:515907/139‑1 (MQ=255)
                                                                                                                            tCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCGTAGGGGACACGCTTTCCCCTGTGCACGCACGACTAAATTTCCCCCCCTTGAAGCTTCTCAATAATGACCGCAATTTTCTGGAAACCCCATAATGCCTTTccc           >  2:530144/1‑139 (MQ=255)
                                                                                                                                   gCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt    <  2:1063822/139‑1 (MQ=255)
                                                                                                                                     aTTAACCCCCCCCTCACAACGTCCTGCATATGGAACACGTTTTCCCCTGTGCACTCACGAATAAATTTCCCCCCCTTTGAACTTTTCACTAATGTCACCAATTTTCTGGAATTTCTATGATGCCTTTCCAGGTAAATgt  >  2:28235/1‑138 (MQ=255)
                                                                                                                                     aTTAACCCACCCCTCACAACGTCCAGCATAGGGAAAACGCTTTCCCCCGTGCACCCACGACTAAATATCCCCCCCTTTAAATTTCCCACTATTGTCACCAATTTACTGTACATCCTATGATGCCTGTCCCGGTAAGTgg  >  1:182484/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                         
AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGG  >  NC_000913/2767301‑2767571

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