Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 705,916 T→C 12.7% intergenic (+26/‑177) nagE → / → glnS N‑acetylglucosamine specific PTS enzyme IIABC component/glutamine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913705,9160TC12.7% 94.1 / 5.3 44intergenic (+26/‑177)nagE/glnSN‑acetylglucosamine specific PTS enzyme IIABC component/glutamine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base C (4/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

CGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGC  >  NC_000913/705774‑706058
                                                                                                                                              |                                                                                                                                              
cGCTCGATGATTAGCCCGGTAGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCtt                                                                                                                                                <  1:51938/143‑1 (MQ=255)
  atCGATGATGAGGCCGGTGGTGTGCAGAAAAATCGAAGAAAGCAGCGGCTTGAACATAAAAGCACAAAGAAAAATTGTGGCGGGTCAAACATCGAAAAAAAAAAAAAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGa                                                                                                                                              <  2:348866/142‑1 (MQ=255)
     gatgatTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCAACGCTTGAGCg                                                                                                                                           >  1:503714/1‑143 (MQ=255)
            aGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTCAAGCGTCGCATc                                                                                                                                    <  2:340606/143‑1 (MQ=255)
                  gtggtTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCcaacaa                                                                                                                              <  1:782813/143‑1 (MQ=255)
                             caATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGt                                                                                                                   <  1:530600/143‑1 (MQ=255)
                                 atCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGgag                                                                                                               <  1:674628/143‑1 (MQ=255)
                                 atCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGgag                                                                                                               <  1:323677/143‑1 (MQ=255)
                                  tCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGCGCGTCGCGTCCAACACTGCCAAGCGGCGggga                                                                                                              >  1:353660/1‑143 (MQ=255)
                                  tCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCCACAAGGACAAGCGGTGgaga                                                                                                              >  1:135067/1‑143 (MQ=255)
                                       gaTTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGACAAAAACCGCTGTAAGAAATAAAAAAGTAATCTGCATTATGCCTGATGCGACACATGAACGACGCAAACAACAATAACCAGCGGTAGAAAacttc                                                                                                         >  2:235320/1‑143 (MQ=255)
                                         tttCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTtctc                                                                                                       >  1:486096/1‑143 (MQ=255)
                                          ttCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTctct                                                                                                      >  1:759514/1‑143 (MQ=255)
                                                   ttGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCtttt                                                                                             >  1:489065/1‑143 (MQ=255)
                                                      aTCATTAAAGCTCAGCGCCATATTGTGGCGGGGCAAACACCGCTGTATGAAATCAAAAAGTACTCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGctttcttt                                                                                          <  1:795930/143‑1 (MQ=255)
                                                           tAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTtcat                                                                                     >  1:545783/1‑143 (MQ=255)
                                                               gCTCAGGGCCCTATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCatcaat                                                                                 <  1:306218/143‑1 (MQ=255)
                                                                 tCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGAACTTCTCTGGCGCTTTTTTTTTCatcaatga                                                                               >  1:70541/1‑141 (MQ=255)
                                                                  cAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAatcat                                                                              >  1:748666/1‑143 (MQ=255)
                                                                               tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTGATCAATCATACCCACAACCCtt                                                                 >  2:349608/1‑143 (MQ=255)
                                                                               tgGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCtt                                                                 >  1:5900/1‑143 (MQ=255)
                                                                                   gggTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATGTTCTCTGCCGCTCTTTTTTTCATCAATCATCCCCATAATACTTGtca                                                             >  2:726570/1‑141 (MQ=255)
                                                                                   gggTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCATGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTa                                                             <  1:616126/143‑1 (MQ=255)
                                                                                    ggTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTGGTTAg                                                            >  1:55/1‑143 (MQ=255)
                                                                                     gTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGCGATCTTCTCTGCCGCTTTTTTTTCCATCGATCATCCCCATAATCCTTGTTCGa                                                           >  1:353736/1‑143 (MQ=255)
                                                                                        aaaCACGGCTGTGTGAAATCAAAAAGTGATCTGTTTTAGGGCTGAGGCGACGCTCGACTGTCGCATCCAACAATGTCAAGCTGTGTAGATTTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTa                                                        <  2:266640/143‑1 (MQ=255)
                                                                                        aaaCACCGCTGTATGAAATCAAAAAGTAATCTGCGGTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTa                                                        <  1:437649/143‑1 (MQ=255)
                                                                                             ccGTTGTGTGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAt                                                   <  2:415049/143‑1 (MQ=255)
                                                                                               gccgtgtGAAATAAAAAAGTGATCTGCGTTAGGCCTGATGCGACGCGTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAttt                                                 <  1:217185/137‑1 (MQ=255)
                                                                                                cTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGGTTATCAAtttt                                                >  1:519902/1‑143 (MQ=255)
                                                                                                         aTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAGTAAAcga                                       >  2:705935/1‑141 (MQ=255)
                                                                                                          tCAAAAAGTAATCTGCTTTATGGCGGATGCGACGCGTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTaa                                      <  2:759514/143‑1 (MQ=255)
                                                                                                          tCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGCTTATCAATTTTACAAAACTaa                                      >  1:247335/1‑143 (MQ=255)
                                                                                                               aaGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGtt                                 <  2:5900/143‑1 (MQ=255)
                                                                                                                aGTACTCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCACCCAACAATGACAAGCGGTTGATATCTTCTCTGCCGCTTTTTTTTTCATCCATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGGTg                                <  2:790094/143‑1 (MQ=255)
                                                                                                                  tAATGTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTc                              <  2:387448/143‑1 (MQ=255)
                                                                                                                            attgttggATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGCCCTTCTCTGCCGCTTTTTTTTTCCTCAATAATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTTTCAGCCTGTccc                    >  1:40373/8‑143 (MQ=255)
                                                                                                                            attgttggATGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCCTGTTATATTGTCAATTTCAACTTCCTATCAGATGTCAGCATGTgcc                    >  1:733999/8‑143 (MQ=255)
                                                                                                                              tGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTGGATTATCAATTTTAAAAAACTAACAGATGTCAGCCTGTccccc                  >  1:635523/1‑141 (MQ=255)
                                                                                                                                 cTGGTGCGACGCTTGGGCGTCGCATCCAACGATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTta               <  1:548864/143‑1 (MQ=255)
                                                                                                                                   gATGCGACGCTTGAGCGTCGCATACAACTATGACATGCGGGGGCGTTCTTCGTGGCGGCTTTTTTTATCATCAATCCGCCCCATAATCCTTGTTAGATTAGCAATATTAAAAAACTAACAGGTGTCAGCCTGTCCCGCTtata             <  2:249126/143‑1 (MQ=255)
                                                                                                                                    aTGCGACGCTCAAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCTTCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTACAAAACTTCCAGTTGTCAGCCTGTCCCGCTTTTaa            >  2:109888/1‑143 (MQ=255)
                                                                                                                                    aTGCGACGCTCAAGCGTAGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATGATCCACCTAATACTTGTTAGATTATCAATTTTAAAACACTAACAGTTGTCAGGCTGTCCCGCTGATaa            >  1:545479/1‑143 (MQ=255)
                                                                                                                                     tGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCCTCCCCATAATCCGTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAg           <  2:698935/143‑1 (MQ=255)
                                                                                                                                          cGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCAt      <  2:543836/143‑1 (MQ=255)
                                                                                                                                              tGAGCGGCGCATCCAACAATGACCAGCGGTGGAGATCTTCTCCGCCGCTTTTTTTTTCATCAATCATCCCCATAAGCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGc  <  2:217196/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CGCTCGATGATTAGCCCGGTGGTTTGCAGCAATATCGACGATTTCAGTGGCTTGATCATTAAAGCTCAGGGCCATATTGTGGCGGGTCAAACACCGCTGTATGAAATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGC  >  NC_000913/705774‑706058

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