Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,637,646 +T coding (5/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.T96.2% 146.8 / NA 53L2L (CTA→TTA) pitAmetal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base T (28/22);  minor base G (0/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCCGCGCGAT  >  NC_000913/3637504‑3637784
                                                                                                                                              |                                                                                                                                            
aaTTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGACGTAACGTCCTATGcg                                                                                                                                              >  2:309398/1‑142 (MQ=255)
   tagataAGAGAATTATTTTTGGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCttac                                                                                                                                          <  1:427258/140‑4 (MQ=255)
     tataAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCttacat                                                                                                                                        >  1:911903/1‑138 (MQ=255)
       taAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAttt                                                                                                                                      <  1:1044362/143‑1 (MQ=255)
           aTAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTACGCTTACCtttgttt                                                                                                                                  >  2:370653/1‑143 (MQ=255)
             aaTTATTTTTTTAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGc                                                                                                                                <  1:328819/143‑1 (MQ=255)
             aaTTATTTTTTGAGTGCACTCCATGCACGGGGCAGAACAAAAAAAGTCTATATTTCGCTTTGCCCGCGCCGCGAAAGTCACTTATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGc                                                                                                                                <  2:225867/143‑1 (MQ=255)
                 attTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCACGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCCAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGc                                                                                                                            >  1:174102/1‑143 (MQ=255)
                  ttgtttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGGCGCGAAAGTGACTGAGAGTGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGcc                                                                                                                           >  1:482363/4‑143 (MQ=255)
                   tttttGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCt                                                                                                                          <  1:501746/143‑1 (MQ=255)
                         ggTGAAATCCCTGCAGGGGGCGAGTCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGAttt                                                                                                                    <  2:376814/142‑1 (MQ=255)
                         aGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGAttt                                                                                                                    <  1:806354/143‑1 (MQ=255)
                           tGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTgg                                                                                                                  >  1:417471/1‑142 (MQ=255)
                            gAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCGCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTGTGCTGGCCTGGCTTTGCa                                                                                                                 >  2:1150413/1‑143 (MQ=255)
                             aaaTCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCAt                                                                                                                >  1:384139/1‑143 (MQ=255)
                                   aTACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTACTGGCCTGGATTTGCATACCggg                                                                                                          >  2:29730/1‑143 (MQ=255)
                                     gcAGGGGGTCAATCAAAAAAAGTCTCTATTTCACTTCGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTAATTTTGTTTGCTTGTCTGGATTTGCATACCGGGCt                                                                                                        <  2:140990/142‑1 (MQ=255)
                                      cAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTAATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTg                                                                                                       >  1:968111/1‑143 (MQ=255)
                                       aGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGt                                                                                                      >  1:1052584/1‑143 (MQ=255)
                                         ggggCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTGCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACAGTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGtta                                                                                                    >  1:21336/1‑143 (MQ=255)
                                          gggCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAATATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGttat                                                                                                   >  1:897194/1‑143 (MQ=255)
                                                  cAAAAAAAGTCTATATTTCACTTTGTCCGCGCCGCGAAAGTCACTGCTAATGCGCAGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGCTGTTAttgctt                                                                                           >  1:1082107/1‑143 (MQ=255)
                                                   aaaaaaaTTCTATATTTCACGTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTAttgcttg                                                                                          <  2:365148/143‑1 (MQ=255)
                                                   aaaaaaaGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATAACGGGCTGTTATTAttgcttg                                                                                          >  2:963895/1‑143 (MQ=255)
                                                         aGTCTATATTTCACTTTGCCCGCGCCGCGACAGTCACGGATAACGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTAGGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTgg                                                                                    <  2:516864/143‑1 (MQ=255)
                                                          gTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGCCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGc                                                                                   <  1:539170/143‑1 (MQ=255)
                                                            cTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtt                                                                                 >  1:771639/1‑143 (MQ=255)
                                                             tataTTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGCCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTGATTATTGCTTGCACTGGtttt                                                                                >  2:159920/1‑143 (MQ=255)
                                                                 tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                                            >  2:895822/1‑143 (MQ=255)
                                                                 tttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTCACATTTGTTTGCTGGCCTGGGTTTGCATACCGGGCTGTTATTATTGCGTGCACTGGATTTtgtg                                                                            >  2:380302/1‑143 (MQ=255)
                                                                    cACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTg                                                                         >  1:340628/1‑143 (MQ=255)
                                                                      cTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATTGCGGTCTGCTATTATTTCTTGCACTGGCTTTTGTGCTGtt                                                                       >  1:307456/1‑143 (MQ=255)
                                                                              gcgccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTGTGCTGGCCTGGATTTGCATACCGGGCGGGTATTATTGCTTGCACGGGCATATGTGCTCTTCTACGAAg                                                               >  2:475972/1‑143 (MQ=255)
                                                                                 ccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCa                                                            <  1:619547/143‑1 (MQ=255)
                                                                                 ccgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCa                                                            <  2:844940/143‑1 (MQ=255)
                                                                                  cgcgAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACCTTTGTTTGCTGGCCTGGATTTGCATAACGGGCTGCTATTATTGCTTGCACTGGCTTTTGTACTGAACAAAGAAGCCAt                                                           >  1:470602/1‑143 (MQ=255)
                                                                                      aaaGGCACGGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAAt                                                       <  1:351505/143‑1 (MQ=255)
                                                                                       aaGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATGCCGGGCTGTTATTATTGCTTGCACTGGCtttt                                                                               <  2:408166/118‑1 (MQ=255)
                                                                                       aaGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtttt                                                                               >  1:408166/1‑118 (MQ=255)
                                                                                          tCACTGAAACTCGGCCGCGTGCAGGTCGTCAAAATGGCGTAACGGCCTATGCGTACATTTGTTTGCTGGCCTGGAGTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGtt                                                   <  2:952703/143‑1 (MQ=255)
                                                                                            aCTGATAATGCGCCGCGTTCATTTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTc                                                 <  1:630190/143‑1 (MQ=255)
                                                                                              tGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCa                                               <  1:691628/143‑1 (MQ=255)
                                                                                                   aTGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGaca                                          <  1:466317/143‑1 (MQ=255)
                                                                                                              tCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGcc                               <  2:897194/143‑1 (MQ=255)
                                                                                                                aTGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATTGCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGt                             <  2:551816/143‑1 (MQ=255)
                                                                                                                aTGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTCCTGGCCTGGCTTTGCATACCGGACTGTTATTATTGCTTGCACTGGCTTATGTGCTGTTCTACGAAGCCAGCAGTTGTTTCCATGCCACAGCCACCGcccg                             >  2:866368/1‑141 (MQ=255)
                                                                                                                   tCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCTGGCTCATATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTgcg                          >  2:816973/1‑141 (MQ=255)
                                                                                                                   tCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCgcgcc                          >  2:166073/1‑139 (MQ=255)
                                                                                                                      tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGTTTTCCATGACACAGCCACCGCCGTGGCAAc                       >  2:339574/1‑143 (MQ=255)
                                                                                                                       cAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAAcc                      <  1:269644/143‑1 (MQ=255)
                                                                                                                          aaTGGCGTGATGTCCTATTCGTCCATTTGTTTGCTGGCCTGGTTTTTCATACCGGGTTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGtt                   <  2:125020/143‑1 (MQ=255)
                                                                                                                            tGGCGTATCGTCCTATGCTTACATTTGTTTGCTGGCTTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACTAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTAt                 <  2:454378/143‑1 (MQ=255)
                                                                                                                                gTAACGGCCTATGCTTACATTTGTTTGCTGGCCTGGACTTGCATACCGGGCTGTTATTATTGCTAGCGCTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCtat             <  2:930638/143‑1 (MQ=255)
                                                                                                                                 tAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCtata            >  1:719530/1‑143 (MQ=255)
                                                                                                                                         tATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCACTGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTAGCTATACCcgcgcg    >  2:93465/1‑143 (MQ=255)
                                                                                                                                         tATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGGCGTGGCAACCGTTATCTATAcgcgcggg    >  1:401350/1‑141 (MQ=255)
                                                                                                                                          tgccTACATTTGTTTGCTGGGCTGGATTTGCATACCGGTCTGTTATTATTGCTTGCACTGGCTTTGGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCCGCGCGAt  <  2:568980/140‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
AATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATACCCGCGCGAT  >  NC_000913/3637504‑3637784

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: