Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,665 T→G 9.0% intergenic (+312/‑43) 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,6650TG 9.0% 91.4 / 2.7 42intergenic (+312/‑43)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (22/15);  new base G (2/3);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

ACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACC  >  NC_000913/2767527‑2767806
                                                                                                                                          |                                                                                                                                             
aCCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCCGATAGCTATCATTCTCCTGAACAATCCGCTGAACGGCAGAACCTCGCGGTTCTg                                                                                                                                     >  2:68703/1‑149 (MQ=255)
aCCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                                                                     >  1:235188/1‑149 (MQ=255)
    aTTTTCTGGAAATCCTATTATGCCTTTCCCGGTAAATGGTGCGGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                                                                 <  1:380669/149‑1 (MQ=255)
    aTTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                                                                 <  1:380688/149‑1 (MQ=255)
           ggAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATCCTCCTGGACAATCCGATGACCGGCAGAACGTCGCGGTTCTGCCGTTTTTGtt                                                                                                                          >  1:320229/1‑149 (MQ=255)
            gAAATCCTATGATGCCTTTCCCGGTAACTGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGCCAGAACCGCGAGGGTCTGTCATCTTCGttt                                                                                                                         >  2:782236/1‑149 (MQ=255)
                  cTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGACAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAAttt                                                                                                                   >  1:329769/1‑149 (MQ=255)
                    atgatgCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                                                                 <  1:385450/149‑1 (MQ=255)
                        tggcTTTCCCGGTAAATGGTGCTGTGAGTTCAGGAAGATTCTGCAATTTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCagag                                                                                                             <  2:10087/146‑1 (MQ=255)
                        tgCCTTCGCCCGTAATTCGGGCTGTGAGTTTAGATAGTTTCTTCAATTTTTGTGCATTCAATAAATTTTCAGATCTCTATCCTTCTTCTGTACAATCTGTTTCACGGTAGAACATCGCGGTTCTGCCGTTTTTGTTTTCCtttctcgag                                                                                                             <  2:15026/149‑7 (MQ=255)
                             ttCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtg                                                                                                        >  2:3166/1‑149 (MQ=255)
                                   gTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgatg                                                                                                  >  1:611392/1‑149 (MQ=255)
                                          gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATAATCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGCGGTGATGAAAAAtt                                                                                           >  1:529324/1‑149 (MQ=255)
                                            gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                                                         <  1:353232/149‑1 (MQ=255)
                                                    ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAAACCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGGGATGAAAAATTTTGACGTTCt                                                                                 >  2:24134/1‑149 (MQ=255)
                                                     tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTg                                                                                <  1:497689/149‑1 (MQ=255)
                                                     tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTg                                                                                >  1:500254/1‑149 (MQ=255)
                                                     tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCCTGAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTg                                                                                <  1:346472/149‑1 (MQ=255)
                                                         aaaGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCCCTCTGCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGc                                                                            >  1:157913/1‑149 (MQ=255)
                                                                cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCAcgcca                                                                     >  1:73560/1‑145 (MQ=255)
                                                                              cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg                                                       >  2:493698/1‑149 (MQ=255)
                                                                              cATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTGCTGCAGCCACTACAGAACTCCCTTTCCg                                                       <  1:480955/149‑1 (MQ=255)
                                                                                   aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTa                                                  >  1:707341/1‑149 (MQ=255)
                                                                                        aTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGACCTCCCTTTCCGGTCTACCGCt                                             >  2:701682/1‑149 (MQ=255)
                                                                                        aTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCt                                             >  2:743880/1‑149 (MQ=255)
                                                                                        aTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAAACTCGCGGCTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCTCTACATAACTCCCTTTCCGGTCTACCGcg                                             >  2:503183/1‑148 (MQ=255)
                                                                                           ttCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATg                                          <  2:235188/149‑1 (MQ=255)
                                                                                            acAGATAGATATCATTCTCCTGAACAACCCGATGAACCGCAGAACCTCGCGCTTCTGCCGTTTTTCTTTTCCCTTCCAGAGGAGTGGTGATGAAAAACTTTGAAGTTCTGCAGCCACTACAGAACGCCCTTTCCGGTCTACCGCTATgg                                         <  2:334443/148‑1 (MQ=255)
                                                                                               gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTGTGGGtc                                      >  1:27856/1‑148 (MQ=255)
                                                                                                      atCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAAc                               >  1:687589/1‑149 (MQ=255)
                                                                                                      atCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAAc                               <  1:3166/149‑1 (MQ=255)
                                                                                                          ttCTCCGGAACAATCCTATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGtat                           <  1:118258/149‑1 (MQ=255)
                                                                                                              ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTg                       <  1:101324/149‑1 (MQ=255)
                                                                                                              ccTGAACAATCCGATGAACGGCAGAACCGCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTg                       <  2:191472/149‑1 (MQ=255)
                                                                                                               cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTgc                      >  1:276425/1‑149 (MQ=255)
                                                                                                                  aaCAATCCGATGAACGGCAGAACCTTGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTgcagc                   <  1:367823/149‑1 (MQ=255)
                                                                                                                       tCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGGGGTGAAAAAATTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTCTGGGGAGCGGAACGTATATTGCAGCAGATa              >  1:131235/1‑149 (MQ=255)
                                                                                                                              aaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAGAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCAATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGt       >  2:475708/1‑149 (MQ=255)
                                                                                                                                 ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTTTTTTCAGAGGAGGGGTGATGAAAAATTTTGAAGTTATGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTAGGGGTATCTGAACGTCGATTGCGGCAGATAAATCAGTTaa    >  2:633457/1‑149 (MQ=255)
                                                                                                                                 ggCAGAACCTCGCCGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGCAGTTCTGCAGCCCCTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTGaa    >  2:492819/1‑149 (MQ=255)
                                                                                                                                  gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAAc   <  1:493698/149‑1 (MQ=255)
                                                                                                                                   cAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAAcc  <  2:647142/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                             
ACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACC  >  NC_000913/2767527‑2767806

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑