Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,771,806 A→T intergenic (+191/+196) yfjT → / ← yfjU CP4‑57 prophage, putative periplasmic protein/CP4‑57 prophage, conserved protein,Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,771,8060AT96.5% 196.0 / ‑0.8 58intergenic (+191/+196)yfjT/yfjUCP4‑57 prophage, putative periplasmic protein/CP4‑57 prophage, conserved protein,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (24/32);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGATCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGATATGAATGAATACTGCGGATATCAAC  >  NC_000913/2771676‑2771923
                                                                                                                                  |                                                                                                                     
ttCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACtttt                                                                                                               <  2:435849/139‑1 (MQ=255)
 tctcCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACttttt                                                                                                              >  1:891223/1‑139 (MQ=255)
     cTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTc                                                                                                          >  2:421719/1‑139 (MQ=255)
       tAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTa                                                                                                        <  1:906098/139‑1 (MQ=255)
           ctcCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCATTTTTTGGTCTAAGCg                                                                                                    <  1:740748/139‑1 (MQ=255)
           ctcCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCg                                                                                                    <  2:935379/139‑1 (MQ=255)
           ctcCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCg                                                                                                    <  2:891223/139‑1 (MQ=255)
                   gTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTc                                                                                            <  2:863977/139‑1 (MQ=255)
                      ggCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACttt                                                                                                                >  2:102784/1‑116 (MQ=255)
                      ggCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACttt                                                                                                                <  1:102784/116‑1 (MQ=255)
                      ggCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGt                                                                                               <  2:1088564/133‑1 (MQ=255)
                      ggCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGt                                                                                               >  1:1088564/1‑133 (MQ=255)
                      ggCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCata                                                                                         >  1:126428/1‑139 (MQ=255)
                        ccACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATg                                                                                       <  2:126428/139‑1 (MQ=255)
                           cGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTc                                                                                            <  1:370050/131‑1 (MQ=255)
                           cGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTc                                                                                            >  2:370050/1‑131 (MQ=255)
                             tACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCaa                                                                                  <  1:421719/139‑1 (MQ=255)
                                         aTGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATa                                                                      <  1:117949/139‑1 (MQ=255)
                                            aGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAAc                                                                   >  2:1078526/1‑139 (MQ=255)
                                            aGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGATCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAAc                                                                   >  1:426468/1‑139 (MQ=255)
                                                   tGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATgag                                                            >  2:425382/1‑139 (MQ=255)
                                                   tGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATgag                                                            >  1:440545/1‑139 (MQ=255)
                                                   tGGTCAGTTTATTACAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATgag                                                            >  1:743689/1‑139 (MQ=255)
                                                    ggTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGt                                                                                                           <  2:280232/91‑1 (MQ=255)
                                                    ggTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGt                                                                                                           >  1:280232/1‑91 (MQ=255)
                                                    ggTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCat                                                                                          <  2:467288/108‑1 (MQ=255)
                                                    ggTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCat                                                                                          >  1:467288/1‑108 (MQ=255)
                                                    ggTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATgaga                                                           <  2:919683/139‑1 (MQ=255)
                                                      tCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTTCTGTTCACTTTTTGGTCTAAGCg                                                                                                    >  1:338431/1‑96 (MQ=255)
                                                      tCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCg                                                                                                    <  2:338431/96‑1 (MQ=255)
                                                        aGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGa                                                       >  2:885259/1‑139 (MQ=255)
                                                        aGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGa                                                       <  2:241161/139‑1 (MQ=255)
                                                                                ttttttGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCAttaaaattaaaa                               <  1:363149/139‑1 (MQ=255)
                                                                                 tttttGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTaaaaa                              >  2:519509/1‑139 (MQ=255)
                                                                                   tttGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGATCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGa                            <  2:426468/139‑1 (MQ=255)
                                                                                     tGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGata                          <  1:885259/139‑1 (MQ=255)
                                                                                     tGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGata                          <  1:620110/139‑1 (MQ=255)
                                                                                           tAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCg                                                                                                    >  1:92732/1‑59 (MQ=255)
                                                                                           tAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCg                                                                                                    <  2:92732/59‑1 (MQ=255)
                                                                                                  cGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGt                                                                                               >  2:517796/1‑57 (MQ=255)
                                                                                                  cGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGt                                                                                               <  1:517796/57‑1 (MQ=255)
                                                                                                   gCTGTGAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGATATGAATGAATACTGc            <  1:1078526/139‑1 (MQ=255)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCgtt                                                                                                  <  1:748548/49‑3 (MQ=255)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCgtt                                                                                                  >  2:748548/1‑47 (MQ=38)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCg                                                                                                    >  1:840845/1‑47 (MQ=255)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCg                                                                                                    <  2:840845/47‑1 (MQ=255)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGata                          >  2:504996/1‑121 (MQ=255)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGata                          <  1:504996/121‑1 (MQ=255)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGATatgaat                     >  2:797886/1‑126 (MQ=255)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGATatgaat                     <  1:797886/126‑1 (MQ=255)
                                                                                                       tgAGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGATATGAATGAATACTGCGGat        <  2:743689/139‑1 (MQ=255)
                                                                                                         aGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATgaga                                                           >  2:898611/1‑86 (MQ=255)
                                                                                                         aGCGCAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATgaga                                                           <  1:898611/86‑1 (MQ=255)
                                                                                                             cAGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGATATGAATGAATACTGCGGATATCAAc  <  1:519509/139‑1 (MQ=255)
                                                                                                              aGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCata                                                                                         <  1:1099709/51‑1 (MQ=255)
                                                                                                              aGGCTGTGTTAAAACTCCTGTTCACTTTTTGGTCTAAGCGATGGTTTCata                                                                                         >  2:1099709/1‑51 (MQ=255)
                                                                                                                              ccTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGata                          >  2:345340/1‑98 (MQ=255)
                                                                                                                              ccTGTTCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGata                          <  1:345340/98‑1 (MQ=255)
                                                                                                                                  |                                                                                                                     
TTCTCCTTAATCTCCAGAGGTAGGCCACGTACTTTCGAACGATGAGCGCAGTGGTCAGTTTTTTTCAGATATATTTATCTTTTTTTGGCTGTAAGCTCCGCTGTGAGCGCAGGCTGTGTTAAAACTCCTGATCACTTTTTGGTCTAAGCGATGGTTTCATATGATCAACCATTCTTGATAAACCCATGAGACGGACTAACTTTTTCATTAAAATTAAAAAGATATGAATGAATACTGCGGATATCAAC  >  NC_000913/2771676‑2771923

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