Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,845 G→A 7.3% G46R (GGG→AGG)  yfjP → 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,8450GA 7.3% 186.1 / 3.8 55G46R (GGG→AGG) yfjPCP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base A (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAACACTTTATGACCATCGTGGATCTGCC  >  NC_000913/2767713‑2767984
                                                                                                                                    |                                                                                                                                           
taaaaTTTTGAAGTTCTGCAGCCACTTCAGAACTCCCTTTAGGTTCTACAGCTATGGGTATCTGAACGTATATTGCAGCTGATAAATCAGTTAACCCACTACGAGCCAGTGAACGGCATCATGGGTAAAACTGGGGCggg                                                                                                                                      <  2:926804/139‑1 (MQ=255)
aaaaaTTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTCCCGCGATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCggg                                                                                                                                      <  1:1004085/140‑1 (MQ=255)
aaaaaTTATGAAGTTCTGCACCCACTACAGAACTCCCTTTCCGGTATACCGCTAGGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCgtg                                                                                                                                      <  2:926834/140‑3 (MQ=39)
                 gCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCa                                                                                                                     <  1:1068576/140‑1 (MQ=255)
                       aCTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGccc                                                                                                               >  2:763496/1‑140 (MQ=255)
                            aGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGttt                                                                                                          >  2:757824/1‑140 (MQ=255)
                              aaCTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTAGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGc                                                                                                        >  1:374765/1‑140 (MQ=255)
                              aaCTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTAGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGc                                                                                                        >  1:374495/1‑140 (MQ=255)
                               aCTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGcc                                                                                                       >  1:736829/1‑140 (MQ=255)
                                   ccTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTg                                                                                                   >  1:283750/1‑140 (MQ=255)
                                   ccTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTg                                                                                                   >  2:431495/1‑140 (MQ=255)
                                                ccGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGt                                                                                      <  1:1155227/140‑1 (MQ=255)
                                                   cTATGGGTATCTGAACGTCTCTTGCAGCAGATAAATCAGTTACCCCACTCCGTGCCTGTGATCGGCTTCATGGGTAAAACTGGGGCGGGAATGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTTAAGTATCGCCGGTCAg                                                                                   <  2:253938/140‑1 (MQ=37)
                                                   cTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAg                                                                                   <  2:831654/140‑1 (MQ=255)
                                                      tGGGTATCTGAACGTATATTGCAGCCGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGa                                                                                >  2:1130383/1‑140 (MQ=255)
                                                        ggTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATg                                                                              >  2:737276/1‑140 (MQ=255)
                                                        ggTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATg                                                                              >  2:319269/1‑140 (MQ=255)
                                                        ggTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATg                                                                              >  2:319271/1‑140 (MQ=255)
                                                          tATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTa                                                                            <  1:178396/140‑1 (MQ=255)
                                                          tATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTa                                                                            <  1:852619/140‑1 (MQ=255)
                                                          tATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTa                                                                            <  2:996823/140‑1 (MQ=255)
                                                           aTCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGAGGTAg                                                                           >  2:524236/1‑140 (MQ=255)
                                                             cTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCg                                                                         <  2:654498/140‑1 (MQ=255)
                                                               gAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGc                                                                       <  2:259686/140‑1 (MQ=255)
                                                               gAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGc                                                                       <  2:709047/140‑1 (MQ=255)
                                                               gAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGc                                                                       <  1:229217/140‑1 (MQ=255)
                                                               gAACGTATATTGCAGCAGAGAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGc                                                                       <  1:990164/140‑1 (MQ=255)
                                                                                gATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCtt                                                      >  1:475323/1‑140 (MQ=255)
                                                                                gATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCtt                                                      >  1:475271/1‑140 (MQ=255)
                                                                                gATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCtt                                                      >  1:475239/1‑140 (MQ=255)
                                                                                gATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCtt                                                      >  1:475234/1‑140 (MQ=255)
                                                                                                cACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCa                                      <  1:1019150/140‑1 (MQ=255)
                                                                                                cACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCa                                      <  1:851615/140‑1 (MQ=255)
                                                                                                cACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCa                                      <  1:278118/140‑1 (MQ=255)
                                                                                                  cTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGa                                    <  1:781864/140‑1 (MQ=255)
                                                                                                   tACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGAt                                   >  2:94826/1‑140 (MQ=255)
                                                                                                   tACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGAt                                   >  2:214930/1‑140 (MQ=255)
                                                                                                   tACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGAt                                   >  2:214953/1‑140 (MQ=255)
                                                                                                       aGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGga                               >  1:1025782/1‑140 (MQ=255)
                                                                                                       aGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGga                               >  1:762592/1‑140 (MQ=255)
                                                                                                        gCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGgag                              >  1:1184166/1‑140 (MQ=255)
                                                                                                          cAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTGTCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGaa                            <  1:731043/140‑1 (MQ=255)
                                                                                                          cAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGaa                            <  1:775279/140‑1 (MQ=255)
                                                                                                          cAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGaa                            <  1:171705/140‑1 (MQ=255)
                                                                                                          cAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGaa                            <  1:171719/140‑1 (MQ=255)
                                                                                                          cAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGaa                            <  1:427110/140‑1 (MQ=255)
                                                                                                          cAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGaa                            <  1:775100/140‑1 (MQ=255)
                                                                                                           aGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAac                           <  1:214953/140‑1 (MQ=255)
                                                                                                           aGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAac                           <  1:214930/140‑1 (MQ=255)
                                                                                                             tGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAacac                         >  2:522047/1‑140 (MQ=255)
                                                                                                             tGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAacac                         >  1:971378/1‑140 (MQ=255)
                                                                                                                  ggCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAACACTTTAt                    >  2:146877/1‑140 (MQ=255)
                                                                                                                  ggCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAACACTTTAt                    >  2:146801/1‑140 (MQ=255)
                                                                                                                    catcatGGGTAAAACTAGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAACACTTTATGa                  <  2:374765/140‑1 (MQ=255)
                                                                                                                    catcatGGGTAAAACTAGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAACACTTTATGa                  <  2:374495/140‑1 (MQ=255)
                                                                                                                                    ggggCGGGTAAGAGTAGGCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGAAGTTTGGAGAACACTTTATGACCATCGTGGATCTGcc  <  1:689133/140‑1 (MQ=255)
                                                                                                                                    ggggCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAACACTTTATGACCATCGTGGATCTGcc  <  1:698680/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                           
AAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGAGCCAGTGATCGGCATCATGGGTAAAACTGGGGCGGGAAAGAGTAGCCTTTGCAATGCCCTGTTTGCCGGTGAAGTATCGCCGGTCAGCGATGTAGCGGCCTGTACACGTGATCCTTTACGCTTTCGCCTGCAGATTGGAGAACACTTTATGACCATCGTGGATCTGCC  >  NC_000913/2767713‑2767984

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