Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,967,695 G→T 12.3% P18P (CCC→CCA amtB ← ammonium transporter AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,967,6950GT12.3% 105.3 / 3.8 39P18P (CCC→CCAamtBammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base G (18/16);  new base T (3/2);  total (21/18)
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.99e-01

GACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGC  >  NC_002947/5967546‑5967832
                                                                                                                                                      |                                                                                                                                         
gACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCg                                                                                                                                           <  2:186008/150‑1 (MQ=255)
   aGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGc                                                                                                                                        >  1:263686/1‑150 (MQ=255)
     cACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCggg                                                                                                                                         >  1:761309/1‑147 (MQ=255)
     cACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCggg                                                                                                                                         <  2:761309/147‑1 (MQ=255)
           cTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                                <  1:14620/150‑1 (MQ=255)
           cTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                                <  2:263686/150‑1 (MQ=255)
            ttGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAgg                                                                                                                               >  1:208236/1‑150 (MQ=255)
               gATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGCGGCACTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGAc                                                                                                                            >  1:225722/1‑150 (MQ=255)
                          gccgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTcc                                                                                                                 <  2:531003/150‑1 (MQ=255)
                               cATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGcc                                                                                                            >  2:757680/1‑150 (MQ=255)
                                 tAGAACAGGGCCAAGCCCGGAATAGTCCAGAACAGCTCCTACACCGTCGCTGTCAGCATTCAGGCGGTGGCGCCGATGTCCCACCTCGACCCTGCTGCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGCCACTAAGGCTCCTATCCCt                                                                                                          >  1:196048/1‑150 (MQ=255)
                                      cAGGGCCAAGGCCTGAATCGTCACGAACCGCCACAGCGACTCTGCTGTCAGTATTCGGTCCCTGCCGCCGGAATTAAGCGCTGTGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTAGGAGCCCtattgt                                                                                                     <  2:141302/150‑6 (MQ=255)
                                           ccAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTAc                                                                                                <  2:381184/150‑1 (MQ=255)
                                                 ccGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCagag                                                                                          >  2:562743/1‑150 (MQ=255)
                                                    gAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTAGGAGCCCtatt                                                                                                       >  2:670704/1‑131 (MQ=255)
                                                          tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTa                                                                                                 >  1:593111/1‑134 (MQ=255)
                                                          tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTa                                                                                                 <  2:593111/134‑1 (MQ=255)
                                                            aTGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGa                                                                                                                                 >  1:308834/1‑100 (MQ=255)
                                                            aTGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGa                                                                                                                                 <  2:308834/100‑1 (MQ=255)
                                                             tGAACAGCACCAACGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGtttt                                                                              >  2:463767/1‑150 (MQ=255)
                                                                     aCCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTggg                                                                      >  2:765207/1‑150 (MQ=255)
                                                                        aGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGCGCCGCCTGGCAGTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCg                                                                   >  2:457775/1‑150 (MQ=255)
                                                                          cgcCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGtt                                                                 <  2:208236/150‑1 (MQ=255)
                                                                                tGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGttt                                                           <  2:403695/150‑1 (MQ=255)
                                                                                    gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTCCGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTg                                                       >  2:246717/1‑150 (MQ=255)
                                                                                    gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTggg                                                                      >  2:802292/1‑135 (MQ=255)
                                                                                    gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTggg                                                                      <  1:802292/135‑1 (MQ=255)
                                                                                     tCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGa                                                      <  1:833554/150‑1 (MQ=255)
                                                                                      cAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAg                                                     <  1:528825/150‑1 (MQ=255)
                                                                                                         ccGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTggg                                                                      >  2:356739/1‑114 (MQ=255)
                                                                                                         ccGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTggg                                                                      <  1:356739/114‑1 (MQ=255)
                                                                                                                                 ttCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTTTTTTTTTTATCGCGTCGGTATCGGTTTCGCCgg           >  2:359743/1‑150 (MQ=255)
                                                                                                                                 ttCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCgg           >  2:588040/1‑150 (MQ=255)
                                                                                                                                       tGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGa     >  2:175706/1‑150 (MQ=255)
                                                                                                                                          cAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGc  <  2:196048/150‑1 (MQ=255)
                                                                                                                                          cAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGc  <  1:562743/150‑1 (MQ=255)
                                                                                                                                            gggcttgCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGt                   >  2:252801/7‑131 (MQ=255)
                                                                                                                                            gggcttgCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGt                   <  1:251307/125‑1 (MQ=255)
                                                                                                                                            gggcttgCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGt                   <  1:252801/125‑1 (MQ=255)
                                                                                                                                            gggcttgCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGt                   >  2:251307/7‑131 (MQ=255)
                                                                                                                                                  gCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgtt                                                                                   >  2:134168/1‑61 (MQ=255)
                                                                                                                                                  gCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgtt                                                                                   <  1:134168/61‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                         
GACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGC  >  NC_002947/5967546‑5967832

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