Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,631,337	G→A	W112* (TGG→TGA)	pitA →	inorganic phosphate transporter PitA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,631,337	0	G	A	100.0%	198.0 / NA	53	W112* (TGG→TGA)	pitA	inorganic phosphate transporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/22);  total (31/22)

GGATCGTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGGTACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACC  >  NZ_CP009273/3631266‑3631418                                                                        |                                                                                  ggATCGTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCa                                                                 >  1:101087/1‑90 (MQ=255) ggATCGTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCa                                                                 >  1:880459/1‑90 (MQ=255) ggATCGTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCa                                                                 >  2:550600/1‑90 (MQ=255) ggATCGTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCa                                                                 >  1:389494/1‑90 (MQ=255)      gTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAg                                                            >  2:703948/1‑90 (MQ=255)      gTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAg                                                            >  1:88947/1‑90 (MQ=255)      gTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAg                                                            >  1:503324/1‑90 (MQ=255)      gTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAg                                                            >  1:314498/1‑90 (MQ=255)      gTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAg                                                            >  1:276098/1‑90 (MQ=255)        ctcATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGct                                                          <  1:446254/90‑1 (MQ=255)           aTGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCtctc                                                       >  1:168046/1‑90 (MQ=255)           aTGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCtctc                                                       >  2:825903/1‑90 (MQ=255)            tGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCa                                                      <  1:97118/90‑1 (MQ=255)                 ttGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGc                                                 >  2:581077/1‑90 (MQ=255)                 ttGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGc                                                 <  2:503324/90‑1 (MQ=255)                 ttGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGc                                                 >  1:750374/1‑90 (MQ=255)                 ttGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGc                                                 >  2:860593/1‑90 (MQ=255)                 ttGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGc                                                 >  1:431890/1‑90 (MQ=255)                  tGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCt                                                >  2:741209/1‑90 (MQ=255)                  tGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCt                                                >  2:663924/1‑90 (MQ=255)                  tGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCt                                                >  2:839475/1‑90 (MQ=255)                   gCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTg                                               <  1:581077/90‑1 (MQ=255)                   gCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTg                                               <  1:860593/90‑1 (MQ=255)                        ggTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGc                                                 <  2:766135/83‑1 (MQ=255)                        ggTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGc                                                 >  1:766135/1‑83 (MQ=255)                           gtTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGgcgc                                       >  1:753681/1‑90 (MQ=255)                           gtTCTCTATATTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGgcgc                                       <  1:863217/90‑1 (MQ=255)                             tctcTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGa                                     >  2:781655/1‑90 (MQ=255)                             tctcTATGTTGCTGGCGGCGATTATCTAGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGa                                     >  1:693490/1‑90 (MQ=255)                                  aTGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGatcatc                                >  1:369109/1‑90 (MQ=255)                                     ttGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTTGGTGCATGCAGGTCTCTTACGCTGGGTGGGGCGATCATTggg                             >  2:770274/1‑90 (MQ=18)                                     ttGCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCggg                             >  1:69040/1‑90 (MQ=255)                                       gCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           <  1:550600/90‑1 (MQ=255)                                       gCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           <  2:750374/90‑1 (MQ=255)                                       gCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           <  1:598080/90‑1 (MQ=255)                                       gCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           <  2:431890/90‑1 (MQ=255)                                       gCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           <  2:369109/90‑1 (MQ=255)                                       gCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           <  2:88947/90‑1 (MQ=255)                                       gCTGGCGGCGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           <  1:831304/90‑1 (MQ=255)                                          ggcggcGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGatcatc                                <  2:198985/82‑1 (MQ=255)                                          ggcggcGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGatcatc                                >  1:198985/1‑82 (MQ=255)                                          ggcggcGATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTgg                        >  1:306625/1‑90 (MQ=255)                                            cggcgATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGtt                      >  1:556470/1‑90 (MQ=255)                                              gcgATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTa                    <  1:781655/90‑1 (MQ=255)                                                gATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           <  1:562763/81‑1 (MQ=255)                                                gATTATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGAt                           >  2:562763/1‑81 (MQ=255)                                                   tATCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCaa               <  2:276098/90‑1 (MQ=255)                                                     tCTGGAACCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATg             <  1:703948/90‑1 (MQ=255)                                                           aCCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTtga       >  2:697514/1‑90 (MQ=255)                                                           aCCTGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTtga       <  2:880459/90‑1 (MQ=255)                                                              tGGGTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTtgatga    <  2:607257/90‑1 (MQ=255)                                                               gggTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGAc   >  2:106965/1‑90 (MQ=255)                                                                ggTACCTGATACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGAcc  <  1:202827/90‑1 (MQ=255)                                                                        |                                                                                  GGATCGTCTCATGGCCTTGCCATGGTGTTCTCTATGTTGCTGGCGGCGATTATCTGGAACCTGGGTACCTGGTACTTTGGTTTACCTGCATCCAGCTCTCATACGCTGATTGGCGCGATCATCGGGATTGGTTTAACCAATGCGTTGATGACC  >  NZ_CP009273/3631266‑3631418

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑