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.C85.7% 56.5 / 3.6 21intergenic (+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 (6/12);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 6.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAT  >  NC_000913/2767301‑2767569
                                                                                                                                           |                                                                                                                                     
aGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAc                                                                                                                                      >  2:89809/1‑139 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                   >  1:159942/1‑139 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccc                                                                                                                                >  1:879641/1‑139 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                 >  2:851262/1‑138 (MQ=255)
    tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                 <  1:851262/138‑1 (MQ=255)
           tgcctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                         <  2:189122/136‑1 (MQ=255)
             tctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccc                                                                                                                               <  2:965466/131‑1 (MQ=255)
             tctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccc                                                                                                                               >  1:965466/1‑131 (MQ=255)
             tctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTGGAACCTCTCTTTCTAAAAGAACCAGTTTAGAATTCCCCCTTTTTGCATGATACTTAATTGTCCAAACAGATCAAATAGCATTAAccccccccccccaaa                                                                                                                      >  1:189122/1‑131 (MQ=255)
                 ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                    <  2:545053/139‑1 (MQ=255)
                 ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                    <  1:89809/139‑1 (MQ=255)
                              ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                       <  2:159942/139‑1 (MQ=255)
                                 aGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGCTATTAACTTGTCCAATATGATCAAATAGCATTAAAGCCCTCCTCACAACGCACTGCATAGGGAACACg                                                                                                    >  1:1038072/1‑139 (MQ=255)
                                      tGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                         <  1:881506/113‑1 (MQ=255)
                                      tGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                         >  2:881506/1‑113 (MQ=255)
                                        tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                             <  2:879641/139‑1 (MQ=255)
                                                      aaaCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGa                                                                               <  2:927157/139‑1 (MQ=255)
                                                                            cAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAAc                                                         <  2:919135/139‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGACATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCa                                                                                      <  1:85670/97‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGACATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCa                                                                                      >  2:85670/1‑97 (MQ=255)
                                                                                           ccTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                             >  2:1030627/1‑88 (MQ=255)
                                                                                           ccTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                             <  1:1030627/88‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACCACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                    >  2:615419/1‑90 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                    <  1:615419/90‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACACCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTTATCAATTTTCTGGAAATCCCatga                   >  2:1093042/1‑139 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCTCCTCTACCAACGGCGTGAATAAGGCACACGGTTTACCCACTCCAACGAAGACTACGTTGCCCCCCCTTTCAAGCTGTCAGTCATCGCAGGAATCTTCTGGTAATCCGataa                   >  2:1134491/1‑137 (MQ=255)
                                                                                                                                   gCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTATAAACTTCTCAATAATGTAACCAATTTTCTGGAAATCCTATGATACCTTTCCCGGTAAAt  >  1:754171/1‑139 (MQ=255)
                                                                                                                                     aTTAACCCCCCCCTCACAACGTCCTGGATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                    >  2:94179/1‑71 (MQ=255)
                                                                                                                                     aTTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                    <  1:94179/71‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                     
AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAT  >  NC_000913/2767301‑2767569

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: