Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 704,083 C→A A47E (GCG→GAG)  nagE → N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913704,0830CA95.9% 137.1 / ‑5.1 50A47E (GCG→GAG) nagEN‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (24/23);  minor base T (1/1);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAACCCAGAA  >  NC_000913/703948‑704222
                                                                                                                                       |                                                                                                                                           
atatatTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGAttttt                                                                                                                                      <  1:575108/143‑1 (MQ=255)
atagTTTAGGTTTTTTCCAGCGACTCGGTCGGGCGTTACAGCGCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGAttttt                                                                                                                                      <  2:564284/143‑1 (MQ=255)
 tatTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTg                                                                                                                                     <  1:726654/143‑1 (MQ=255)
   ttttAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGAt                                                                                                                                   <  1:676754/143‑1 (MQ=255)
      tAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAAc                                                                                                                                <  1:354720/143‑1 (MQ=255)
         gTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGAGGCACTGTGGCTGCGATTCGGTCAGCCCGATTTCATTAGCGTTGCGTTTATTGCCCACGCGGGAGGTGAGAATTTTAGTAAcatg                                                                                                                             >  2:155810/1‑140 (MQ=255)
            ttgtCCGGCGAGTCGCTACGGCGTACAGCGCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCAt                                                                                                                         <  2:803108/140‑1 (MQ=255)
                ccAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTaa                                                                                                                      <  2:645405/143‑1 (MQ=255)
                               ggCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACATCGCATTAATCTTCGCAATCGgtg                                                                                                       >  2:549907/1‑143 (MQ=255)
                               ggCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCCGTGAGATTTTTGATAACCTCGCACTAATCTTCGCCATCGctg                                                                                                       >  2:677861/1‑143 (MQ=255)
                                gCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGgtgt                                                                                                      <  1:113275/143‑1 (MQ=255)
                                   ttACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTTTGATTTTTGATAATTTAGCATTCATCGTCGCAATCGGTGTGGc                                                                                                   >  1:489307/1‑143 (MQ=255)
                                    gactggtcacTATCGCTGTGGTGCGTGGGGCGGCACTGTTGCTGCGATGCGGTCAGCCAGAGGTACTGAACGTGGCGTCTATTGCCCAGTCGGATTATTGTATTTTGTATAACCTCGCATTAATCTTCGCCATCGGTGTGGCa                                                                                                  <  2:912749/134‑1 (MQ=255)
                                         cTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCCGCCAg                                                                                             >  1:574241/1‑143 (MQ=255)
                                          tCCCTCTCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGGCTATTGCCCAGGCGGGCGGTGAGATTTTTGCTAACCTCGCATTAATCTTCGCCATCGGTGCGGCATCCacc                                                                                            >  2:775951/1‑141 (MQ=255)
                                           cccTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCt                                                                                           <  1:59481/143‑1 (MQ=255)
                                           cccTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTGGCGTTTATTGCCCAGGCGGGCGGTGAGACGTTTGATAACCGCGCATCAGTCCTCGCCCTCGGGGTGGCCTCCCGcc                                                                                           >  2:874913/1‑142 (MQ=255)
                                            ccTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTg                                                                                          >  1:180320/1‑143 (MQ=255)
                                            ccTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTg                                                                                          >  1:818550/1‑143 (MQ=255)
                                                   cgGTGCTGGCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGaaa                                                                                   <  2:194559/143‑1 (MQ=255)
                                                   cgGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGaaa                                                                                   <  1:687683/143‑1 (MQ=255)
                                                   cgGGGCTGCGGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGAGTTACTTAACGTCGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGaaa                                                                                   <  2:772204/143‑1 (MQ=255)
                                                       gctgcCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGAGTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCCCGCCTCAATCGTCGCCGTCGGAGTGGCATCCAGCTgcgcgcaagcga                                                                               >  2:345205/1‑132 (MQ=255)
                                                        ctgcCTGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAg                                                                              <  1:553667/143‑1 (MQ=255)
                                                         tgcCGGGGGCGGCACTGTTGCTGCGAGTGGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAgc                                                                             <  1:283730/143‑1 (MQ=255)
                                                                 gcggcACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTgc                                                                     >  1:101836/1‑143 (MQ=255)
                                                                   ggcACTGTTGCTGCGGTTCGCTCAGCCAGATTTACTTGACTTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTTATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTgcgg                                                                   <  2:674715/143‑1 (MQ=255)
                                                                   ggcACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAATACAGCGCTGGTgcgc                                                                   >  2:118389/1‑142 (MQ=255)
                                                                   ggcACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGAGTTTTGATAACCTCGCAGTAATCTTCGCCATCGGTGTGGCATCCAGATGGTCGAAAGACAGCGCAGGTgcgg                                                                   >  2:873925/1‑143 (MQ=255)
                                                                    gcACGGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTgcggc                                                                  <  2:167443/143‑1 (MQ=255)
                                                                       cctgtgctgcGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCCGTGAGATTTTTTATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTgcggcggc                                                               <  1:329605/139‑1 (MQ=255)
                                                                            gctgcGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGCTGCGGCGGCGCTgg                                                          >  1:681376/1‑143 (MQ=255)
                                                                                      ggccatccAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGt                                                <  1:376867/137‑1 (MQ=255)
                                                                                            ccAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTa                                          >  1:391861/1‑143 (MQ=255)
                                                                                            ccAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCGCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGCTGCGCCGGCGCTGTCGGGTGGCGTAGGTTa                                          >  2:706692/1‑143 (MQ=255)
                                                                                                 tttACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCCTCCAGCTGCTCGAAAGACAGCGCCGGGGCGGCGGCGCTGGCGGGTGccggagttcaccctg                                     >  2:293626/1‑129 (MQ=255)
                                                                                                   tACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTtgtg                                   <  1:405995/143‑1 (MQ=255)
                                                                                                        aaCGTTGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAAc                              <  2:428248/143‑1 (MQ=255)
                                                                                                        aaCGTTGCGTTTATTGCACAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAAc                              <  1:775690/143‑1 (MQ=255)
                                                                                                            ttGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGAGCGGTAGGTCACCCTGTGTTAACCcaa                          >  2:231479/1‑143 (MQ=255)
                                                                                                            ttGCGTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGGGCGGCGGCGCTGGGGTGAGCGGGAGGTTACTTTGTGTTAACCaaa                          >  2:9058/1‑143 (MQ=255)
                                                                                                                gTTTATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGa                      >  1:564761/1‑143 (MQ=255)
                                                                                                                  ttATTGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCTTTAATCTTCGCCATCGGTGTTGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGTGTGCGGTAGGTTGCTCTGTGTTGACGACAGCGacg                    >  2:565959/1‑141 (MQ=255)
                                                                                                                     ttGCCCAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTg                 <  1:565676/143‑1 (MQ=255)
                                                                                                                        cccAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGTCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCCGGCGGGTCCGGTACGATACCTTGCGTTAACCACAGCGaggtcgacc              >  1:643529/1‑135 (MQ=255)
                                                                                                                         ccAGGCGGGCGGTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCa             >  2:593157/1‑143 (MQ=255)
                                                                                                                             gcgggcggTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGGTACTTTGTGTTCACCAAAGCGATGGTGACCATCaa         >  1:67432/1‑143 (MQ=255)
                                                                                                                               gggcggGGAGATTCTTGATCCCCTCGCAGTAATCTTCGCCATCGGTGTGGCATCCCGCCGGCCGAAAGACAGCGCTGGTGCGGCGGCGCTTGCGGGTGCGTTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAAcc       <  1:345205/143‑1 (MQ=255)
                                                                                                                                  gggTGTGAGTGTTGCTAGCGTCGCATGAAGGTGCGCCGTCGGTGTGGCATCCAGCTGGTCGAAAGACAGGGCTGGTGCGGCGGCGGTGGCGGGTTCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAAcccgg    <  2:170593/142‑3 (MQ=255)
                                                                                                                                    gTGAGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAGCCAAAGCGATGGGGACCAGCAACCCAGaa  >  2:461648/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                           
ATATTTTAGGTTTTTTCCAGCGACTCGGTAGGGCGTTACAGCTCCCTATCGCGGTGCTGCCGGTGGCGGCACTGTTGCTGCGATTCGGTCAGCCAGATTTACTTAACGTTGCGTTTATTGCCCAGGCGGGCGGTGCGATTTTTGATAACCTCGCATTAATCTTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAACCCAGAA  >  NC_000913/703948‑704222

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