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.T87.8% 105.0 / 10.0 41coding (4/1500 nt)pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base T (15/21);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
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.

TATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATC  >  NC_000913/3637509‑3637770
                                                                                                                                          |                                                                                                                              
tataAGATAATTATTTTTTGAGTGACATCCATCCAGGGGGCAAACCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCttacat                                                                                                                          <  1:551161/143‑6 (MQ=255)
tataAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAACTGGCGTAACGTCCTATGCttacat                                                                                                                          >  2:285522/1‑138 (MQ=255)
     gATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttgtt                                                                                                                     >  1:625765/1‑143 (MQ=255)
      aTAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAGGCCTATATTCCACTTTGCCCGCGCCGCGAAAGTCCCGGATAATGGGCCGCGTCCATGCCCCCAAAAGGGCGAACGGCCCAAGGCTACCAtttgttt                                                                                                                    >  1:311408/1‑143 (MQ=255)
       tAATTATTTTTTGGGTTGAATCCATTCCGGGGGGCAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAtttgtttg                                                                                                                   <  2:311408/143‑1 (MQ=255)
        aaTTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGc                                                                                                                  <  1:138228/143‑1 (MQ=255)
        aaTTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGc                                                                                                                  <  1:153764/143‑1 (MQ=255)
          ccacttcttGAGCGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTg                                                                                                                <  1:512477/136‑1 (MQ=255)
           tattTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCTCTGGTAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAACTGGCGTAACGTCCTACGCCTACCTTTGGTTGCggg                                                                                                               >  2:324509/1‑143 (MQ=255)
               ttttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATAGTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCCTACATTTGTTTGCGGGCCTg                                                                                                           >  1:67482/1‑143 (MQ=255)
                tttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTgg                                                                                                          >  1:61677/1‑143 (MQ=255)
                    aGTGAACTCCATACAGGACGCGAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGAttt                                                                                                      <  2:380522/143‑1 (MQ=255)
                      tGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGc                                                                                                    <  2:284137/143‑1 (MQ=255)
                      tGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATATCCTCAAAATGGCGTACCGTCCTATGCTTACATTCGTTTGCTGGACTGGATTTGc                                                                                                    >  1:45260/1‑143 (MQ=255)
                           tCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATAcc                                                                                               <  1:508856/143‑1 (MQ=255)
                            ccATACAGGGGGCAAATCAAAAAAAGTCTATGTTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGACGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCg                                                                                              <  1:336620/143‑1 (MQ=255)
                            ccATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGCTTGCTGCCCTGGACTCGCTTACCg                                                                                              >  2:96088/1‑143 (MQ=255)
                             cATACAGGGGGAAAATCAAAAAAATTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCgg                                                                                             <  1:133901/143‑1 (MQ=255)
                                aCAGGCGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCt                                                                                          <  2:522788/143‑1 (MQ=255)
                                     gggCAAATCAAAAAAAGTCTATGTTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCGTGTCCTCAAAAGGGCGTAGCCTCCTATGCTTACATTTGTTTCCTGGCCTGGATTTGCCAACCGGGCTGttct                                                                                     >  2:38570/1‑141 (MQ=255)
                                      ggCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGtgatt                                                                                    >  2:433890/1‑143 (MQ=255)
                                             cAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAGAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgctt                                                                             <  2:263770/143‑1 (MQ=255)
                                                  aaaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATCATGCGCAGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATGACGGGCTGGTATTATTGCGTGCCCt                                                                        >  1:645576/1‑143 (MQ=255)
                                                   caGTCTATGTTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAGTGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCTTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTg                                                                       <  1:20029/142‑1 (MQ=255)
                                                     gTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGGTATTATTGCTTGCAATGGc                                                                     >  2:425694/1‑143 (MQ=255)
                                                      tCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCt                                                                    >  2:478989/1‑143 (MQ=255)
                                                            tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                              <  1:479292/143‑1 (MQ=255)
                                                            tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                              >  1:589805/1‑143 (MQ=255)
                                                              tCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGACACGTCCTATGCCTACACTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCt                                                            >  2:623462/1‑143 (MQ=255)
                                                                       ccgcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGa                                                   >  1:368293/1‑143 (MQ=255)
                                                                         gcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAg                                                 <  2:625765/143‑1 (MQ=255)
                                                                            ccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATCGCGGGCTGTGAGTCTTGCTTGCACTGGCTTTTGTGCTGTTCTCCGAAGCCa                                              >  2:504567/1‑143 (MQ=255)
                                                                            ccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTCCGAAGccc                                              >  1:321945/1‑142 (MQ=255)
                                                                                gAAAGTCACTGATAATGCGCCGCGTTAATGTCCTCAAAATGGCGTAACGTACTATGCTTACATTTGTTCGCTGGCCTGGATTTGCATACCGGGCAGTTATTATTGCTTGCACTGTCCTTTGTGCTGTTTTACGAAGCCATCaa                                          >  2:242949/1‑143 (MQ=255)
                                                                                     tcccTGATAATGGGCCGCGTTCGTGTCCTCAACATGGCGTAACGTCCTATGCCTACATTTGTTGGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGtt                                     <  1:68876/140‑1 (MQ=255)
                                                                                       aCTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTc                                   <  2:635845/143‑1 (MQ=255)
                                                                                        cTGATAATTCGCCGCGTTCATGTCCTCAGAATGGCGGAACGGCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCAGACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTcc                                  <  2:106151/143‑1 (MQ=255)
                                                                                            tAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTCCTGGCCTGGATTTGCCTACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACTAAGCCATCAATGGTTTCCATGa                              >  2:63533/1‑143 (MQ=255)
                                                                                                   ccgcGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGcc                       <  2:61677/143‑1 (MQ=255)
                                                                                                           aTGTCCTCAAAAGGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGt               <  1:365681/143‑1 (MQ=255)
                                                                                                                   aaaaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCg       <  1:509353/143‑1 (MQ=255)
                                                                                                                        ggCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATc  <  1:261965/143‑1 (MQ=255)
                                                                                                                        ggCGCACCGTCCTACGCTTCCATTTCTTTGCTGGCCTGGATTTGCATACCGGGCTGGTATTATTGCTTGCACTGGCGTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATc  <  2:526133/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                              
TATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATC  >  NC_000913/3637509‑3637770

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: