Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	1,712,421	C→A	intergenic (+263/‑348)	nth → / → dtpA	DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	1,712,421	0	C	A	100.0%	160.2 / NA	55	intergenic (+263/‑348)	nth/dtpA	DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/34);  total (21/34)

ATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAA  >  NC_000913_3_bme_pgi/1712276‑1712559                                                                                                                                                  |                                                                                                                                           aTAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACa                                                                                                                                        <  1:592451/150‑1 (MQ=255)  tAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTAAAAATCATTTAAATTTACACACGAAACAAATATTGACATaaaa                                                                                                                                       >  2:765186/1‑150 (MQ=255)     ccAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAAc                                                                                                                                    <  2:732454/150‑1 (MQ=255)      cAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACa                                                                                                                                   >  1:229903/1‑150 (MQ=255)      cAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACa                                                                                                                                   <  2:265702/150‑1 (MQ=255)        aaaaTGCAGTTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACAtt                                                                                                                                 <  2:548955/150‑1 (MQ=255)        aaaaTGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACAtt                                                                                                                                 <  2:115073/150‑1 (MQ=255)          aaTGCAGGTTAATTGTTTTTTTAATAGCGTAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTac                                                                                                                               <  1:272854/150‑1 (MQ=255)          aaTGCAGGTTAATTGTTTTTTTAATAGCCATATTTTCTATTCATTCGTGATCAAGAACACGAAGCTGTCGCATTGTAAAAACATTGCAAGAATAAGGTTAAAATCATTTAAATTTACACCCGCAACAAATAATGACATACAAAACATTac                                                                                                                               >  1:687622/1‑150 (MQ=255)          aaTGCAGGGTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTac                                                                                                                               <  1:823465/150‑1 (MQ=255)             ccagctaaaTTGTTTTTTAAATAGCGAAATTTTCTAATCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGAAAGCATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACAAACAAAACATTACACt                                                                                                                            >  1:567542/8‑150 (MQ=255)              cAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTg                                                                                                                           >  2:932778/1‑150 (MQ=255)                ggTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACTCAACAAAAATTGACATACAAAACATTACACtcga                                                                                                                         >  1:459145/1‑147 (MQ=255)                ggTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTgac                                                                                                                         >  1:528562/1‑148 (MQ=255)                   tAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTAGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTAt                                                                                                                      >  2:451236/1‑150 (MQ=255)                     aTTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTAttt                                                                                                                    <  1:823822/150‑1 (MQ=255)                     aTTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTAttt                                                                                                                    <  1:827650/150‑1 (MQ=255)                        gTTTTTTTAATAGCGAATTTTTCTATTCATTCGTCATCAAGATCACGAAGCTGTCGCAATGTAAAAACATTGCAAGAATCACGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTc                                                                                                                 <  2:320521/150‑1 (MQ=255)                         tttttttAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACAAACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCa                                                                                                                >  1:828212/1‑150 (MQ=255)                         tttttttAATAGCGAAATTATCTATTCATTCGTGATCAATATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACAAACAAAACATTACACTGGCTATTTTTCa                                                                                                                <  1:812423/150‑1 (MQ=255)                          ttttttAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAg                                                                                                               <  1:834760/150‑1 (MQ=255)                          ttttttAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAg                                                                                                               <  1:782202/150‑1 (MQ=255)                          ttttttAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAg                                                                                                               <  2:257843/150‑1 (MQ=255)                                aaTAGCGAAATTTTCTATTCATTCTTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATCTTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTg                                                                                                         <  1:709518/150‑1 (MQ=255)                                aaTAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTg                                                                                                         <  1:788631/150‑1 (MQ=255)                                   aGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGAc                                                                                                      <  2:103705/150‑1 (MQ=255)                                   aGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGAc                                                                                                      <  2:80733/150‑1 (MQ=255)                                    gCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACt                                                                                                     >  2:695432/1‑150 (MQ=255)                                            ttCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACt                                                                                             <  1:476855/150‑1 (MQ=255)                                              ttATTCATTCTTGATCAAGATCACGATGCTGTCGAAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAAAAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAg                                                                                           <  2:623554/149‑1 (MQ=255)                                              cTATTCATTCGTGATCAAGATCACGATGATGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAg                                                                                           <  2:495116/150‑1 (MQ=255)                                              atATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAg                                                                                           <  1:767659/149‑1 (MQ=255)                                                attcattcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTa                                                                                         <  2:787471/150‑1 (MQ=255)                                                       cGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGc                                                                                  >  1:238423/1‑150 (MQ=255)                                                        gTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCt                                                                                 >  1:914647/1‑150 (MQ=255)                                                          gATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATTTTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTaa                                                                               >  2:881131/1‑150 (MQ=255)                                                              aaGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAg                                                                           <  2:950498/150‑1 (MQ=255)                                                                gATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCa                                                                         >  1:316410/1‑150 (MQ=255)                                                                 aTCACGATGCTGTCGAAATGTAAAAACATTGCAAGATTCTGGTTAATATCATTTAAATTTAGAGAAGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTATATAGCAg                                                                        <  2:567542/150‑1 (MQ=255)                                                                      gATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACat                                                                   <  2:246310/150‑1 (MQ=255)                                                                                      aaaaaCATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACCCTGGATATTATTAACAAACTGTACTATATCACTAGTTAAAACGCTAAATCGCAGAACATATCGCCGTTAAGCCAt                                                   >  2:846989/1‑150 (MQ=255)                                                                                      aaaaaCATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAt                                                   >  2:246490/1‑150 (MQ=255)                                                                                       aaaaCATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATc                                                  >  2:148549/1‑150 (MQ=255)                                                                                               gCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAATCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCct                                          <  2:291315/150‑1 (MQ=255)                                                                                                cAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATTTTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCctc                                         >  2:136745/1‑150 (MQ=255)                                                                                                cAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCctc                                         <  2:528562/150‑1 (MQ=255)                                                                                                aaaGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACAATACACTGGCTATTTTTCTGAATCTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCctc                                         <  2:459145/149‑1 (MQ=255)                                                                                                           ttAAAATCATTTAAATTTACACACGAAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCATAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTcagtgcagtg                              <  1:925929/150‑1 (MQ=255)                                                                                                              aaaTCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGaaa                           <  2:130715/150‑1 (MQ=255)                                                                                                                          tttACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGtt               <  2:914647/150‑1 (MQ=255)                                                                                                                                cacGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCTTTACGttt         >  2:770431/1‑150 (MQ=255)                                                                                                                                    caacaaATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCTCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg     <  2:740230/150‑1 (MQ=255)                                                                                                                                     gacaaATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGACCATATCGCCGTTATGCCATCACCCCCCTCTCATTGCAGTGAAAAAATCTGCCGTTCCGTTTTTTGa    >  1:383901/2‑150 (MQ=255)                                                                                                                                       caaATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaa  >  2:572877/1‑150 (MQ=255)                                                                                                                                       aaaaTATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaa  <  2:910634/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                           ATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAA  >  NC_000913_3_bme_pgi/1712276‑1712559

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑