Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,967,628 A→C 7.7% L22R (CTG→CGG)  amtB ← ammonium transporter AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,967,6280AC 7.7% 137.3 / 4.2 52L22R (CTG→CGG) amtBammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (21/27);  new base C (2/2);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

CATCATCACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGG  >  NC_002947/5967480‑5967772
                                                                                                                                                    |                                                                                                                                                
catcatCACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCa                                                                                                                                                  <  4:152069/149‑1 (MQ=255)
  tcatcaCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAgg                                                                                                                                                <  8:7787/149‑1 (MQ=255)
  tcatcaCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAgg                                                                                                                                                >  5:280611/1‑149 (MQ=255)
         cGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGcc                                                                                                                                         <  2:26667/149‑1 (MQ=255)
           aCAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGc                                                                                                                                          <  6:299877/146‑1 (MQ=255)
           aCAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGc                                                                                                                                          >  5:299877/1‑146 (MQ=255)
             acgaCGTGCTTGGATCGCACCATGCCGCCATAGAACAGGCCCAATCCCGGAATAGTCATGAACAGCACCAGCGCCGTGGCGGTCAGCATCCAGGGGGTGTCGACGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGCTAGCCCtgg                                                                                                                                     <  7:142300/146‑1 (MQ=255)
             aGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCggg                                                                                                                                     <  4:312283/149‑1 (MQ=255)
               cACGTTCTTGGTTCGCACCATGCCGCCATAGATCAGGGCCAAGCCCGGAATTGTCATGAACAGCTCCAGCGCCGTTGCGGTCAGCATCCAGGGGGTGTCGCCGGATTTAAGCGCTGGGGCAGCTTCCTGTGCCAGTGCAAGCCCGGGCa                                                                                                                                   <  2:18212/149‑1 (MQ=255)
               cACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCa                                                                                                                                   >  1:326331/1‑149 (MQ=255)
                     cTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCAGCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACg                                                                                                                              <  5:267229/148‑1 (MQ=255)
                      ttGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                            >  3:33060/1‑149 (MQ=255)
                        ggATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGa                                                                                                                             <  5:332365/146‑1 (MQ=255)
                        ggATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGa                                                                                                                             <  6:280611/146‑1 (MQ=255)
                        ggATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                            <  6:230806/147‑1 (MQ=255)
                        ggATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                            <  5:296501/147‑1 (MQ=255)
                         gATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGa                                                                                                                         <  2:326331/149‑1 (MQ=255)
                            cGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACaa                                                                                                                      >  5:179275/1‑149 (MQ=255)
                             gCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAgg                                                                                                                           >  2:102277/1‑143 (MQ=255)
                             gCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAgg                                                                                                                           <  1:102277/143‑1 (MQ=255)
                                     ccgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTcc                                                                                                             <  4:33060/149‑1 (MQ=255)
                                     ccgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTc                                                                                                              <  6:179275/148‑1 (MQ=255)
                                      cgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTcc                                                                                                             >  5:183930/1‑148 (MQ=255)
                                             gAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGCTTGCCCTGGCaga                                                                                                                                 >  2:234245/1‑119 (MQ=255)
                                             gAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGCTTGCCCTGGCaga                                                                                                                                 <  1:234245/121‑3 (MQ=255)
                                             gAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCt                                                                                                      <  6:14086/148‑1 (MQ=255)
                                              aaCAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGc                                                                                                    >  7:312172/1‑149 (MQ=255)
                                                           ccGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCa                                                                                         >  6:95184/1‑147 (MQ=255)
                                                                  aGTCATGAACAGCACCAGCGCAGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgt                                                                                >  8:47413/1‑149 (MQ=255)
                                                                    tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgttg                                                                              >  4:239598/1‑149 (MQ=255)
                                                                    tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgttg                                                                              >  2:276441/1‑149 (MQ=255)
                                                                         aaCAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtc                                                                         >  4:265727/1‑149 (MQ=255)
                                                                         aaCAGCACCAGCGCCGTTGCCGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATTGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtc                                                                         <  4:292314/149‑1 (MQ=255)
                                                                           cAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGTGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGTGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtctc                                                                       <  6:25772/149‑1 (MQ=255)
                                                                           cAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtctc                                                                       <  5:95184/149‑1 (MQ=255)
                                                                              cACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCt                                                                                                      <  2:279219/115‑1 (MQ=255)
                                                                              cACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCt                                                                                                      >  1:279219/1‑115 (MQ=255)
                                                                              cACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTg                                                                    >  1:191892/1‑149 (MQ=255)
                                                                                     gcCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGtt                                                             <  3:239598/149‑1 (MQ=255)
                                                                                              gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGt                                                    <  6:183930/149‑1 (MQ=255)
                                                                                               tCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTg                                                   >  4:125019/1‑149 (MQ=255)
                                                                                               tCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTg                                                   <  2:261967/149‑1 (MQ=255)
                                                                                                    aTCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCt                                              >  4:163915/1‑149 (MQ=255)
                                                                                                        aGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCttttt                                          >  2:6476/1‑149 (MQ=255)
                                                                                                            ggTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTtgct                                      <  3:163915/149‑1 (MQ=255)
                                                                                                                                cTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATc                  <  8:312172/149‑1 (MQ=255)
                                                                                                                                 tGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCg                 >  3:280077/1‑149 (MQ=255)
                                                                                                                                   gggCAGCTTCCTCTGCCAGGGCCAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGt               >  8:269430/1‑149 (MQ=255)
                                                                                                                                        gCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGc          >  7:152392/1‑149 (MQ=255)
                                                                                                                                            ccTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCgg       <  6:287376/148‑1 (MQ=255)
                                                                                                                                             ctctGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTa     <  1:6476/149‑1 (MQ=255)
                                                                                                                                             ctctGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTa     <  8:152392/149‑1 (MQ=255)
                                                                                                                                                tGCCCGGGCTTGCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACgg  >  1:39563/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CATCATCACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGG  >  NC_002947/5967480‑5967772

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