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.C90.4% 160.1 / 8.4 52intergenic (+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/1);  new base C (20/27);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAA  >  NC_000913/2767305‑2767568
                                                                                                                                       |                                                                                                                                      
tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccc                                                                                                                                  >  2:1690626/1‑139 (MQ=255)
tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccc                                                                                                                                  >  1:1341166/1‑139 (MQ=255)
tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcaacc                                                                                                                                   >  2:1326195/1‑135 (MQ=255)
 ctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAActcaca                                                                                                                                  >  2:839399/1‑134 (MQ=255)
      ttGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccccctca                                                                                                                           <  1:886362/139‑5 (MQ=255)
        gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑caca                                                                                                                        <  2:1341166/139‑1 (MQ=255)
        gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCCTTCTAAACGAACCAGTATAGAATCCGCCCTTTTGGCATGCTATTAACTTGCCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑caca                                                                                                                        <  2:1141065/139‑1 (MQ=255)
             ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTc                                                                                                                   <  2:449631/139‑1 (MQ=255)
              cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCT‑‑CACAACGTCCt                                                                                                                 >  1:660878/1‑139 (MQ=255)
                       aCGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAggg                                                                                                         >  2:1185531/1‑139 (MQ=255)
                       aCGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAggg                                                                                                         <  2:1771511/139‑1 (MQ=255)
                         gCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCT‑‑CACAACGTCCTGCAAAGGGAac                                                                                                      <  2:660878/139‑1 (MQ=255)
                              gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAac                                                                                                      >  1:852370/1‑135 (MQ=255)
                              gTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAac                                                                                                      <  2:852370/135‑1 (MQ=255)
                                ccTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑caca                                                                                                                        <  2:331790/115‑1 (MQ=255)
                                ccTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑caca                                                                                                                        >  1:331790/1‑115 (MQ=255)
                                 cTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGtttt                                                                                               <  1:310649/139‑1 (MQ=255)
                                    tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                            <  1:1690626/139‑1 (MQ=255)
                                    tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                            <  1:1185531/139‑1 (MQ=255)
                                         ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAACCCCCCCT‑‑CACAACGTCCTGAATAGGGAACACGCTTTCCCCtata                                                                                       >  1:358263/1‑136 (MQ=255)
                                         ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGtt                                                                                                 <  1:1556177/129‑1 (MQ=255)
                                         ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGtt                                                                                                 >  2:1556177/1‑129 (MQ=255)
                                         ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCtgtg                                                                                       >  2:1134128/1‑139 (MQ=255)
                                              tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGGCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAccc                                                                                  <  1:1134128/139‑1 (MQ=255)
                                                      tcttctTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACt                                                                            <  1:73396/139‑1 (MQ=255)
                                                             tctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                   <  1:1326195/139‑1 (MQ=255)
                                                                aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTc                                                                                                                   >  1:425546/1‑88 (MQ=255)
                                                                aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTc                                                                                                                   <  2:425546/88‑1 (MQ=255)
                                                                  aCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCTCT‑‑CACAACATCCTACATAGCGACCACATATGACACTGTGCACCCACGACTACATTTcccacaca                                                              >  2:790449/1‑134 (MQ=255)
                                                                     aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa                                                           <  1:839399/139‑1 (MQ=255)
                                                                      aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                            >  2:982277/1‑105 (MQ=255)
                                                                      aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                            <  1:982277/105‑1 (MQ=255)
                                                                                aaTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACAAGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCaat                                                >  1:1149142/1‑139 (MQ=255)
                                                                                aaTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                   >  1:896646/1‑58 (MQ=255)
                                                                                aaTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                   <  2:896646/58‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCtgt                                                                                        <  1:1579542/94‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCtgt                                                                                        >  2:1579542/1‑94 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGa                                                                              >  1:1164757/1‑104 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGa                                                                              <  2:1164757/104‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTTAATAATGt                                           >  1:1521809/1‑139 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCACCCCT‑‑CACAACGTCCCGCATAGGGAACAAGTTTTCGCCTGCGCACCCCCGACTAAATTTTCCCCCCTTAACACTACCCaaaaatct                                           >  1:599/1‑137 (MQ=255)
                                                                                                   atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACCACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                             <  2:599/139‑1 (MQ=255)
                                                                                                             gTCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATATGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGACACCAATTTTCTGGAAATCCTatg                   >  2:114412/1‑139 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCa                                                                                 >  2:326181/1‑76 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCa                                                                                 <  1:326181/76‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa                                                           >  1:807842/1‑98 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTa                                                           <  2:807842/98‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                  >  1:46295/1‑139 (MQ=255)
                                                                                                                   tatGATCAAATAGCATTAACCCCCCCCT‑‑CCCAACGTCCTCCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa                                                         >  2:1392651/1‑95 (MQ=255)
                                                                                                                   tatGATCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTaaa                                                         <  1:1392651/95‑1 (MQ=255)
                                                                                                                       aTCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAac                                                                                                      >  2:604199/1‑46 (MQ=255)
                                                                                                                       aTCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAac                                                                                                      <  1:604199/46‑1 (MQ=255)
                                                                                                                        tCAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTCACCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTTTCGAAATCCTATGATGCCTTTccc        >  2:1671390/1‑139 (MQ=255)
                                                                                                                         cAAATAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg       <  1:790449/139‑1 (MQ=255)
                                                                                                                             tAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGtttt                                                                                               >  2:1037398/1‑47 (MQ=255)
                                                                                                                             tAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGtttt                                                                                               <  1:1037398/47‑1 (MQ=255)
                                                                                                                             tAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCtgtg                                                                                       >  1:163909/1‑55 (MQ=255)
                                                                                                                             tAGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCtgtg                                                                                       <  2:163909/55‑1 (MQ=255)
                                                                                                                              aGCATTAACCCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTTAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa  <  2:1521809/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCT‑‑CACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAA  >  NC_000913/2767305‑2767568

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