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

AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGA  >  NC_000913/2767301‑2767552
                                                                                                                                           |                                                                                                                     
aGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTaaa                                                                                                                       >  1:1422389/1‑138 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑acc                                                                                                                  >  2:743919/1‑139 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑ccccc                                                                                                               <  1:109423/139‑1 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑ccccc                                                                                                               >  1:1000212/1‑139 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGGCAATGGATTTATGAAACTCTTCTTTCGAAACGAACCAGTATAGAATTCGCACTTTTGGCATGATATTAACCTGTCCAATATGATCAAATAGCATTAA‑‑ccccc                                                                                                               >  1:976492/1‑139 (MQ=255)
       gTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccccccc                                                                                                            <  2:322691/139‑2 (MQ=255)
            gTCTCCCGATGATGTACGCCAAGTCCTGTCGATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                        <  1:613161/139‑1 (MQ=255)
            gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccccccctcgca                                                                                                       <  2:819473/139‑7 (MQ=255)
            gTCTCCCGATGATGTACGCCAAGTCCAGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACCTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                        <  2:976492/139‑1 (MQ=255)
                 ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                   <  2:1000212/139‑1 (MQ=255)
                              ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                      <  2:1516529/139‑1 (MQ=255)
                                        tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTTACAACGTCCTGCATAGGGAACACGTTTTccc                                                                            <  1:1407764/139‑1 (MQ=255)
                                            tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑c                                                                                                                     >  1:261953/1‑95 (MQ=255)
                                            tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑c                                                                                                                     <  2:261953/95‑1 (MQ=255)
                                                  tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc                                                                                                                 >  2:651362/1‑92 (MQ=255)
                                                  tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc                                                                                                                 <  1:651362/92‑1 (MQ=255)
                                                        aCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc                                                                                                                 >  1:1251569/1‑86 (MQ=255)
                                                        aCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑cccc                                                                                                                 <  2:1251569/86‑1 (MQ=255)
                                                          tcttctTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTaa                                                          <  2:924612/139‑1 (MQ=255)
                                                                         aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa                                           <  2:1422389/139‑1 (MQ=255)
                                                                                         tgtataagagacagATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGt                           <  1:925705/126‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                      <  2:616195/81‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCCTCACAACGTCCTCCATAGGGaaa                                                                                      >  1:616195/1‑80 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCAc                                                                <  2:641780/102‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCAc                                                                >  1:641780/1‑102 (MQ=255)
                                                                                                                tGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTGTCAATAATGTCACCAATTTTCTGGAAATCCTat    >  1:524146/1‑139 (MQ=255)
                                                                                                                 gTCCAATATGATCAAATAGCATTAAA‑‑CCCCGCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAACTTCACCAATATTCTGGAAATCCTatg   >  1:1279786/1‑139 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACTTCCTGCATAGGTAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga  >  1:625888/1‑139 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATa                                                                                            <  2:1087858/49‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATa                                                                                            >  1:1087858/1‑49 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa                                         <  1:464737/100‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTGTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAACCTTCTCAATAATGTCACCAATTTTCTGGAAATcccctgt  >  2:194899/1‑134 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAA‑‑CACCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa                                         >  2:464737/1‑100 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAA‑‑CCCCTCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTTCCCCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCCGGAAATCCTatga  >  2:370148/1‑139 (MQ=255)
                                                                                                                         tttgatcaTAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGAc                                                             >  2:46497/8‑73 (MQ=255)
                                                                                                                         tttgatcaTAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGAc                                                             <  1:46497/66‑1 (MQ=255)
                                                                                                                                    cATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                   <  2:811556/56‑1 (MQ=255)
                                                                                                                                    cATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                   >  1:811556/1‑56 (MQ=255)
                                                                                                                                           |                                                                                                                     
AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAA‑‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGA  >  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: