Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,419,159 G→A *163* (TAG→TAA
S51S (AGC→AGT
sieB →
ydaF ←
phage superinfection exclusion protein, Rac prophage
uncharacterized protein, Rac prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,419,1590GA95.6% 128.2 / ‑1.6 45*163* (TAG→TAA
S51S (AGC→AGT
sieB
ydaF
phage superinfection exclusion protein, Rac prophage
uncharacterized protein, Rac prophage
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (24/19);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-01

AAAAATTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCG  >  NC_000913/1419026‑1419294
                                                                                                                                     |                                                                                                                                       
aaaaaTTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTAc                                                                                                                                    >  2:670331/1‑139 (MQ=255)
 aaaaTTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACt                                                                                                                                   <  2:907616/139‑1 (MQ=255)
     ttAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgt                                                                                                                               <  2:465376/139‑1 (MQ=255)
         aaCGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCtt                                                                                                                           >  2:355568/1‑139 (MQ=255)
                   catcatTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCAt                                                                                                                 >  2:575272/1‑139 (MQ=255)
                                cAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAg                                                                                                    >  2:1023620/1‑139 (MQ=255)
                                cAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAg                                                                                                    >  1:64826/1‑139 (MQ=255)
                                     tGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCg                                                                                               <  1:1023620/139‑1 (MQ=255)
                                       aTAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                                              >  2:153860/1‑106 (MQ=255)
                                       aTAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                                              <  1:153860/106‑1 (MQ=255)
                                             cAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTaaa                                                                                       >  2:119620/1‑139 (MQ=255)
                                             cAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTaaa                                                                                       >  2:10206/1‑139 (MQ=255)
                                               aaCTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATg                                                                                     >  1:636413/1‑139 (MQ=255)
                                                     tCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                                              >  2:937920/1‑92 (MQ=255)
                                                     tCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                                              <  1:937920/92‑1 (MQ=255)
                                                       ttattaTCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCtt                                                                             <  2:636413/139‑1 (MQ=255)
                                                            atCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCtt                                                                             <  1:1241722/134‑1 (MQ=255)
                                                            atCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCtt                                                                             >  2:1241722/1‑134 (MQ=255)
                                                            atCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGt                                                                        >  2:399702/1‑139 (MQ=255)
                                                             tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                                       <  1:119620/139‑1 (MQ=255)
                                                             tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                                       <  1:364559/139‑1 (MQ=255)
                                                              ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGtat                                                                      >  1:1014578/1‑139 (MQ=255)
                                                                cGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATAc                                                                    >  1:474416/1‑139 (MQ=255)
                                                                  tCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGt                                                                  >  2:1109057/1‑139 (MQ=255)
                                                                  tCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGt                                                                  >  1:444821/1‑139 (MQ=255)
                                                                              aCTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGaa                                                      <  1:670331/139‑1 (MQ=255)
                                                                                      tttATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCAt                                              <  2:64826/139‑1 (MQ=255)
                                                                                        tATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATcc                                            <  2:1014578/139‑1 (MQ=255)
                                                                                         aTGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCt                                           <  1:10206/139‑1 (MQ=255)
                                                                                                tATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                            >  1:3253/1‑99 (MQ=255)
                                                                                                tATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                            <  2:3253/99‑1 (MQ=255)
                                                                                                  ttCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGa                                  >  1:652531/1‑139 (MQ=255)
                                                                                                   tCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGaa                                 <  2:444821/139‑1 (MQ=255)
                                                                                                          cagcagGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg                          <  1:1109057/139‑1 (MQ=255)
                                                                                                          cagcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                                              <  1:5093/87‑1 (MQ=255)
                                                                                                          cagcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                                              >  2:5093/1‑87 (MQ=255)
                                                                                                             cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                            >  1:1097237/1‑86 (MQ=255)
                                                                                                             cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                            <  2:1097237/86‑1 (MQ=255)
                                                                                                                     gACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGa                                                       <  1:118963/99‑1 (MQ=255)
                                                                                                                     gACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGa                                                       >  2:118963/1‑99 (MQ=255)
                                                                                                                     gACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGa               >  1:99688/1‑139 (MQ=255)
                                                                                                                        gTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTc            >  2:817325/1‑139 (MQ=255)
                                                                                                                          tttAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCg                                                                                               >  2:181091/1‑54 (MQ=255)
                                                                                                                          tttAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCg                                                                                               <  1:181091/54‑1 (MQ=255)
                                                                                                                                  ttAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCg  >  2:185476/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                       
AAAAATTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCG  >  NC_000913/1419026‑1419294

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