Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,882 T→G 27.7% I163S (ATC→AGC)  yagG → CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913285,8820TG27.7% 81.1 / 41.7 47I163S (ATC→AGC) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base T (17/17);  new base G (6/7);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

GACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGT  >  NC_000913/285745‑286023
                                                                                                                                         |                                                                                                                                                
gACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCAGCGCGATGCCg                                                                                                                                        <  4:897245/148‑1 (MQ=255)
  cccTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATcg                                                                                                                                                >  1:227359/1‑138 (MQ=255)
  cccTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATcg                                                                                                                                                <  2:227359/138‑1 (MQ=255)
                 tCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCgcgc                                                                                                                                             <  4:1075341/126‑1 (MQ=255)
                 tCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCgcgc                                                                                                                                             >  3:1075341/1‑126 (MQ=255)
                                                   catcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTggg                                                                                               <  1:354337/135‑1 (MQ=255)
                                                   catcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTggg                                                                                               >  2:354337/1‑135 (MQ=255)
                                                       aCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgt                                                                              >  3:926430/1‑148 (MQ=255)
                                                                   aaaGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGCGCATCATCGGCaaa                                                                                                                >  2:170249/1‑102 (MQ=255)
                                                                   aaaGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCaaa                                                                                                                <  1:170249/102‑1 (MQ=255)
                                                                            cACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑Gt                                                                                                                                <  4:188702/77‑1 (MQ=255)
                                                                            cACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑Gt                                                                                                                                >  3:188702/1‑77 (MQ=255)
                                                                                    gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTg                                                                                                                               >  1:1029415/1‑70 (MQ=255)
                                                                                    gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTg                                                                                                                               <  2:1029415/70‑1 (MQ=255)
                                                                                     cAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctgc                                                <  4:926430/148‑1 (MQ=255)
                                                                                      aGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCaa                                                                                                                 >  1:124626/1‑82 (MQ=255)
                                                                                      aGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCaa                                                                                                                 <  2:124626/82‑1 (MQ=255)
                                                                                            ggCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTg                                                                                                 <  4:646381/92‑1 (MQ=255)
                                                                                            ggCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTg                                                                                                 >  3:646381/1‑92 (MQ=255)
                                                                                                 ttcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                              >  4:344661/1‑90 (MQ=255)
                                                                                                 ttcttcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                              <  3:344661/90‑1 (MQ=255)
                                                                                                         ggcggcggcGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGcatca                                                                                                                        <  3:235788/56‑1 (MQ=255)
                                                                                                         ggcggcggcGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGcatca                                                                                                                        >  4:235788/1‑56 (MQ=255)
                                                                                                               ggcggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgtg                                                                             >  4:682003/1‑93 (MQ=255)
                                                                                                               ggcggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCAtgtg                                                                             <  3:682003/93‑1 (MQ=255)
                                                                                                                           cgctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTg                                                                >  1:332689/1‑94 (MQ=255)
                                                                                                                           cgctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTg                                                                <  2:332689/94‑1 (MQ=255)
                                                                                                                            gctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCt         <  4:575220/148‑1 (MQ=255)
                                                                                                                            gctATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCt         >  3:575220/1‑148 (MQ=255)
                                                                                                                               accAGCGGCAGCGCGATGCCGCTG‑‑‑ATGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcg                                                                           >  4:1025917/3‑79 (MQ=37)
                                                                                                                               accAGCGGCAGCGCGATGCCGCTG‑‑‑ATGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcg                                                                           <  3:1025917/77‑1 (MQ=38)
                                                                                                                                ccAGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTc                                        <  4:193549/112‑1 (MQ=255)
                                                                                                                                ccAGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTc                                        >  3:193549/2‑113 (MQ=255)
                                                                                                                                  aGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtggtg                                                    <  2:744373/99‑1 (MQ=255)
                                                                                                                                  aGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtggtg                                                    >  1:744373/1‑99 (MQ=255)
                                                                                                                                  aGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAgg   >  3:171459/1‑148 (MQ=255)
                                                                                                                                   gCGGCAGCGCGATGCCGCTGATAGCGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCt                                                                 <  3:633498/88‑1 (MQ=25)
                                                                                                                                   gCGGCAGCGCGATGCCGCTGATAGCGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCt                                                                 >  4:633498/1‑88 (MQ=25)
                                                                                                                                   gCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGt  <  4:171459/148‑1 (MQ=255)
                                                                                                                                    cGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCt               <  4:154627/134‑1 (MQ=255)
                                                                                                                                    cGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCt               >  3:154627/1‑134 (MQ=255)
                                                                                                                                    cGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCaa                       >  3:412329/1‑126 (MQ=255)
                                                                                                                                    cGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCaa                       <  4:412329/126‑1 (MQ=255)
                                                                                                                                       cATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCg                                                                                       >  3:273177/1‑59 (MQ=255)
                                                                                                                                       cATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCg                                                                                       <  4:273177/59‑1 (MQ=255)
                                                                                                                                        aTCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACAcc          >  1:789384/1‑135 (MQ=255)
                                                                                                                                        aTCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACCACCAAAGAGCGCTAcccc          <  2:789384/135‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
GACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGT  >  NC_000913/285745‑286023

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