Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 702,066 G→C 54.8% P229R (CCT→CGT)  nagA ← N‑acetylglucosamine‑6‑phosphate deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913702,0660GC54.8% 0.1 / 53.8 31P229R (CCT→CGT) nagAN‑acetylglucosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base C (8/9);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

CAATGTTCGACCAGATTACGCACGCCTTCAATCATGGTTAAGGATGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGT  >  NC_000913/701781‑702328
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                      
cAATGTTCGACCAGATTACGCACGCCTTCAATCATGGTTAAGGATGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGt                                                                                                                                                                                                                                                                      >  2:251794/1‑288 (MQ=255)
      tCGACCAGATTACGCACCCCTTCAATCATGGTTAAGGCTGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGAc                                                                                                                                                                                                                                                               <  2:149598/289‑1 (MQ=255)
        gACCAGATTACGCACGCCTTCAATCATGGTTAAGGATGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCg                                                                                                                                                                                                                                                             >  2:114027/1‑289 (MQ=255)
         caagcattccacacgcCTTCAATCAGGGTTAAGGATGACCCGCTTACCGTACCTTTACCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCgg                                                                                                                                                                                                                                                            <  2:52435/278‑1 (MQ=255)
         aCCAGATTACGCACGCCTTCAATCATGGTTAAGGATGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCgg                                                                                                                                                                                                                                                            <  2:19124/289‑1 (MQ=255)
                        ccTTCAATCATGGTTAAGGATGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATc                                                                                                                                                                                                                                              >  2:179922/1‑288 (MQ=255)
                        ccTTCAATCATGGTTAAGGATGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGACCGCCAGGCCAGGTTCACGACCGGTAATATAAGGccttt                                                                                                                                                                                                                                             >  1:52435/1‑285 (MQ=255)
                          tacAATCATGGTTAATGATGAACCACTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATcgcg                                                                                                                                                                                                                                           <  2:88699/287‑1 (MQ=255)
                                   ggTTAAGGATGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGa                                                                                                                                                                                                                                  >  2:107813/1‑289 (MQ=255)
                                             gAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCaa                                                                                                                                                                                                                        <  1:251794/289‑1 (MQ=255)
                                              aaCCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGTCGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCAAGGCCACTGTCACGACCGCTTATATCCGGCATCACGTGGTAGAGAGGGGTGGCgaa                                                                                                                                                                                                                       >  2:66088/1‑289 (MQ=255)
                                               aCCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATcgcg                                                                                                                                                                                                                                           <  1:239995/268‑1 (MQ=255)
                                               aCCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATcgcg                                                                                                                                                                                                                                           >  2:239995/1‑268 (MQ=255)
                                                   cTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGtt                                                                                                                                                                                                                                         >  1:182259/1‑266 (MQ=255)
                                                   cTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGtt                                                                                                                                                                                                                                         <  2:182259/266‑1 (MQ=255)
                                                                                                    ttttACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATCCGGCATCGCGTTGTACAGATGGGTGGCAACGGTAATCCCCGAGCCGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTg                                                                                                                                                                 >  2:53494/1‑289 (MQ=255)
                                                                                                                          ttCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCa                                                                                                                                                                                                                                                >  2:96423/1‑188 (MQ=255)
                                                                                                                          ttCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCa                                                                                                                                                                                                                                                <  1:96423/188‑1 (MQ=255)
                                                                                                                               tGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCa                                                                                                                                                          >  2:167782/1‑269 (MQ=255)
                                                                                                                               tGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCa                                                                                                                                                          <  1:167782/269‑1 (MQ=255)
                                                                                                                                    gCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATcgcg                                                                                                                                                                                                                                           <  2:190388/183‑1 (MQ=255)
                                                                                                                                    gCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATcgcg                                                                                                                                                                                                                                           >  1:190388/1‑183 (MQ=255)
                                                                                                                                                              aaccagtaccaGTTTTTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGa                                                                                                          <  1:254421/277‑1 (MQ=255)
                                                                                                                                                                                          cGTTGAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGt                                                                             <  1:114027/287‑1 (MQ=255)
                                                                                                                                                                                              tAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAg                                                                                                            >  2:12953/1‑252 (MQ=255)
                                                                                                                                                                                              tAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAg                                                                                                            <  1:12953/252‑1 (MQ=255)
                                                                                                                                                                                                       gAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGt                                                               <  1:179922/288‑1 (MQ=255)
                                                                                                                                                                                                                                 ccATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAgcg                                       <  1:107813/286‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaaTGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTCGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAAtt              <  1:53494/289‑1 (MQ=255)
                                                                                                                                                                                                                                                              aGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAACTCGACCAGCGAGGCATCAGGCTTACGcacacactttggtata       >  2:278740/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                   cgtcgAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGt  >  1:103930/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                      
CAATGTTCGACCAGATTACGCACGCCTTCAATCATGGTTAAGGATGAACCGCTTAACGTACCGTTCTCATCCACACAAAGTCCGTTACGGTAGTATATTGTTTTACCCGCAAAAATGAACTGTTCAATGTTGGCACCTGCTGGCGCGGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTTGGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTTTGCTGATGACTTCCGCAGGAACCATTTCCGGTGCCAGGGTCACTTTGGTAATGACGTCGGCGTTTTCACACAGGAAATCGACCAGCGCGGCATCAGGCTTACGCACAAAATTCGGATTATGGGT  >  NC_000913/701781‑702328

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