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.T94.7% 117.2 / ‑0.8 38coding (4/1500 nt)pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (15/21);  total (16/22)
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: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑T‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATAC  >  NC_000913/3637504‑3637775
                                                                                                                                              |                                                                                                                                   
aaTTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCTCGTTCATGTCCTCAACATGGCGTAACGTCCTATGC‑t                                                                                                                                    >  1:831663/1‑143 (MQ=255)
        aaGATAATTACTTCTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACAtttg                                                                                                                            <  1:66736/143‑1 (MQ=255)
             aaTTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCACAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGc                                                                                                                       >  1:317591/1‑143 (MQ=255)
             aaTTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGc                                                                                                                       <  1:797494/143‑1 (MQ=255)
              attattTTTGGAGTGGAATCCAGCCAGGGGGCAAATCAAAAAAAGTCTATATTGCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCt                                                                                                                      <  2:613283/143‑1 (MQ=255)
              atgagtttttGAACGAAAGCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCt                                                                                                                      <  2:294828/138‑1 (MQ=255)
                   tttttGATTGAAATTCATACAGGGGGCAAATCAAATAAAGTCTATATTTCTCTTTGCTCGCGCCGTGAAAGTCACTGATAATGCGCCGCGTTCATGTTCTCAAAATCTCGTAACGTCTTATGCTT‑TCATTTGTTTGTTGTCCt                                                                                                                 >  2:783098/1‑143 (MQ=255)
                                tCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATAcc                                                                                                    >  2:634829/1‑143 (MQ=255)
                                 ccATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCg                                                                                                   >  1:192944/1‑143 (MQ=255)
                                   aTACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCggg                                                                                                 <  1:497735/143‑1 (MQ=255)
                                    tACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGc                                                                                                >  1:397438/1‑143 (MQ=255)
                                         ggggCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGtta                                                                                           <  1:263109/143‑1 (MQ=255)
                                             gaaaGCAAAAAAGGTCTGTATGTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTtattat                                                                                       <  1:177006/142‑1 (MQ=255)
                                                ctCAAAAAAAGTCTGTATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgc                                                                                    <  1:731763/142‑1 (MQ=255)
                                                ataaaaaaaaGTCTATATTTCACTTTGCCCTCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgc                                                                                    <  2:832995/140‑1 (MQ=255)
                                                aTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgc                                                                                    <  2:827429/143‑1 (MQ=255)
                                                    aaaaaaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAACATGGCGTAACGTCCTATGCTT‑GCATGTGTTTGCTGGCCTGGATTGGCATACTGGGCTGGcattatgtctgct                                                                                >  1:321429/1‑136 (MQ=255)
                                                     caaaaGTCTATATTTCACTTGGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCa                                                                               <  2:116526/142‑1 (MQ=255)
                                                       aacgTCTATGTTTCACGTTGCCCGCGCCGCGAAAGTCGACGATACTGCGCCGCGTTCATGTCCTCACAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACt                                                                             <  2:308881/140‑1 (MQ=255)
                                                        aaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTg                                                                            <  1:794197/143‑1 (MQ=255)
                                                                 tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                                   >  1:797599/1‑143 (MQ=255)
                                                                     aCTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGt                                                               >  1:43471/1‑143 (MQ=255)
                                                                              gcgccgcgAAAGTCACTGATAAGGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑TAACATTAGTCTGCTGGCCTGGATTGGCATACCGGGCCGATATTATTACTTACCCTGGCTTTGGTGCTGTCCTACGAAg                                                      >  2:735401/1‑143 (MQ=255)
                                                                                  cgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCGTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCAt                                                  >  1:738474/1‑143 (MQ=255)
                                                                                    cgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCa                                                <  1:425364/143‑1 (MQ=255)
                                                                                          tCACTGATAATGCGCCGCGTTCATGGCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGtt                                          <  1:234040/143‑1 (MQ=255)
                                                                                            aCTGAGAATGCGCCGCGATCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTc                                        <  1:233268/143‑1 (MQ=255)
                                                                                             cTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCCAGGAAGCCAGCATTGGTTTcc                                       >  1:678190/1‑143 (MQ=255)
                                                                                              tGAGAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCa                                      <  2:724119/143‑1 (MQ=255)
                                                                                                   aTGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGaca                                 <  2:824165/143‑1 (MQ=255)
                                                                                                      cgccgcGTTCATGTCCTCAAAATGGCGTAACGGCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAg                              <  2:397438/143‑1 (MQ=255)
                                                                                                         cgcgTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCa                           <  2:831663/143‑1 (MQ=255)
                                                                                                            gTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGGTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACg                        >  1:445957/1‑143 (MQ=255)
                                                                                                                  gTCCTCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTgg                  >  1:680058/1‑143 (MQ=255)
                                                                                                                      tCAAAATGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGGTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGTCACAGCCAACGCCGTGGCAAc              >  1:720880/1‑143 (MQ=255)
                                                                                                                      tCAAAATGGCGTAACGTCCTATGCTT‑ACATGTGTTTGCTGGCCTGGATTTGCAGACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAAc              <  1:249548/143‑1 (MQ=255)
                                                                                                                        aaaaTGGCGTAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACCGCCACCGCCGTCGCAACCg            >  2:621690/1‑143 (MQ=255)
                                                                                                                                 tAACGTCCTATGCTT‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGCCAACCGTTATCtata   >  2:791208/1‑143 (MQ=255)
                                                                                                                                  aaCGTCCTATGC‑TGACGTTGGTTGGCTGGCCTGGAGTGGCATACCGGGCTGTTATTATTGCTTGCACTGGCTGTTGCGCTGTGCTACGAGGCCATCAAGGGTTGCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATAc  <  2:116980/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                   
AATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑T‑ACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATAC  >  NC_000913/3637504‑3637775

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