Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,088,867	G→A	Q331* (CAA→TAA)	pgaB ←	poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine N‑deacetylase and beta‑1,6 glycoside hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,088,867	0	G	A	100.0%	175.2 / NA	50	Q331* (CAA→TAA)	pgaB	poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine N‑deacetylase and beta‑1,6 glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/28);  total (22/28)

GGTAGCAAACGATTTGGAAACCAGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTGCTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAATGAGGGATTATTGGCGATTAATACCCGCGGAA  >  NC_000913/1088731‑1088998                                                                                                                                         |                                                                                                                                    ggTAGCAAACGATTTGGAAACCAGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTg                                                                                                                                  >  2:886678/1‑140 (MQ=255)         aCGATTTGGAAACCAGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGtttt                                                                                                                          <  1:287354/140‑1 (MQ=255)          cGATTTGGAAACCAGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTc                                                                                                                         <  1:518170/140‑1 (MQ=255)            aTTTGGAAACCAGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGt                                                                                                                       >  1:1164516/1‑140 (MQ=255)                   aaCCAGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAAc                                                                                                                >  2:788297/1‑140 (MQ=255)                      cAGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTa                                                                                                             <  1:819298/140‑1 (MQ=255)                       aGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTaa                                                                                                            >  2:66512/1‑140 (MQ=255)                       aGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTaa                                                                                                            >  1:184485/1‑140 (MQ=255)                        gACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAAt                                                                                                           <  1:836085/140‑1 (MQ=255)                         aCCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATc                                                                                                          <  1:66512/140‑1 (MQ=255)                           ctctTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCaa                                                                                                        <  1:570962/140‑1 (MQ=255)                            tctTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAg                                                                                                       <  1:788297/140‑1 (MQ=255)                             ctTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGa                                                                                                      <  2:1104565/140‑1 (MQ=255)                               ttGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATc                                                                                                    <  1:37269/140‑1 (MQ=255)                                    cAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGatat                                                                                               <  1:680518/140‑1 (MQ=255)                                      gCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGc                                                                                             >  1:918539/1‑140 (MQ=255)                                       cccATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCa                                                                                            >  2:966215/1‑140 (MQ=255)                                       cccATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCa                                                                                            >  2:849200/1‑140 (MQ=255)                                       cccATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCa                                                                                            <  1:966215/140‑1 (MQ=255)                                          aTCGCCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAg                                                                                                       <  2:754698/126‑1 (MQ=255)                                          aTCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAg                                                                                                       >  1:754698/1‑126 (MQ=255)                                                 tCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTtg                                                                                  >  1:559727/1‑140 (MQ=255)                                                 tCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTtg                                                                                  >  2:249889/1‑140 (MQ=255)                                                   gggATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTtgtg                                                                                <  2:184485/140‑1 (MQ=255)                                                   gggATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTtgtg                                                                                <  2:1164516/140‑1 (MQ=255)                                                   gggATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTtgtg                                                                                <  2:304528/140‑1 (MQ=255)                                                   gggATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTtgtg                                                                                <  2:231515/140‑1 (MQ=255)                                                      aTCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTgtggtg                                                                             <  1:945395/140‑1 (MQ=255)                                                      aTCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTgtggtg                                                                             >  2:1150473/1‑140 (MQ=255)                                                       tCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGa                                                                            >  2:239769/1‑140 (MQ=255)                                                        cAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGAt                                                                           >  2:959395/1‑140 (MQ=255)                                                        cAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGAt                                                                           <  1:959395/140‑1 (MQ=255)                                                         aGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGAtt                                                                          >  2:994639/1‑140 (MQ=255)                                                                gCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTg                                                                   <  1:249889/140‑1 (MQ=255)                                                                 ccTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGt                                                                  <  1:849200/140‑1 (MQ=255)                                                                 ccTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGt                                                                  <  2:96862/140‑1 (MQ=255)                                                                  cTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTa                                                                 <  1:1150473/140‑1 (MQ=255)                                                                    gCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACg                                                               >  2:970878/1‑140 (MQ=255)                                                                                 tGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTa                                                                 <  2:1022478/125‑1 (MQ=255)                                                                                 tGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTa                                                                 >  1:1022478/1‑125 (MQ=255)                                                                                      tttGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCa                                             <  1:886678/140‑1 (MQ=255)                                                                                             atCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCtt                                      <  2:1079256/140‑1 (MQ=255)                                                                                                  tCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAAt                                 >  2:260998/1‑140 (MQ=255)                                                                                                  tCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAAt                                 >  2:790283/1‑140 (MQ=255)                                                                                                  tCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAAt                                 >  2:1091096/1‑140 (MQ=255)                                                                                                         cTGAATTAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAATGAGGGat                          >  1:1124847/1‑140 (MQ=255)                                                                                                               tAGCACATCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAATGAGGGATTATTgg                    <  2:559727/140‑1 (MQ=255)                                                                                                                      tCAATATTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAATGAGGGATTATTGGCGATTaa             <  1:239769/140‑1 (MQ=255)                                                                                                                         atatTGCGATCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAATGAGGGATTATTGGCGATTAATAc          <  2:1124847/140‑1 (MQ=255)                                                                                                                                 aTCCATTTACTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAATGAGGGATTATTGGCGATTAATACCCGCggaa  >  1:715594/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                    GGTAGCAAACGATTTGGAAACCAGACCTCTTTGACCAGCCCATCACCATCGGGATCAGCAAATGCCTGCAAATACACGGTTGATATTTGCATATCTTTCACCCGCTGAATTAGCACATCAATATTGCGATCCATTTGCTGGAGGTTTTCGTCATAAACGTAATCAAGATCGATATGCATTATCCGTTGTGGTGATTTTTCCTGTACGGTAATAATTTGCTGGGCAAACTCTTTTAATGAGGGATTATTGGCGATTAATACCCGCGGAA  >  NC_000913/1088731‑1088998

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑