Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,887 C→A 28.9% L165M (CTG→ATG)  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,8870CA28.9% 82.8 / 37.6 45L165M (CTG→ATG) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base A (6/7);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

GACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGG  >  NC_000913/285745‑286032
                                                                                                                                              |                                                                                                                                                    
gACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCAGCGCGATGCCg                                                                                                                                                 <  4:897245/148‑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                                                             >  1:744373/1‑99 (MQ=255)
                                                                                                                                  aGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtggtg                                                             <  2:744373/99‑1 (MQ=255)
                                                                                                                                  aGCGGCAGCGCGATGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAgg            >  3:171459/1‑148 (MQ=255)
                                                                                                                                   gCGGCAGCGCGATGCCGCTGATAGCGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCt                                                                          >  4:633498/1‑88 (MQ=25)
                                                                                                                                   gCGGCAGCGCGATGCCGCTGATAGCGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCt                                                                          <  3:633498/88‑1 (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                                                                                                <  4:273177/59‑1 (MQ=255)
                                                                                                                                       cATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCg                                                                                                >  3:273177/1‑59 (MQ=255)
                                                                                                                                        aTCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACAcc                   >  1:789384/1‑135 (MQ=255)
                                                                                                                                        aTCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACCACCAAAGAGCGCTAcccc                   <  2:789384/135‑1 (MQ=255)
                                                                                                                                          cgcgCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGGGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTg                                               <  4:76931/105‑1 (MQ=255)
                                                                                                                                          cgcgCTGCCGCTG‑‑‑ATGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTg                                               >  3:76931/1‑105 (MQ=255)
                                                                                                                                            cgcTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCgggg  >  1:873290/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTG‑‑‑GTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGG  >  NC_000913/285745‑286032

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: