Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,169,170 T→C 8.9% M263T (ATG→ACG)  PP_4547 → glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,169,1700TC 8.9% 123.8 / 3.1 45M263T (ATG→ACG) PP_4547glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base C (2/2);  total (23/22)
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.83e-01

GTGGAAATCACCCTGGAGCATGGCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCGCTGCTGCTGTTCTGCC  >  NC_002947/5169022‑5169307
                                                                                                                                                    |                                                                                                                                          
gTGGAAATCACCCTGGAGCATGGCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGcat                                                                                                                                            <  2:42890/149‑1 (MQ=255)
       tCACCCTGGAGCATGGCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGcatgcacatg                                                                                                                                     >  4:301704/1‑149 (MQ=255)
         aCCCTGGAGCATGGCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGcatgcacatgca                                                                                                                                   >  2:56064/1‑149 (MQ=255)
          cccTGGAGCATGGCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCatgcacatgcac                                                                                                                                  <  4:350888/149‑1 (MQ=255)
             tGGAGCATGGCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCt                                                                                                                                >  5:108393/1‑148 (MQ=255)
               gAGCATGGCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAAGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGa                                                                                                                              <  8:158184/148‑1 (MQ=255)
                    tGGCGATGTGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTgg                                                                                                                        <  5:347211/149‑1 (MQ=255)
                     ggCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGc                                                                                                                       >  6:137786/1‑149 (MQ=255)
                     ggCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCACGCACATGCACCTGAGCCTGGc                                                                                                                       >  1:10560/1‑149 (MQ=255)
                     ggCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCACGCACATGCACCTGAGCCTGGc                                                                                                                       >  1:4605/1‑149 (MQ=255)
                                        tGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCaacaac                                                                                                   <  3:301704/149‑1 (MQ=255)
                                        tGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCaacaac                                                                                                   <  1:264262/149‑1 (MQ=255)
                                          gACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAAcc                                                                                                  <  6:108393/148‑1 (MQ=255)
                                             cAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGcaac                                                                                                      >  2:144199/1‑141 (MQ=255)
                                             cAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGcaac                                                                                                      <  1:144199/141‑1 (MQ=255)
                                             cAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGtt                                                                                              >  7:85324/1‑149 (MQ=255)
                                             cAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCACGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGtt                                                                                              <  2:4605/149‑1 (MQ=255)
                                             cAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCACGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGtt                                                                                              <  2:10560/149‑1 (MQ=255)
                                              aGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTc                                                                                             <  7:34793/149‑1 (MQ=255)
                                              aGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTc                                                                                             <  4:174321/149‑1 (MQ=255)
                                                        cTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCg                                                                                      >  3:261907/1‑146 (MQ=255)
                                                         tACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACa                                                                                                                  <  4:213312/118‑1 (MQ=255)
                                                         tACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACa                                                                                                                  >  3:213312/1‑118 (MQ=255)
                                                            aaGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGa                                                                                  >  6:7060/1‑146 (MQ=255)
                                                            aaGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGa                                                                                  <  5:7060/146‑1 (MQ=255)
                                                             aGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAg                                                                              >  3:136902/1‑149 (MQ=255)
                                                             aGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAg                                                                              >  4:286852/1‑149 (MQ=255)
                                                              gcgcTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACACGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAg                                                                              <  1:56064/149‑1 (MQ=255)
                                                               cgcTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGc                                                                            <  8:85324/149‑1 (MQ=255)
                                                                       aaaaGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCAccc                                                                    >  7:195828/1‑149 (MQ=255)
                                                                                   ccATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCAGCTGCTgcgc                                                         >  6:317980/1‑148 (MQ=255)
                                                                                          cACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGt                                                 >  1:367845/1‑149 (MQ=255)
                                                                                            cGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTgg                                               <  2:295477/149‑1 (MQ=255)
                                                                                              gCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGcc                                             <  4:136902/149‑1 (MQ=255)
                                                                                               cATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCg                                            <  5:137786/149‑1 (MQ=255)
                                                                                                    aGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGCCCGTGGCCGGCAtg                                       <  4:261907/149‑1 (MQ=255)
                                                                                                              cATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCAcc                             >  3:220071/1‑149 (MQ=255)
                                                                                                              cATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGATCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCAcc                             >  8:285673/1‑149 (MQ=255)
                                                                                                                tGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTg                           <  8:195828/149‑1 (MQ=255)
                                                                                                                  gCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTgcg                         >  2:21001/1‑149 (MQ=255)
                                                                                                                   ccAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTgcgc                        <  3:286852/149‑1 (MQ=255)
                                                                                                                   ccAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTgcgc                        >  1:159151/1‑149 (MQ=255)
                                                                                                                          aTTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCgc                 >  3:153922/1‑149 (MQ=255)
                                                                                                                           ttCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCgct                >  2:268740/1‑149 (MQ=255)
                                                                                                                           ttCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCgct                >  5:228789/1‑149 (MQ=255)
                                                                                                                             cGCCGAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCgctg               <  4:220071/148‑1 (MQ=255)
                                                                                                                                 cAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCgct                >  6:151683/1‑143 (MQ=255)
                                                                                                                                 cAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCgct                <  5:151683/143‑1 (MQ=255)
                                                                                                                                         cggcaccggcaTGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCGCTGCTGCTGTTCTGcc  >  5:86333/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GTGGAAATCACCCTGGAGCATGGCGATGCGCTGCAGGCGATGGACCAGGCGGTGCGCTACAAGCGCTTGGTAAAAGGTGTGGCCCATGCCCACGGCATGCAGGCCTGCTTCATGGCCAAGCCATTCGCCCAACTGGCCGGCACCGGCATGCACATGCACCTGAGCCTGGCCGACAGCGCCGGCAACAACCTGTTCGCCAGCGAAGACAAGGCCGGCACCCCGCTGCTGCGCCAGGCCGTGGCCGGCATGCTGCGCCACCTGCGCGACTCGCTGCTGCTGTTCTGCC  >  NC_002947/5169022‑5169307

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: