Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,169,742 G→T 7.5% A454S (GCC→TCC)  PP_4547 → glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,169,7420GT 7.5% 166.0 / 2.9 53A454S (GCC→TCC) PP_4547glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (19/30);  new base T (2/2);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

CACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTGGATGTGG  >  NC_002947/5169596‑5169886
                                                                                                                                                  |                                                                                                                                                
cACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGcc                                                                                                                                                <  7:180343/149‑1 (MQ=255)
 aCTCGACGCCATGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGCCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCt                                                                                                                                               <  5:61192/149‑1 (MQ=255)
     gACGCCCTGCAACATTCCGATTGGGCCCGGGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGa                                                                                                                                             <  7:118429/147‑1 (MQ=255)
       cGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGccc                                                                                                                                         <  6:155847/149‑1 (MQ=255)
            tGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTaaa                                                                                                                                    >  8:125073/1‑149 (MQ=255)
              cAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAt                                                                                                                                   <  1:286835/148‑1 (MQ=255)
                aCATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCCAAATaa                                                                                                                                 >  2:251669/1‑148 (MQ=255)
                  aTTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACaa                                                                                                                              <  2:52442/149‑1 (MQ=255)
                  aTTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCCAAATAACaa                                                                                                                              <  1:251669/149‑1 (MQ=255)
                   ttCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAg                                                                                                                             <  8:274336/149‑1 (MQ=255)
                     ccGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTATCTGCACCAGGCCTGAGCCCCTAAATAACAAgg                                                                                                                            <  2:269674/148‑1 (MQ=255)
                     ccGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAgg                                                                                                                            >  3:27946/1‑148 (MQ=255)
                        aTTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCt                                                                                                                                               >  5:282331/1‑126 (MQ=255)
                        aTTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCt                                                                                                                                               <  6:282331/126‑1 (MQ=255)
                        aTTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAaca                                                                                                                        <  6:79289/149‑1 (MQ=255)
                            ggCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACAccc                                                                                                                     <  2:23098/148‑1 (MQ=255)
                            ggCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCa                                                                                                                    <  3:239270/149‑1 (MQ=255)
                            ggCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCa                                                                                                                    >  5:150780/1‑149 (MQ=255)
                             gCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCAt                                                                                                                   <  3:305123/149‑1 (MQ=255)
                              cccGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATg                                                                                                                  <  5:291478/149‑1 (MQ=255)
                                 gTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAt                                                                                                                                   <  4:77148/129‑1 (MQ=255)
                                 gTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAt                                                                                                                                   >  3:77148/1‑129 (MQ=255)
                                  tGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAGTCCt                                                                                                                                               >  7:294171/1‑116 (MQ=38)
                                  tGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAGTCCt                                                                                                                                               <  8:294171/116‑1 (MQ=38)
                                       cATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACgcagca                                                                                                         >  1:217458/1‑149 (MQ=255)
                                           gggCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTg                                                                                                      >  8:251800/1‑148 (MQ=255)
                                              cGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACgca                                                                                                            >  7:249440/1‑139 (MQ=255)
                                              cGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACgca                                                                                                            <  8:249440/139‑1 (MQ=255)
                                               gAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTg                                                                                                      <  7:195608/144‑1 (MQ=255)
                                               gAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTg                                                                                                      >  8:195608/1‑144 (MQ=255)
                                                 aGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACg                                                                                               <  7:55702/149‑1 (MQ=255)
                                                     ttCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAGTCCt                                                                                                                                               >  7:210606/1‑97 (MQ=37)
                                                     ttCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAGTCCt                                                                                                                                               <  8:210606/97‑1 (MQ=37)
                                                           ggCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGgcg                                                                                       <  7:251800/147‑1 (MQ=255)
                                                              gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACgcagc                                                                                                          <  6:218586/125‑1 (MQ=255)
                                                              gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACgcagc                                                                                                          >  5:218586/1‑125 (MQ=255)
                                                              gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACg                                                                                               >  2:110338/1‑136 (MQ=255)
                                                              gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACg                                                                                               <  1:110338/136‑1 (MQ=255)
                                                                    cTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGcc                                                                            >  5:84060/1‑149 (MQ=255)
                                                                      gAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACaa                                                                                                                              >  6:229627/1‑97 (MQ=255)
                                                                      gAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACaa                                                                                                                              <  5:229627/97‑1 (MQ=255)
                                                                       aaGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTc                                                                         >  1:266100/1‑149 (MQ=255)
                                                                                  cgcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAgcg                                                              <  7:125073/149‑1 (MQ=255)
                                                                                   gcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCTAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCAt                                                                                                        <  3:205503/106‑1 (MQ=255)
                                                                                   gcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCTAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCAt                                                                                                        >  4:205503/1‑106 (MQ=255)
                                                                                         tACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCt                                                       <  6:150780/149‑1 (MQ=255)
                                                                                                           gAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTa                                     <  6:84060/149‑1 (MQ=255)
                                                                                                                  gCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCaa                              <  2:217458/149‑1 (MQ=255)
                                                                                                                          gACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAggg                      >  4:301410/1‑149 (MQ=255)
                                                                                                                            cTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTc                                                                         <  5:208466/96‑1 (MQ=255)
                                                                                                                            cTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTc                                                                         >  6:208466/1‑96 (MQ=255)
                                                                                                                                             aCCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTGGAtgtg   >  5:67714/1‑149 (MQ=255)
                                                                                                                                              ccAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTGGATGTgg  <  6:67714/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTGGATGTGG  >  NC_002947/5169596‑5169886

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: ‑