Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 702,066 G→C 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,0660GC97.8% 127.0 / NA 45P229R (CCT→CGT) nagAN‑acetylglucosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (25/19);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTT  >  NC_000913/701927‑702204
                                                                                                                                            |                                                                                                                                          
ggTGGCCGCAGTAACCTGACAAAAAATGACGACTTTCAGACGTTTAGCGTTGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGtt                                                                                                                                         >  1:639869/1‑143 (MQ=255)
        cggtaattatctccaGTTATTCTCCTATAAGAGGTTTAGCGTGGCGAATGTT‑GGCGTAAGCAACATGCAGCCCATCAGCAATACCAGCGCAATAAATGTCAGCTTCGTCGAGGATCGAGCCCGCCAGGCCACGTTCACGACCg                                                                                                                                 <  1:405497/130‑1 (MQ=255)
        cAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCg                                                                                                                                 >  1:205869/1‑143 (MQ=255)
              ccAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCCGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTGCACGACCGGTAata                                                                                                                           >  2:19200/1‑143 (MQ=255)
               cctccacaGTTTGTCGCCTTTCAGCCGTGTACCTTTGCGAATGTT‑GCCGTAATCAACATGCAGGCCATTAGCAATGATACCGCAATACATGTCAGCTTCGTAGAGGATCGCGCCCGCCAGCCCACGTTCACGACCGGTAatat                                                                                                                          <  1:686591/139‑1 (MQ=255)
                       gttcgctgCCTCTCAGCCGTTGAGCGTCG‑GTACGTTGGGCGGAATCAACATGCCTGCCATCAGCAATAATACCGCAATCATTGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCTACGTTCACGACCGGTAATATACGGCATc                                                                                                                  <  1:194359/136‑1 (MQ=255)
                          tGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATcgcg                                                                                                               >  2:525034/1‑143 (MQ=255)
                               ccTTTCAGACGTTTAGCGTTGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCAGCGCGTTGTa                                                                                                          >  2:478313/1‑143 (MQ=255)
                                  ttCAGACGTTTAGCGTTGCGAAGGTG‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAg                                                                                                       <  2:34369/143‑1 (MQ=255)
                                    cAGACGTTTAGCGTTGCGAATGTT‑TGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCTAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGAt                                                                                                     <  2:329256/143‑1 (MQ=255)
                                        cGTTTAGCGTTGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGt                                                                                                 <  2:771491/143‑1 (MQ=255)
                                              gCGTGGCGAGTGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCaaa                                                                                           <  2:27556/143‑1 (MQ=255)
                                               cGTTGCAAATGTT‑GACTTAATCAACATGCAGGTCATACGCGATAATACCGCAAGAAATGGGAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATCCGGCATCGCGTTGTACAGATGGGTGGCAAAg                                                                                          <  2:288653/143‑1 (MQ=255)
                                                 ttGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGt                                                                                        <  2:754762/143‑1 (MQ=255)
                                                  tccGAATGTT‑GGCGTAATCAACATGCAGGCCATCATCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTa                                                                                       <  2:237842/141‑1 (MQ=255)
                                                         tgtggtcGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGGCGAGGATCGCGCGCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATcccc                                                                                 <  1:174401/137‑1 (MQ=255)
                                                             ggCGTCATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCcgcgc                                                                             <  2:298796/143‑1 (MQ=255)
                                                                 tAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATCTAAGACAGCACGTAGTACCGATGGGTGGGAAAGGTGATCCCCGCGCGGaa                                                                         >  1:522158/1‑143 (MQ=255)
                                                                   aTCAACATGCAGGCCATCAGCAATAATACCGCAATAGATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAGTCCCCGCGCGGAAAc                                                                       >  2:571011/1‑143 (MQ=255)
                                                                      aaCATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGGTCACGACCGGTAATATACGGCATCGCGAGCTAGAGATGCGGGGCCAAGCTCCTCCCCCCGCGGCAAccag                                                                    >  2:812774/1‑141 (MQ=255)
                                                                      aaCATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCAGCCAGGCCACGTTCACGACCGGTAATATAAGGCATCGCGTTGTACAGGTGGGTGGAAAAGGTAATCAGCGCAGGGaaagtgg                                                                    >  2:320228/1‑139 (MQ=255)
                                                                       aCATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATGGCGCCCGCCAGGCCACGCTCACGACCGGTAATATAAGGCAGCGCCTGGGCCAGATGGGGGGCCAGCGCAATCCGCGCGGCGGAACCGGc                                                                   >  1:877856/1‑143 (MQ=255)
                                                                             aGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCAAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGGACATGTGGGTGGCAAAGGTACCCCCCGCGCGGAAACCGGCTTTTGc                                                             >  2:862494/1‑143 (MQ=255)
                                                                                  atcatcAATAATACCTATATAATTTTTTGTTTCGTTGAGGTTCTCGCCCGCCAGGCCACGTTCACGACCGGTAACATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGcttctt                                                        <  2:79800/143‑1 (MQ=255)
                                                                                   tCAGCACTAATGCCGTGAGGAATGGCAGCTGCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCttt                                                       <  1:805557/143‑1 (MQ=255)
                                                                                   tCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCttt                                                       >  2:882420/1‑143 (MQ=255)
                                                                                         ataataCCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGt                                                 >  2:113520/1‑143 (MQ=255)
                                                                                              aCCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGAAATCCCCGGGCGGGAACCGGCTTTTGCTTCTTACAACGACGCGt                                            >  2:479489/1‑143 (MQ=255)
                                                                                                cGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGCTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCGCTTGCTTCTTTCGACGGCGCGtgg                                          >  1:282330/1‑141 (MQ=255)
                                                                                                 gCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATCTACGGCATCGCGTTGTACAGATGGGCGGCAAAGGTAAT‑CGCGCGCGGAAACCGCCCTTTGCTCCTCTCacagtccccgtggc                                        >  1:40131/1‑130 (MQ=255)
                                                                                                  gaagaagggTCAGCGTCGTCGAGGCTCGGGCCCGGCGGGCCAGGTTCACGGCCGGTAATATACGGCATGGCGTTGTACAGATGGGTGGCAAAGGTAATGCCGGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGa                                        <  1:166847/135‑1 (MQ=255)
                                                                                                  ctttaatggTCAGCTTCGTCGACGTTCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAAGCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGa                                        <  1:213013/135‑1 (MQ=255)
                                                                                                  cgaTACATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGa                                        <  1:876895/141‑1 (MQ=255)
                                                                                                  cAATAACTGTCAGCTCCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGGCCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCCTCGTCGCGTTGGa                                        <  1:537460/143‑1 (MQ=255)
                                                                                                  cAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGCTGGGGGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGCTGGa                                        >  2:719351/1‑143 (MQ=255)
                                                                                                   aaTAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGCGGCAACTGGAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTGGGAg                                       >  1:101512/1‑143 (MQ=255)
                                                                                                            cAGCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGTCCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGATAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCa                              <  2:555588/143‑1 (MQ=255)
                                                                                                              gCTTCGTCGAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGGTGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGGCCGGCAGa                            >  2:70711/1‑143 (MQ=255)
                                                                                                                   gtcgAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACca                       <  2:369940/143‑1 (MQ=255)
                                                                                                                    tcgAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACcac                      >  1:32235/1‑143 (MQ=255)
                                                                                                                      gAGGATCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGGAACCACaa                    >  2:63562/1‑143 (MQ=255)
                                                                                                                           tCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCACCGTCGCGTTGGAGTGACCGGCTGacgcccacatcccg               >  2:32164/1‑130 (MQ=255)
                                                                                                                           tCGCGCCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGATGGCAAAGGTCATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCg               >  2:354010/1‑143 (MQ=255)
                                                                                                                             gcgcCCGCCAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCgtg             >  2:276697/1‑141 (MQ=255)
                                                                                                                                 ccgccAGGCCACGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCAttt         >  2:255366/1‑143 (MQ=255)
                                                                                                                                     cagggcacgTTCACGGCCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCACAGGCAATCCCCGCGCGGCAACCCGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCa     <  1:785338/143‑1 (MQ=255)
                                                                                                                                        gccacgTTCACGGCCGGTAATATACCGCATCGCGTTAAACAGTTTTGGGTCAAAGGTAAGCCCCGCGCGGAACCCGGCCTTTGCTTCTTTCCACGGCGCGTTGGAGTGACCGGCAGAAACCACAAGCCCGGCATTTGCCAGtt  <  1:619499/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
GGTGGCGTCAGTAACCAGACACAGTTTGTCGCCTTTCAGACGTTTAGCGTTGCGAATGTT‑GGCGTAATCAACATGCAGGCCATCAGCAATAATACCGCAATAAATGTCAGCTTCGTCGAGGATCGCGCCCGCCAGGCCAGGTTCACGACCGGTAATATACGGCATCGCGTTGTACAGATGGGTGGCAAAGGTAATCCCCGCGCGGAAACCGGCTTTTGCTTCTTTCAACGTCGCGTTGGAGTGACCGGCAGAAACCACAATCCCGGCATTTGCCAGTT  >  NC_000913/701927‑702204

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