Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,637,645:1 +T 100% coding (4/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.T89.2% 92.6 / 5.3 37coding (4/1500 nt)pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (13/19);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATA  >  NC_000913/3637503‑3637774
                                                                                                                                                |                                                                                                                                   
gAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTCACGTCCTATGc                                                                                                                                      >  1:597371/1‑143 (MQ=255)
         aaGATAATTATTTTTTGAGTAAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTATGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttg                                                                                                                           <  1:263769/143‑1 (MQ=255)
               attattTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCt                                                                                                                     <  2:304810/143‑1 (MQ=255)
                ttattTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTg                                                                                                                    <  1:365457/143‑1 (MQ=255)
                    tttttGAGTGAAATCCAGACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCt                                                                                                                <  2:426337/143‑1 (MQ=255)
                     ttttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTCTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGCCCTCAAAATGGCGTAACGTCCTATGCCTACACTTGTTCGCTGGCCTg                                                                                                               >  2:85011/1‑143 (MQ=255)
                       ttGAGTAAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGa                                                                                                             <  2:232468/143‑1 (MQ=255)
                          aGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATGTGTTTGCTGGCCTGGAttt                                                                                                          <  1:526199/143‑1 (MQ=255)
                                    atgcggggggCGACTCAAACAAGGTCTCTATTTCGCTTTGTCCGCGCCGCGAAATGCACTGATAATGCGCCGCGTTCATGTCCGCAAAATGGCGTAGAGTCGTGTCCAACAATATATGGGCGGGCCTGGCTTTTCATACCggg                                                                                                <  2:1084532/138‑1 (MQ=255)
                                     tACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATACTGCGCCGCGTTCATGTCCTCAAAATGGCGCAACGTCCTATGCTTACATTTGCTTGCTGGCCTGCATTTGCATACGGGGc                                                                                               >  2:700727/1‑143 (MQ=255)
                                     tACAGGGGGCAAATCAAAAAAAGTCTATAGTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGc                                                                                               <  1:26820/143‑1 (MQ=255)
                                        aGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCGATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGt                                                                                            >  2:326171/1‑143 (MQ=255)
                                             gCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGttatta                                                                                       >  1:765642/1‑143 (MQ=255)
                                               aaaTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTattatt                                                                                     <  2:375217/143‑1 (MQ=255)
                                                  tCAAAAAACGTCTATATTTCACTTTGCCCGCGCCGCGCAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTTCTGACCCGGATTCACCTACCGCACTGTTAGTAttgct                                                                                  >  2:760686/1‑143 (MQ=255)
                                                                  tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                                  >  2:682606/1‑143 (MQ=255)
                                                                  tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCAGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGTCTTtgtg                                                                  >  2:472978/1‑143 (MQ=255)
                                                                     cacgtggcccgcccCGCGAAAGTCGCTTGTAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTg                                                               <  2:452635/137‑1 (MQ=255)
                                                                     cACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTg                                                               <  1:734658/143‑1 (MQ=255)
                                                                          tGCCCGCGCCGCGAACGTCACGGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTGTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTa                                                          <  2:1005851/143‑1 (MQ=255)
                                                                              cgcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTGCTGCGaa                                                      >  2:573438/1‑143 (MQ=255)
                                                                                  ccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCGATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCa                                                  <  1:326171/143‑1 (MQ=255)
                                                                                      caaaGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCaa                                              <  1:502636/142‑1 (MQ=255)
                                                                                             aCTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTc                                       >  1:710590/1‑143 (MQ=255)
                                                                                             aCTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGCGCTGGTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTc                                       >  2:956904/1‑143 (MQ=255)
                                                                                              cTGATAATGCGCCTCCTTCATGTCCTAAAAATTGCTCAACGTCCTATGTTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTcc                                      <  2:744800/143‑1 (MQ=255)
                                                                                                gATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCAt                                    <  1:67244/143‑1 (MQ=255)
                                                                                                gATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCAt                                    <  2:573487/143‑1 (MQ=255)
                                                                                                  tAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCCATCCGCACATTGGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGa                                  <  2:597212/143‑1 (MQ=255)
                                                                                                        gccgcGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTTTTATTTCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGc                            <  1:466426/143‑1 (MQ=255)
                                                                                                        gccgcGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCCATGGTTTCCATGACACAGc                            >  1:3215/1‑143 (MQ=255)
                                                                                                                   gTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTgg                 <  1:584092/143‑1 (MQ=255)
                                                                                                                      cTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCaa              >  1:963960/1‑143 (MQ=255)
                                                                                                                       tCAAAATGGCGTAACGACCTACGCTGACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCCTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAAc             <  1:245723/143‑1 (MQ=255)
                                                                                                                           aaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCCTGACTGGTTTCCATGCCACAGCCAACGCCGTGGCCACCGtt         >  2:286275/1‑143 (MQ=255)
                                                                                                                           aaTGGCGTAACATCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCACTGGTTTCCATGACACGGCCAACGTCGTGGCGACCGtt         >  1:844575/1‑143 (MQ=255)
                                                                                                                             tGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCACCCGGGGGATGAATCATGGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTAt       <  1:1097869/143‑1 (MQ=255)
                                                                                                                                  tAACGTCATATGCTTACATTTGTTTGCTGGCCTGGATTTGCCTACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCtata  <  1:536299/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                   
GAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATA  >  NC_000913/3637503‑3637774

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: