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.T99.5% 103.4 / ‑7.0 37L2L (CTA→TTA) pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base T (18/17);  minor base A (1/0);  total (19/18)
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: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCC  >  NC_000913/3637511‑3637777
                                                                                                                                       |                                                                                                                                     
tacgATAATGATTTTTTGAGTGGACACCGTACCGGGGGCAAAGCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAttt                                                                                                                               <  1:59482/140‑1 (MQ=255)
   gATAATTATTTTTTGAGTGAAATCCATACAGGGGGCGAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttgtt                                                                                                                            >  1:528799/1‑143 (MQ=255)
                   gTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTg                                                                                                            <  1:538253/143‑1 (MQ=255)
                   gTGAAATCCATACAGGCGGCAAGTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGTGTATCGCCCTCTGCTTACATTTGTTTGCTGGCCTGGACTTg                                                                                                            <  2:714539/143‑1 (MQ=255)
                            aTACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCggg                                                                                                   <  1:493047/143‑1 (MQ=255)
                              aCAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCt                                                                                                 <  1:79305/143‑1 (MQ=255)
                               cAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTg                                                                                                >  1:538025/1‑143 (MQ=255)
                                  ggggCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCCTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATAGCGGGCCGtca                                                                                             >  2:409846/1‑141 (MQ=255)
                                     gCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATTCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTCGTTTGCTGGCCTGGATTTGCATACCGGGCTGttatta                                                                                          >  1:678045/1‑143 (MQ=255)
                                     gCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGttatta                                                                                          <  1:300392/143‑1 (MQ=255)
                                       aaaTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTCTGCTTACATTTGTTTGCTGGCCTGGATTTGCATATCGCCCTGTTactatt                                                                                        >  1:208654/1‑143 (MQ=255)
                                                      tataTTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCttt                                                                         >  1:629109/1‑143 (MQ=255)
                                                          tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTACTTGCACTGGCttttctt                                                                     >  2:313921/1‑140 (MQ=255)
                                                          tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCACGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCTTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                                     <  1:323832/143‑1 (MQ=255)
                                                          tttCACTTTGCCCGCGCCGCGAAAGGCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGACCTAAGCAAACATTTGTTTGCTGGCCTGGATTTGCATCCCGGGCTGGTCTTATTGCTTGCACTGGCTTTtgtg                                                                     >  1:812027/1‑143 (MQ=255)
                                                            tCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCt                                                                   >  2:241822/1‑143 (MQ=255)
                                                            tCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCt                                                                   <  1:317593/143‑1 (MQ=255)
                                                             cACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTg                                                                  <  1:623095/143‑1 (MQ=255)
                                                               cTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCGCAAAATGGCGTAACGTCCTATGCTTACATTTGTTGGCTGGCCTGGATTTGCATACGGGGCTGTTATTATTGCTTGCCCTGGCTTTTGTGCTggt                                                                >  1:669261/1‑141 (MQ=255)
                                                                 ttGCCCGCGCCGACAAAGTCACTCTTAATGCGCCGCGTTCACGTCCTCAAAATGGCGTAACGTCCTATTCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCt                                                              <  1:145337/143‑1 (MQ=255)
                                                                     ccgcgccgcgAAAGTCACTGATAATGCGCCGCGTTCCTGTCCTAAAAATGGCGTAACGTACTATGCTTGCATTTGTTTGCTGGCCTGGATTTGCATACCGGGATAGTATTATTGCTTGGACTTGCTTTTGTCCTTTTCTACGa                                                          >  2:438653/1‑143 (MQ=255)
                                                                     ccgcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTCCGa                                                          >  1:683538/1‑143 (MQ=255)
                                                                     ccgcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGGTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGa                                                          >  1:613529/1‑143 (MQ=255)
                                                                      cgcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGaa                                                         <  1:601150/143‑1 (MQ=255)
                                                                       gcgccgcgAAAGTCACTGATAATGAGAAGCGTTACGGTGCTCAAAATGGCGTAACGGCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACggcg                                                        >  1:346263/1‑140 (MQ=255)
                                                                           cgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCAt                                                    <  1:703569/143‑1 (MQ=255)
                                                                           cgcgAAAGCCCCTGATAATGCGCCGCGTGCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCAt                                                    <  1:613916/143‑1 (MQ=255)
                                                                                   tCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGCTTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGACGCCATCAATgggt                                            >  1:44571/1‑141 (MQ=255)
                                                                                        gATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCGTCAATGGTTTCCAt                                       >  1:314615/1‑143 (MQ=255)
                                                                                          tAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGa                                     <  1:103696/143‑1 (MQ=255)
                                                                                                         aTGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTATACGAAGACATCAATTGTTTCCATGACACAGCCAACGCCGt                      >  1:551146/1‑143 (MQ=255)
                                                                                                            tCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGc                   >  1:253625/1‑143 (MQ=255)
                                                                                                              cgcAAACTGGCGTAACGTCCTATGCTTGCATTTGTTTGCTGGCCTGTATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCaa                 <  2:369045/141‑1 (MQ=255)
                                                                                                                 aaaaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCg              <  1:460134/143‑1 (MQ=255)
                                                                                                                 aaaaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCg              >  2:87401/1‑143 (MQ=255)
                                                                                                                         gTAGCGTCCTATGCTGACATTTGTTTGCTGGCCTGGCTTTGCGTACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCtat      <  2:394805/143‑1 (MQ=255)
                                                                                                                             cGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATAccc  <  2:310059/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
TAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCC  >  NC_000913/3637511‑3637777

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: