Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,545,565	G→A	*228* (TAG→TAA)	gntX →	DNA catabolic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,545,565	0	G	A	100.0%	188.4 / NA	59	*228* (TAG→TAA)	gntX	DNA catabolic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/30);  total (29/30)

CCCGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCC  >  NC_000913/3545431‑3545699                                                                                                                                       |                                                                                                                                       cccGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGcc                                                                                                                                    <  2:7355/139‑1 (MQ=255)     tGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCga                                                                                                                                >  2:595972/1‑139 (MQ=255)      gCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgat                                                                                                                               >  2:798960/1‑139 (MQ=255)         aGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                            <  1:595972/139‑1 (MQ=255)         aGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                            <  2:71175/139‑1 (MQ=255)             cGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  2:1043451/139‑1 (MQ=255)                      gTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                               <  2:690837/139‑1 (MQ=255)                         aTTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             <  1:1017469/122‑1 (MQ=255)                         aTTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             >  2:1017469/1‑122 (MQ=255)                         aTTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCa                                                                                                            >  2:1085645/1‑139 (MQ=255)                                    tgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGCCGTATtata                                                                                                                >  1:313109/1‑124 (MQ=255)                                    tgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGCCGTATtata                                                                                                                <  2:313109/124‑1 (MQ=255)                                             caccGGAAGTACCGTCGCGGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                            <  2:810068/103‑1 (MQ=255)                                             caccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                            >  1:810068/1‑103 (MQ=255)                                                cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     <  1:446370/139‑1 (MQ=255)                                                  gAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   >  1:637697/1‑139 (MQ=255)                                                  gAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   <  1:231280/139‑1 (MQ=255)                                                  gAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCa                                                                                   <  1:702103/139‑1 (MQ=255)                                                    aGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCAtt                                                                                 >  2:566189/1‑139 (MQ=255)                                                      tACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                            >  2:913001/1‑94 (MQ=255)                                                      tACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatgat                                                                                                                            <  1:913001/94‑1 (MQ=255)                                                       aCCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        >  2:175589/1‑97 (MQ=255)                                                       aCCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCAAACCTTGTAAAGCCTCGATGATGGGc                                                                                                                        <  1:175589/97‑1 (MQ=255)                                                       aCCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGGGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACt                                                                              >  2:736027/1‑139 (MQ=255)                                                         cGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             <  2:308881/90‑1 (MQ=255)                                                         cGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             >  1:308881/1‑90 (MQ=255)                                                          gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                >  2:331635/1‑118 (MQ=255)                                                          gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                <  1:331635/118‑1 (MQ=255)                                                          gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATg                                                                           <  1:1085645/139‑1 (MQ=255)                                                           tCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTaaag                                                                                                                                      >  1:1062243/1‑79 (MQ=255)                                                           tCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTaaag                                                                                                                                      <  2:1062243/79‑1 (MQ=255)                                                                        gcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             >  1:334628/1‑75 (MQ=255)                                                                        gcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             <  2:334628/75‑1 (MQ=255)                                                                               cTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                         <  1:280654/72‑1 (MQ=255)                                                                               cTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                         >  2:280654/1‑72 (MQ=255)                                                                                      gCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGt                                                                     >  1:384426/1‑117 (MQ=255)                                                                                      gCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGt                                                                     <  2:384426/117‑1 (MQ=255)                                                                                          tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTCAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACtt                                           >  1:686192/1‑139 (MQ=255)                                                                                               cggcggCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCa                                      <  2:637697/139‑1 (MQ=255)                                                                                                    gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCaa                                                                                                           >  1:930448/1‑65 (MQ=255)                                                                                                    gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCaa                                                                                                           <  2:930448/65‑1 (MQ=255)                                                                                                      tGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             >  1:528581/1‑45 (MQ=255)                                                                                                      tGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                                             <  2:528581/45‑1 (MQ=255)                                                                                                       gTCCAGGTCTGGTGCCTTTGTCGAACTTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                          <  2:643448/79‑1 (MQ=255)                                                                                                       gTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTAt                                                                                          >  1:643448/1‑79 (MQ=255)                                                                                                        tCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                         <  2:141447/47‑1 (MQ=255)                                                                                                        tCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                                         >  1:141447/1‑47 (MQ=255)                                                                                                         ccAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGc                            <  1:798960/139‑1 (MQ=255)                                                                                                                  gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                               <  2:487909/47‑1 (MQ=255)                                                                                                                  gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                                               >  1:487909/1‑47 (MQ=255)                                                                                                                  gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                <  1:733119/62‑1 (MQ=255)                                                                                                                  gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTc                                                                                                >  2:733119/1‑62 (MQ=255)                                                                                                                  gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGga                   >  2:25628/1‑139 (MQ=255)                                                                                                                       tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaaa                                                                                                     >  1:35940/1‑52 (MQ=255)                                                                                                                       tttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTaaaa                                                                                                     <  2:35940/52‑1 (MQ=255)                                                                                                                        ttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     <  2:52487/67‑1 (MQ=255)                                                                                                                        ttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGc                                                                                     >  1:52487/1‑67 (MQ=255)                                                                                                                                ccTTGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACaaa     >  2:303919/1‑139 (MQ=255)                                                                                                                                   tGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAAt    <  1:407296/137‑1 (MQ=255)                                                                                                                                   tGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAAt    >  2:407296/1‑137 (MQ=255)                                                                                                                                   tGTAAAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATcc  >  2:83383/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       CCCGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAACCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGAAGAAGGGACACAAATCC  >  NC_000913/3545431‑3545699

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