Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 264,717 (A)5→6 coding (16/477 nt) ykfG ← CP4‑6 prophage; RadC‑like JAB domain‑containing protein YkfG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913264,7121.A97.4% 133.6 / ‑3.4 38coding (21/477 nt)ykfGCP4‑6 prophage; RadC‑like JAB domain‑containing protein YkfG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (16/21);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CCGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTGGATGT  >  NC_000913/264632‑264787
                                                                                 |                                                                           
ccGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt                                                                              <  1:61383/81‑1 (MQ=255)
ccGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt                                                                              <  1:61382/81‑1 (MQ=255)
ccGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt                                                                              <  1:61385/81‑1 (MQ=255)
ccGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt                                                                              <  1:61384/81‑1 (MQ=255)
 cGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt                                                                              >  2:61427/1‑80 (MQ=255)
 cGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt                                                                              >  2:61436/1‑80 (MQ=255)
 cGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTa                                                                            <  1:61386/81‑1 (MQ=255)
  gggCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTaa                                                                           <  2:61387/81‑1 (MQ=255)
   ggCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTaaa                                                                          <  1:61388/81‑1 (MQ=255)
      tCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTaaaa                                                                         >  2:61447/1‑79 (MQ=255)
      tCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTaaaa                                                                         >  2:61440/1‑79 (MQ=255)
         tGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTaaaaaagg                                                                     >  2:61429/1‑77 (MQ=255)
          gCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTaaaaaagga                                                                    >  2:61437/1‑76 (MQ=255)
            aGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAg                                                                  >  1:61443/1‑81 (MQ=255)
               tGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTg                                                               <  2:61390/81‑1 (MQ=255)
                 gCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGtt                                                             >  1:61442/1‑81 (MQ=255)
                 gCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGt                                                              <  1:61391/80‑1 (MQ=255)
                 gCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAGGGAAGCTGttt                                                            >  1:61458/1‑81 (MQ=255)
                   ggTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGttt                                                            >  1:61456/1‑80 (MQ=255)
                   ggTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGttt                                                            >  2:61453/1‑80 (MQ=255)
                   ggTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTc                                                           <  1:61392/81‑1 (MQ=255)
                    gTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCa                                                          >  2:61454/1‑81 (MQ=255)
                      ccAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATc                                                        >  2:61444/1‑81 (MQ=255)
                       cAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATc                                                        <  2:61393/80‑1 (MQ=255)
                                 aGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCtt                                             <  2:61394/81‑1 (MQ=255)
                                    gcCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCgg                                          >  1:61432/1‑81 (MQ=255)
                                      ccGCTGAATGAGACGCCGGTCCTGTGGCGTCGTCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTc                                        <  2:61395/81‑1 (MQ=255)
                                      ccGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTc                                        <  1:61396/81‑1 (MQ=255)
                                       cGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCa                                       <  1:61397/81‑1 (MQ=255)
                                        gCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAg                                      <  1:61398/81‑1 (MQ=255)
                                         cTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGt                                     >  2:61457/1‑81 (MQ=255)
                                             aTGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGat                                 <  1:61399/81‑1 (MQ=255)
                                              tGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGata                                <  2:61400/81‑1 (MQ=255)
                                              tGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGata                                <  1:61401/81‑1 (MQ=255)
                                              tGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGata                                <  1:61402/81‑1 (MQ=255)
                                                 gaCGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGataata                             >  2:61450/1‑81 (MQ=255)
                                                   cGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACg                           <  2:61403/81‑1 (MQ=255)
                                                    gCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGc                          <  2:61405/81‑1 (MQ=255)
                                                    gCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATAAGc                          <  1:61404/81‑1 (MQ=255)
                                                        gTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGaa                      >  1:61477/1‑81 (MQ=255)
                                                           cTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATgg                   >  1:61472/1‑81 (MQ=255)
                                                             gtgGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCg                 >  2:61448/1‑81 (MQ=255)
                                                                gCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGt                >  1:61465/1‑79 (MQ=255)
                                                                gCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTg               <  2:61407/80‑1 (MQ=255)
                                                                gCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGa              <  1:61406/81‑1 (MQ=255)
                                                                 cGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAg             >  1:61451/1‑81 (MQ=255)
                                                                  gTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGc            >  1:61455/1‑81 (MQ=255)
                                                                     aTCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCctc         <  1:61408/81‑1 (MQ=255)
                                                                          gCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTg       <  2:61409/78‑1 (MQ=255)
                                                                          gCCGGGTAAAAAAGGAAAGCTGTTTCATCTATTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTGGa     <  2:61410/80‑1 (MQ=255)
                                                                            cGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTGGATGt  >  1:61471/1‑81 (MQ=255)
                                                                                 |                                                                           
CCGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTGGATGT  >  NC_000913/264632‑264787

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