Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,799,301 A→G 8.7% A330A (GCT→GCC) ‡ ansB ← glutaminase‑asparaginase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,799,3010AG 8.7% 126.1 / 3.2 46A330A (GCT→GCC) ‡ansBglutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (2/2);  total (23/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 = 8.62e-01

GGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATACCGAACC  >  NC_002947/2799154‑2799447
                                                                                                                                                   |                                                                                                                                                  
ggaatcGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCAGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                                                                                                                                                    <  8:77725/143‑1 (MQ=255)
ggCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                                                                                                                                                    <  6:106158/148‑1 (MQ=255)
      gAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGa                                                                                                                                             >  4:23324/1‑149 (MQ=255)
           ccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAg                                                                                                                                        >  4:30686/1‑149 (MQ=255)
           ccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAg                                                                                                                                        >  8:23625/1‑149 (MQ=255)
            cTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGt                                                                                                                                       >  4:330787/1‑149 (MQ=255)
             tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcct                                                                                                                                          <  1:162553/145‑2 (MQ=255)
             tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcct                                                                                                                                          >  2:162553/1‑144 (MQ=255)
                   ccGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCtt                                                                                                                                <  6:243471/149‑1 (MQ=255)
                    cGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTg                                                                                                                               <  7:277237/149‑1 (MQ=255)
                    cGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTg                                                                                                                               <  2:64240/149‑1 (MQ=255)
                          aCCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATcc                                                                                                                         >  4:205151/1‑149 (MQ=255)
                              gCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCt                                                                                                                     <  3:330787/149‑1 (MQ=255)
                                ccTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGt                                                                                                                   <  5:202572/149‑1 (MQ=255)
                                   gAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCg                                                                                                                >  4:81330/1‑149 (MQ=255)
                                                  aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTAc                                                                                                  <  3:23324/148‑1 (MQ=255)
                                                       tCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcctga                                                                                                                                        >  4:308681/1‑102 (MQ=255)
                                                       tCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcctga                                                                                                                                        <  3:308681/105‑4 (MQ=255)
                                                         ccAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAaccac                                                                                          >  2:318077/1‑149 (MQ=255)
                                                          cAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAaccacc                                                                                         >  7:166043/1‑149 (MQ=255)
                                                          cAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCATGTCATCCGGCTGTTCGGCATTACGCAGTACGAAaccacc                                                                                         >  1:157380/1‑149 (MQ=255)
                                                           agaagaTCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCAccc                                                                                        >  5:291052/1‑149 (MQ=255)
                                                            gaagaTCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCt                                                                                       >  5:38488/1‑149 (MQ=255)
                                                             aagaTCGGCTGCAGCTCTTTGGTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTg                                                                                      <  7:70331/149‑1 (MQ=255)
                                                              agaTCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTgg                                                                                     <  3:30686/149‑1 (MQ=255)
                                                               gaTCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCAtt                                                                                                           >  1:12383/1‑126 (MQ=255)
                                                               gaTCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCAtt                                                                                                           <  2:12383/126‑1 (MQ=255)
                                                                            tctTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGcgac                                                                         >  7:48709/1‑147 (MQ=255)
                                                                                       cTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCt                                                            <  8:166043/149‑1 (MQ=255)
                                                                                        tGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTg                                                           >  6:357565/1‑149 (MQ=255)
                                                                                                ggTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGc                                                      <  3:205151/146‑1 (MQ=255)
                                                                                                ggTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGt                                                   >  1:233279/1‑149 (MQ=255)
                                                                                                 gTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGtt                                                  <  2:233279/149‑1 (MQ=255)
                                                                                                       gCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTgc                                             <  6:38488/148‑1 (MQ=255)
                                                                                                           ccATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCa                                          <  1:318077/147‑1 (MQ=255)
                                                                                                             aTCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGc                                                          >  5:356583/1‑129 (MQ=255)
                                                                                                             aTCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGc                                                          <  6:356583/129‑1 (MQ=255)
                                                                                                                gCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTg                                   <  2:157380/149‑1 (MQ=255)
                                                                                                                   aGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTg                                                           <  7:270789/122‑1 (MQ=255)
                                                                                                                   aGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTg                                                           >  8:270789/1‑122 (MQ=255)
                                                                                                                         cgcgCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                          <  8:48709/149‑1 (MQ=255)
                                                                                                                          gcgcTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCa                                          >  8:95658/1‑132 (MQ=255)
                                                                                                                          gcgcTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCa                                          <  7:95658/132‑1 (MQ=255)
                                                                                                                          gcgcTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                          >  1:364084/1‑148 (MQ=255)
                                                                                                                          gcgcTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                          >  5:20509/1‑148 (MQ=255)
                                                                                                                                tttGCGGGTTCAGGTCGTGAGCCACGAACCAGACGTTCTTGTCATACGGCTGTTCTGCATTACGCAGTACGAAACCACCCTGGTTGACATGAGACGATCGGATGATCTGAACGACGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTAcc                   >  2:213182/1‑149 (MQ=255)
                                                                                                                                           aGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATAc        <  2:364084/149‑1 (MQ=255)
                                                                                                                                                 tGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATACCGAAcc  <  3:81330/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATACCGAACC  >  NC_002947/2799154‑2799447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑