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.T100.0% 89.1 / NA 31L2L (CTA→TTA) pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (19/12);  total (19/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑TACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCC  >  NC_000913/3637503‑3637777
                                                                                                                                               |                                                                                                                                    
gAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGc                                                                                                                                       <  2:263235/143‑1 (MQ=255)
    aaaataAGATAATTATTTTTTGAATACAATCCATACACGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑ttac                                                                                                                                  <  2:622159/140‑4 (MQ=255)
        taAGATAATTATTTTTTGATTGAAATCCATACAGGTGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGAGCCACGAAAGTCACTGATAATGCGCCACGTTCATGTCCCGAAAATGGCGTAACGTCCTCTGCTTACAttt                                                                                                                               >  2:267641/1‑143 (MQ=255)
                          aGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGAttt                                                                                                             >  2:303543/1‑143 (MQ=255)
                                aTCCATACAGGTGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATAc                                                                                                       <  2:894676/143‑1 (MQ=255)
                                aTCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATAc                                                                                                       >  2:374817/1‑143 (MQ=255)
                                aTCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTGGCATAc                                                                                                       >  2:385479/1‑143 (MQ=255)
                                                   cAAAAAATGTCTATATTTCTCTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgctt                                                                                    <  1:116024/143‑1 (MQ=255)
                                                   cAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAGTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgctt                                                                                    >  1:78332/1‑143 (MQ=255)
                                                       aaaaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACCTTTGTTTGCTGGCCTGGATTTGCATACCGGCCCGTTCTTATTGCCTGCAc                                                                                >  2:288227/1‑143 (MQ=255)
                                                              tataTTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCttt                                                                         >  2:626348/1‑143 (MQ=255)
                                                               atatTTCACTTTCCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGATCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtttt                                                                        <  2:341142/143‑1 (MQ=255)
                                                                 atTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgt                                                                      <  1:303543/143‑1 (MQ=255)
                                                                  tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGCCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTCTtgtg                                                                     <  2:742018/143‑1 (MQ=255)
                                                                               gcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAg                                                        <  2:776900/143‑1 (MQ=255)
                                                                                     cgAAAGTCAGTGAGAGTGCGTCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCCGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCa                                                  <  1:727532/143‑1 (MQ=255)
                                                                                     cgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTGCGAAGCCATcc                                                  >  1:49256/1‑142 (MQ=255)
                                                                                       aaaGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTGCATTTGTTTGCTGGCATGGATTTGCCTACCCGGCTGTTATTATTGCTTGCACTAGCCTTTGTGCTGTTTCACGAAGCCaaccag                                                >  2:895005/1‑138 (MQ=255)
                                                                                           tCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCCATGGtt                                            >  1:152631/1‑143 (MQ=255)
                                                                                              cTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATCTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCGGTTCCACGAAGCCATCAATGGTTGcc                                         >  2:533442/1‑143 (MQ=255)
                                                                                               tGATAATGCGCCGCGTTCATGTCCTCAAAATGTCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCa                                        <  1:450082/143‑1 (MQ=255)
                                                                                               tGATAATGCGCCGAGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTGGTTTGCTGGCCTGGATTTGCATACCGGGCGGTTATTATTGCTTGCACTGGCTTCTGTGCCGTTCTGCGAAGCCCTCCCGGGTTTCCa                                        >  1:180319/1‑143 (MQ=255)
                                                                                                gATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCAt                                       >  2:432869/1‑143 (MQ=255)
                                                                                                              ttCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCCTCAATGGTTTCCATGACACAGTCAACGc                         >  1:286543/1‑143 (MQ=255)
                                                                                                               tCATGTCCCCGAAATGGCGTACCGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACCGCCaaatcg                        <  2:596073/143‑5 (MQ=255)
                                                                                                                  tGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTg                     >  2:231174/1‑143 (MQ=255)
                                                                                                                  tGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTg                     <  2:54075/143‑1 (MQ=255)
                                                                                                                       tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAAc                >  1:246042/1‑143 (MQ=255)
                                                                                                                       tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGCTTCCATGACACAGCCAACGCCGTGGCAAc                >  2:48105/1‑143 (MQ=255)
                                                                                                                       tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGGTCTACGAAGCCATCAATGGGTTCCATGACACAGCCACCGCCGTGGCacc                >  1:352519/1‑141 (MQ=255)
                                                                                                                        cAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAAcc               <  1:809758/143‑1 (MQ=255)
                                                                                                                               gCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGGTGGTCTACGAAGCCATCAATGGTTTCCATGAGACAGCCAACGCCGTGGCAACCGTTATCt        >  1:78681/1‑143 (MQ=255)
                                                                                                                                     cGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATAccc  <  1:440274/143‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                    
GAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑TACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCC  >  NC_000913/3637503‑3637777

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