Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,403,310 T→A 9.4% intergenic (‑249/‑64) PP_2105 ← / → PP_2106 hypothetical protein/ammonium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,403,3100TA 9.4% 138.7 / 3.8 51intergenic (‑249/‑64)PP_2105/PP_2106hypothetical protein/ammonium transporter
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base A (2/3);  total (25/26)
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.40e-01

CGGAGCGCCTGTGACCGGTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATTCTCATGGGCGCCATTCTGGTGCTGG  >  NC_002947/2403165‑2403458
                                                                                                                                                 |                                                                                                                                                    
cGGAGCGCCTGTGACCGTTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCg                                                                                                                                                   <  3:68821/149‑1 (MQ=255)
   aGCGCCTGTGACCGGTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGACACCGATTCGGTg                                                                                                                                                <  8:122790/149‑1 (MQ=255)
      gcCTGTGACCGGTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAg                                                                                                                                             <  5:141895/149‑1 (MQ=255)
      gcCTGTGACCGGTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCT‑GCACCGAATCGGTGccag                                                                                                                                            <  6:315523/149‑4 (MQ=255)
       catgtgACCGGTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCTATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGAATCGGTGccag                                                                                                                                            <  8:301916/147‑4 (MQ=255)
        cTGTGACCGGTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAgg                                                                                                                                            >  3:67265/1‑148 (MQ=255)
        cTGTGACCGGTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGAATCGGTGccagc                                                                                                                                           <  5:68302/149‑5 (MQ=255)
                  ttATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTc                                                                                                                                 <  6:318253/149‑1 (MQ=255)
                      gcgcCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTg                                                                                                                              >  6:57055/1‑148 (MQ=255)
                      gcgcCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTg                                                                                                                              >  5:195856/1‑148 (MQ=255)
                        gcCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAAc                                                                                                                           >  8:266306/1‑149 (MQ=255)
                               cAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCt                                                                                                                     <  7:118400/148‑1 (MQ=255)
                                         gcgcAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGCACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCttttt                                                                                                          <  4:48776/149‑1 (MQ=255)
                                           gcAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCttttt                                                                                                          >  2:147504/1‑147 (MQ=255)
                                             aGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTcagc                                                                                                      <  7:266306/149‑1 (MQ=255)
                                                    aCCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTg                                                                                                                 <  2:124994/131‑1 (MQ=255)
                                                    aCCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTg                                                                                                                 >  1:124994/1‑131 (MQ=255)
                                                       gATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGcc                                                                                             >  1:199735/1‑148 (MQ=255)
                                                        aTAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTa                                                                                           >  8:90962/1‑149 (MQ=255)
                                                              cTTAACCCAATACACGACAAACCCAAGAAAAAACCACGCAGGCTTCACACCGCTAACCTGAACAGGGAAGGCGCTGGCACCGATTCGCTGCAGGGATGCAGGTTCTTGAACTGTCATCGTCCCATTTTCAGCAGGTCTGCCTATTTCAt                                                                                     >  4:40936/1‑149 (MQ=255)
                                                                 aaCCCAATCCACGACAAACCCAAGAAAAAACCTCGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATg                                                                                    <  3:40936/147‑1 (MQ=255)
                                                                 aaCCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATgg                                                                                   <  6:20505/148‑1 (MQ=255)
                                                                   cccAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGaaa                                                                                <  2:199735/149‑1 (MQ=255)
                                                                     cAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAAc                                                                              <  1:147504/149‑1 (MQ=255)
                                                                       aTCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACAt                                                                            <  5:279327/149‑1 (MQ=255)
                                                                        tCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAAAAAACCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATg                                                                           >  1:108713/1‑149 (MQ=255)
                                                                           aCGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGAATCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGca                                                                         >  5:200397/1‑148 (MQ=255)
                                                                             gACAAACCCAAGAAAAAACCACGCAGGCTTTACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAg                                                                      >  2:108109/1‑149 (MQ=255)
                                                                             gACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAg                                                                      >  8:113704/1‑149 (MQ=255)
                                                                               cAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCCTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAgcg                                                                    >  5:246589/1‑149 (MQ=255)
                                                                                   cccAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTAt                                                                >  5:126054/1‑149 (MQ=255)
                                                                                    ccAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATg                                                               <  6:195856/149‑1 (MQ=255)
                                                                                      aaGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTcag                                                                                                       >  7:80705/1‑107 (MQ=255)
                                                                                      aaGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTcag                                                                                                       <  8:80705/107‑1 (MQ=255)
                                                                                         aaaaaaCCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGActc                                                          >  3:200024/1‑149 (MQ=255)
                                                                                           aaaaCCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGActctc                                                        >  5:53707/1‑149 (MQ=255)
                                                                                             aaCCACGCAGGCTTTACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTg                                                      <  1:108109/149‑1 (MQ=255)
                                                                                                           cACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCTGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCCTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACAGTTcc                                          <  6:246589/147‑1 (MQ=255)
                                                                                                            aCCCCGCTCAGCTTAAGAGGGTAGGGGCTGGCACCGTTTCGGTGCATGGATTCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATTGAAAACATGCACAGCGCTATGGACTCTCCGGTCCACGGTTCCAac                                       <  3:95343/149‑1 (MQ=255)
                                                                                                            aCCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCaa                                        >  4:197107/1‑148 (MQ=255)
                                                                                                              cccGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAacac                                     >  6:194074/1‑149 (MQ=255)
                                                                                                               ccGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACAcc                                    >  2:296637/1‑149 (MQ=255)
                                                                                                                  cTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTg                                 <  3:110291/149‑1 (MQ=255)
                                                                                                                          aaCAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTcagca                                                                                                     <  2:236460/73‑1 (MQ=255)
                                                                                                                          aaCAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTcagca                                                                                                     >  1:236460/1‑73 (MQ=255)
                                                                                                                          aaCAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATtct                         >  6:228162/1‑149 (MQ=255)
                                                                                                                          aaCAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATtc                          <  6:126054/148‑1 (MQ=255)
                                                                                                                            cAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATTCTCa                       <  3:197107/149‑1 (MQ=255)
                                                                                                                                      ccctgCACCGAATCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATTCTCATGGGCGCCAt             >  2:4617/5‑149 (MQ=255)
                                                                                                                                          gCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATTCTCATGGGCGCCATTc           <  7:113704/147‑1 (MQ=255)
                                                                                                                                              cGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATTCTCATGGGCGCCATTCTGGTGc     >  8:148684/1‑149 (MQ=255)
                                                                                                                                                 ttCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATTCTCATGGGCGCCATTCTGGTGCTgg  >  3:301246/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGGAGCGCCTGTGACCGGTTATGCGCCGATTCAGGGCGCCCGCGCAGCGGATACCGATAAACCTTAACCCAATCCACGACAAACCCAAGAAAAAACCACGCAGGCTTCACCCCGCTCACCTGAACAGGGTAGGGGCTGGCACCGATTCGGTGCAGGGATGCAGGTTCTTGAACTGTCCTGGTGCCTTTTTCAGCAGGTCTGCCTATTTCATGGAAAACATGCACAGCGCTATGGACTCTCTGGTCCACGGTTCCAACACCCTGTTCATTCTCATGGGCGCCATTCTGGTGCTGG  >  NC_002947/2403165‑2403458

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