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_000913	3,637,645	1	.	T	100.0%	129.3 / NA	43	L2L (CTA→TTA)	pitA	metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (18/25);  total (18/25)

TATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑TACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATAC  >  NC_000913/3637509‑3637775                                                                                                                                          |                                                                                                                                   tataAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑ttacgt                                                                                                                              >  1:583735/1‑138 (MQ=255)     aGATGATTATTTTTTGAGTGAAATCCATCCAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTCGCCTGCGCCGTTACAGTGAATTTGAAGTCGCCGCGTTCATGTACTCACAATGGCGTAACGTCCTATGCTTACAtttgt                                                                                                                           <  2:792009/143‑1 (MQ=255)         aaTTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGc                                                                                                                       >  2:74816/1‑143 (MQ=255)             attTTTTGAGTCAAATCCACACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGc                                                                                                                   <  2:547580/143‑1 (MQ=255)               gttctGAGTGAAAGACATACAGAGGGCACATCAAAAGAAGTCTAGAGGTCACTTTGCCTGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAAGGGCGTACCGTCCTATGCTTACATTTGTTTGCTGGCCt                                                                                                                 <  2:220174/139‑1 (MQ=255)                   tGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGAt                                                                                                             >  1:790500/1‑143 (MQ=255)                     aGTGAAACCCATACAGGCGTCAAATCAAAAAAAGTCTATATTCCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGAttt                                                                                                           <  2:49666/143‑1 (MQ=255)                      gTGAAATCCATACAGCGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTg                                                                                                          <  2:468140/143‑1 (MQ=255)                       tGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCCCTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGc                                                                                                         >  2:342793/1‑143 (MQ=255)                          aaTCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATa                                                                                                      >  1:786293/1‑143 (MQ=255)                           agccagccAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATAc                                                                                                     <  1:1092423/136‑1 (MQ=255)                               aTACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCggg                                                                                                 >  2:34912/1‑143 (MQ=255)                                      gggCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGatat                                                                                          >  2:202669/1‑143 (MQ=255)                                        gCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGATTGTATATAAAAATTCTCCGCGTTCTTTTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGttatta                                                                                        <  2:790500/143‑1 (MQ=255)                                                      gTCTACACTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAGCGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGc                                                                          <  2:906992/143‑1 (MQ=255)                                                        cTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCtt                                                                        >  1:318742/1‑143 (MQ=255)                                                         tataTTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCttt                                                                       <  1:728612/143‑1 (MQ=255)                                                           tctttCACTTCGCCCGCCACGCGTAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtg                                                                     <  1:1023286/141‑1 (MQ=255)                                                           tatTTCACTCTGCCCGCGCCGCGAAAGTCACTGACAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtg                                                                     <  2:883825/143‑1 (MQ=255)                                                             tttCACGTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCGATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTtgtg                                                                   <  2:764142/143‑1 (MQ=255)                                                              ttCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGc                                                                  <  1:68727/143‑1 (MQ=255)                                                                   tttGCCCGCGCCGCGAACGTCACTGATAATGCGCCGCGTTCATGTCCTCAAACTGGCGGAACGTCCTATGCTTACATTCGATTGCTGGCCGGGATTTGCATACCGGGCTGTTATTATTCGTTGCACTGGGTTTGGTGCTGTTc                                                             >  2:86768/1‑143 (MQ=255)                                                                      tctcgcgccgcgAATGTCACTGATAATGCGCCGCGGTCATGTCCTCAAAATGGCGTAACGTCCGAGGCTTACATTTGGTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTAc                                                          <  1:887084/140‑1 (MQ=255)                                                                         cgcgccgcgAAAGCCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGaa                                                       <  1:74816/143‑1 (MQ=255)                                                                                 gAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCaa                                               >  1:1062835/1‑143 (MQ=255)                                                                                      tCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGGAACGTCCTATGCTTACATTTGTTTGCTGGCCTAGATTCGCATACCGGGCTGTTATTATTGCCTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGtt                                          <  2:396077/143‑1 (MQ=255)                                                                                       cACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCGATGCTTACATTTGTTTGCTGGCCTGGATTTTCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGttt                                         <  2:81773/143‑1 (MQ=255)                                                                                        aCTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCGGTTCTACGAAGCCATCAATGGTTTc                                        >  2:932584/1‑143 (MQ=255)                                                                                         cTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTCTGTGCTGTTCTACGAAGCCGTCCATGGTTTcc                                       >  1:111517/1‑143 (MQ=255)                                                                                           gATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCAt                                     <  2:174285/143‑1 (MQ=255)                                                                                             ctaTGCGCCGTGTTCATGTCCTCAAAATGGCGTAACGTCCTCTGCTTACATTTGTTTGCTGGCCTGGATTGGCATACCGGGCTGTTATTATTGCTTGCACTTGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGa                                   <  2:903414/141‑1 (MQ=255)                                                                                                tGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGacac                                <  2:289645/143‑1 (MQ=255)                                                                                                   gccgcGTTCATGTCCTCAAAATTGCGTAACGTCCTATGCTTACATGTGTTTGCTGGCCTGCATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCCGTTCTACGAAGCCATCAATGGTTTCCATGACACAGc                             <  1:274946/143‑1 (MQ=255)                                                                                                          tCATGTCCTCAAAATGGCTTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGCTTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGcc                      <  1:725574/143‑1 (MQ=255)                                                                                                           cATGTCCTCGAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCg                     <  2:401448/143‑1 (MQ=255)                                                                                                              gTCCTCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCACACGCCGTgc                  >  2:216402/1‑142 (MQ=255)                                                                                                                  tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGCCACAGCCAACGCCGTGGCCAc              >  1:609491/1‑143 (MQ=255)                                                                                                                  tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCATGGATTTGCATACCGGGACGCTATGATTGCTTGCACTGGCTTGTGTTCTGATCTAACAAACCATCCATGGTTTGCATTCGCCAGCCAACCCCCTgaacaa              >  2:787907/1‑138 (MQ=255)                                                                                                                  tCAAAATGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCATGGATTTACATACGGGGCTGTTATTATTGCTTGCACTGGCTATTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGGCACAGCCAGCGGCGTGGCagc              >  2:149877/1‑141 (MQ=255)                                                                                                                    aaaaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCa                                      >  1:266845/1‑117 (MQ=255)                                                                                                                    aaaaTGGCGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCa                                      <  2:266845/117‑1 (MQ=255)                                                                                                                           cGTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACCGCCAACGCCGTGGCCACCGTTATCta     >  2:682154/1‑143 (MQ=255)                                                                                                                            gTAACGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCtat    >  1:901138/1‑143 (MQ=255)                                                                                                                              aaCGTCCTATGCTTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATAc  <  1:589386/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                   TATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGC‑TACATTTGTTTGCTGGCCTGGATTTGCATACCGGGCTGTTATTATTGCTTGCACTGGCTTTTGTGCTGTTCTACGAAGCCATCAATGGTTTCCATGACACAGCCAACGCCGTGGCAACCGTTATCTATAC  >  NC_000913/3637509‑3637775

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