Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,637,646 +T coding (5/1500 nt) pitA → metal phosphate:H(+) symporter PitA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,637,6451.T96.4% 168.7 / ‑0.3 57coding (4/1500 nt)pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (22/32);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCCGCGC  >  NC_000913/3637504‑3637781
                                                                                                                                               |                                                                                                                                        
aaTTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCAGGGCCTCAAAGTGGCGTACAGTGCTATGCt                                                                                                                                         >  1:1094078/1‑143 (MQ=255)
      agaaGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtt                                                                                                                                    <  1:254400/141‑1 (MQ=255)
         aGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttgt                                                                                                                                 <  1:776934/143‑1 (MQ=255)
           aTAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTCCGTCGCCTCAGAGTGGCGTACCGTCCTATGCTTACAtttgttt                                                                                                                               >  1:860778/1‑143 (MQ=255)
            tAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttgtttg                                                                                                                              <  1:1132848/143‑1 (MQ=255)
                 attTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGc                                                                                                                         >  1:719056/1‑143 (MQ=255)
                  ttttttGAGTGAAATCCATACAGTCGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGcc                                                                                                                        <  2:823909/143‑1 (MQ=255)
                   tttttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTTCCCGTTCCGTGTAATTCACTGTTAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCt                                                                                                                       <  1:399755/143‑1 (MQ=255)
                    ttttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTg                                                                                                                      >  2:254714/1‑143 (MQ=255)
                    ttttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTg                                                                                                                      >  1:965817/1‑143 (MQ=255)
                      ttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCGATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGa                                                                                                                    >  1:802316/1‑143 (MQ=255)
                           tGAAATCCATTCAGTGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGc                                                                                                               <  2:245206/143‑1 (MQ=255)
                           tGAAATCCATCCAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGCTCATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTACTATGCTTACACATGTTTGCTGGCCTGCATATGc                                                                                                               >  2:864761/1‑143 (MQ=255)
                           tGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGc                                                                                                               <  2:395478/143‑1 (MQ=255)
                           tGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTAGCCTGGCTTTGc                                                                                                               >  2:981366/1‑143 (MQ=255)
                            gAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCTGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCa                                                                                                              <  1:800101/143‑1 (MQ=255)
                               atccatccAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCACGCGCCGCAAAAGAAACTGCTAATGCGCCGCGTTCATGTAATCAAAATGGCGTAAAGTCCTATGCTTACATTTGTTTGCTGGCCTGGCTTTGCATAc                                                                                                           >  1:840925/1‑143 (MQ=255)
                                         ggggCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCATTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGtta                                                                                                 <  1:298626/143‑1 (MQ=255)
                                            gCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCATCAAAATGGCGTAACGTCCTATGCTTACATTGGTTTGCTGGCCTGGCTTTGCATACCGGGCTGttcttc                                                                                              >  1:1163360/1‑142 (MQ=255)
                                               aaTCAAAAAAAGTCTATATTTCACTTCGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttg                                                                                           <  2:208695/143‑1 (MQ=255)
                                                aTCAAAAAAAGTCTATATTTCACTTGGGCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCGCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgct                                                                                         <  1:828736/143‑1 (MQ=255)
                                                 tCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCGTGGATTTGCATACCGGGCTGTTATTAttgct                                                                                         >  1:130313/1‑143 (MQ=255)
                                                  cAAAAAAAGTCTAGAGTTCACTTTGCCCGCGCCGCGAAAGTCGCTGATAGGGCGCCGCGTTCATGTCCTCAAAATGGCGTATCGTCCTATGCTTACATTTGTTTGCTGGCCTGGTTTTGCATACCGGGCTGTTATTAttgctt                                                                                        <  2:370605/143‑1 (MQ=255)
                                                    aaaaaaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgcttgc                                                                                      >  2:1125792/1‑143 (MQ=255)
                                                         aGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTgg                                                                                 <  1:532718/143‑1 (MQ=255)
                                                          gTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGc                                                                                >  1:1164224/1‑143 (MQ=255)
                                                           tCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGCTATTATTGCTTGCACTGGCt                                                                               >  2:346133/1‑143 (MQ=255)
                                                            cTATATTTCACTTTGCCTGCCCCGCGAAATTCACTGATACTGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCTGGCTGTTATTATTGCTTGCACTGGCtt                                                                              <  1:204745/143‑1 (MQ=255)
                                                            cTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCCATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtt                                                                              <  1:165797/143‑1 (MQ=255)
                                                              atatTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtttt                                                                            <  1:836289/143‑1 (MQ=255)
                                                                atTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgt                                                                          <  1:881683/143‑1 (MQ=255)
                                                                 tttCACTTTGCCCGCGCCGCGAACGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                                         <  2:421572/143‑1 (MQ=255)
                                                                 gttgaCTTTGCTTGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                                         <  2:414279/139‑1 (MQ=255)
                                                                 cttCACTTTGCGCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCCTGTCCTCAAAATGGCGTAACGTCCTATGCGTCCATTCGTCTGCTGGCCTGGATTTGCCTACCGGGCGGTTATTCGTGCTTGAACCGGCcttgttg                                                                         >  2:832260/2‑139 (MQ=255)
                                                                    cACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTg                                                                      <  1:718943/143‑1 (MQ=255)
                                                                           cccGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACg                                                               >  1:108956/1‑143 (MQ=255)
                                                                           cccGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATGATTGCTTCCACTGGCTTTTGTGCTGGTcggtg                                                               >  2:304806/1‑139 (MQ=255)
                                                                            ccgcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGa                                                              >  1:1070547/1‑143 (MQ=255)
                                                                                 ccgcgAACGTCCCTGATATTGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCa                                                         <  2:460854/143‑1 (MQ=255)
                                                                                      aaaGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAAt                                                    <  1:623582/143‑1 (MQ=255)
                                                                                            acggATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTc                                              <  1:506813/140‑1 (MQ=255)
                                                                                             cTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTcc                                             <  1:453319/143‑1 (MQ=255)
                                                                                               gATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCAt                                           <  1:96598/143‑1 (MQ=255)
                                                                                                 tgaTGCTCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGa                                         <  1:522044/141‑1 (MQ=255)
                                                                                                 tAATGCGCCGCGTTCATGTCCTCCAAATGGCGTAACGTCCTATGCTTACATTTGTTTGTTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGa                                         <  1:101730/143‑1 (MQ=255)
                                                                                                       gccgcGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGc                                   >  1:319152/1‑143 (MQ=255)
                                                                                                       gccgcGGTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTTTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGc                                   <  1:350998/143‑1 (MQ=255)
                                                                                                         cgcgTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCa                                 <  2:108956/143‑1 (MQ=255)
                                                                                                          gcgTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCaa                                <  1:254714/143‑1 (MQ=255)
                                                                                                                aTGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGt                          >  1:372336/1‑143 (MQ=255)
                                                                                                                 tGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGTCGTCAATTGTTGCCATGACACAGCCAGCGCCGTg                         >  2:976295/1‑143 (MQ=255)
                                                                                                                     cTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCaa                     >  1:754355/1‑143 (MQ=255)
                                                                                                                      tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAAc                    >  1:319228/1‑143 (MQ=255)
                                                                                                                           cgcgCGTACCGTTCTATGGTTTCATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACCAAGCCATCAATGTTTTCCATGACACAGCCAACGCCGTGGCAACCGTTa               <  2:821848/140‑1 (MQ=255)
                                                                                                                            tGGCGTAACGCCCTATGCTTACATTTGTTTGCTGGCCTGGCTTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCAGGACACAGCCAACGCCGTGGCAACCGTTAt              <  2:102818/143‑1 (MQ=255)
                                                                                                                                     gTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCcg     >  2:732886/1‑143 (MQ=255)
                                                                                                                                      tcttATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTGGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCcgc    <  2:693363/140‑1 (MQ=255)
                                                                                                                                        cTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCACCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCcgcgc  <  1:180073/143‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                        
AATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCCGCGC  >  NC_000913/3637504‑3637781

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