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.C88.0% 70.3 / 3.5 25intergenic (+83/‑272)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base C (8/14);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 7.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGA  >  NC_000913/2767301‑2767552
                                                                                                                                           |                                                                                                                     
aGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑aaa                                                                                                                       >  1:1422389/1‑138 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑a‑cc                                                                                                                  >  2:743919/1‑139 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑c                                                                                                               <  1:109423/139‑1 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑c                                                                                                               >  1:1000212/1‑139 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGGCAATGGATTTATGAAACTCTTCTTTCGAAACGAACCAGTATAGAATTCGCACTTTTGGCATGATATTAACCTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑c                                                                                                               >  1:976492/1‑139 (MQ=255)
       gTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑cccc                                                                                                            <  2:322691/139‑2 (MQ=255)
            gTCTCCCGATGATGTACGCCAAGTCCTGTCGATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTcaca                                                                                                        <  1:613161/139‑1 (MQ=255)
            gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc‑cccctcgca                                                                                                       <  2:819473/139‑7 (MQ=255)
            gTCTCCCGATGATGTACGCCAAGTCCAGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACCTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTcaca                                                                                                        <  2:976492/139‑1 (MQ=255)
                 ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTc                                                                                                   <  2:1000212/139‑1 (MQ=255)
                              ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAac                                                                                      <  2:1516529/139‑1 (MQ=255)
                                        tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTTACAACGTCCTGCATAGGGAACACGTTTTccc                                                                            <  1:1407764/139‑1 (MQ=255)
                                            tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c                                                                                                                     >  1:261953/1‑95 (MQ=255)
                                            tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c                                                                                                                     <  2:261953/95‑1 (MQ=255)
                                                  tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc                                                                                                                 >  2:651362/1‑92 (MQ=255)
                                                  tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc                                                                                                                 <  1:651362/92‑1 (MQ=255)
                                                        aCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc                                                                                                                 >  1:1251569/1‑86 (MQ=255)
                                                        aCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑c‑ccc                                                                                                                 <  2:1251569/86‑1 (MQ=255)
                                                          tcttctTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTaa                                                          <  2:924612/139‑1 (MQ=255)
                                                                         aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa                                           <  2:1422389/139‑1 (MQ=255)
                                                                                         tgtataagagacagATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGt                           <  1:925705/126‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACCC‑CCC‑CCCTCACAACGTCCTGCATAGGGAac                                                                                      <  2:616195/81‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACCC‑CCC‑CCCTCACAACGTCCTCCATAGGGaaa                                                                                      >  1:616195/1‑80 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCAc                                                                <  2:641780/102‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCAc                                                                >  1:641780/1‑102 (MQ=255)
                                                                                                                tGTCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTGTCAATAATGTCACCAATTTTCTGGAAATCCTat    >  1:524146/1‑139 (MQ=255)
                                                                                                                 gTCCAATATGATCAAATAGCATTAAA‑‑C‑CCC‑GCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAACTTCACCAATATTCTGGAAATCCTatg   >  1:1279786/1‑139 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACTTCCTGCATAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga  >  1:625888/1‑139 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATa                                                                                            <  2:1087858/49‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATa                                                                                            >  1:1087858/1‑49 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa                                         <  1:464737/100‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTGTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAACCTTCTCAATAATGTCACCAATTTTCTGGAAATcccctgt  >  2:194899/1‑134 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AA‑‑CACCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa                                         >  2:464737/1‑100 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AA‑‑C‑CCCTCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTTCCCCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCCGGAAATCCTatga  >  2:370148/1‑139 (MQ=255)
                                                                                                                         tttgatcaTAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGAc                                                             >  2:46497/8‑73 (MQ=255)
                                                                                                                         tttgatcaTAGCATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGAc                                                             <  1:46497/66‑1 (MQ=255)
                                                                                                                                    cATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                   <  2:811556/56‑1 (MQ=255)
                                                                                                                                    cATT‑AAC‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                   >  1:811556/1‑56 (MQ=255)
                                                                                                                                           |                                                                                                                     
AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑‑C‑CCC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGA  >  NC_000913/2767301‑2767552

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: ‑