Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,901,765	G→A	G152G (GGC→GGT)	bglH ←	carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,901,765	0	G	A	100.0%	159.7 / NA	46	G152G (GGC→GGT)	bglH	carbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/21);  total (25/21)

AAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGA  >  NC_000913/3901630‑3901899                                                                                                                                        |                                                                                                                                       aaaCCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGAccc                                                                                                                                     <  1:510701/139‑1 (MQ=255)  aaCAAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCa                                                                                                                                    <  1:523491/139‑1 (MQ=255)       gCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGaa                                                                                                                               >  1:368486/1‑139 (MQ=255)                gTACTGCTGACCAACGTTACCATCCATCATCACAACGGGATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                      <  2:639170/139‑1 (MQ=255)                gTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                      <  2:208345/139‑1 (MQ=255)                 tACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGccc                                                                                                                     >  2:745226/1‑139 (MQ=255)                     gctgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                          <  2:458809/130‑1 (MQ=255)                     gctgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                          >  1:458809/1‑130 (MQ=255)                            aaCGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                          <  2:84323/139‑1 (MQ=255)                                         atcatcACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAgg                                                                                             <  1:343524/139‑1 (MQ=255)                                           catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                           >  2:681872/1‑139 (MQ=255)                                               acaACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTcc                                                                                       <  1:482400/139‑1 (MQ=255)                                                 aaCGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTcccc                                                                                     >  1:477056/1‑139 (MQ=255)                                                  gcGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTccccc                                                                                    >  2:66792/2‑139 (MQ=255)                                                      cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                      >  1:127337/1‑101 (MQ=255)                                                      cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                      <  2:127337/101‑1 (MQ=255)                                                      cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAcc                                                                                >  2:803228/1‑139 (MQ=255)                                                      cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAcc                                                                                >  1:270223/1‑139 (MQ=255)                                                        tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGa                                                                                                        >  1:711729/1‑113 (MQ=255)                                                        tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGa                                                                                                        <  2:711729/113‑1 (MQ=255)                                                         cAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGa                                                                             <  1:803228/139‑1 (MQ=255)                                                          aaCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                            <  1:730064/139‑1 (MQ=255)                                                          aaCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                            <  1:154509/139‑1 (MQ=255)                                                           aCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAAc                                                                           >  2:439311/1‑139 (MQ=255)                                                            cccGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACg                                                                          <  2:574358/139‑1 (MQ=255)                                                                tttGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTaa                                                                      >  2:423572/1‑139 (MQ=255)                                                                cttGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTaa                                                                      >  2:502931/2‑139 (MQ=255)                                                                    ccGGTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGc                                                                  >  2:435771/1‑139 (MQ=255)                                                                                        ttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAtt                                              >  2:134341/1‑139 (MQ=255)                                                                                          gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                      <  1:734732/65‑1 (MQ=255)                                                                                          gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                      >  2:734732/1‑65 (MQ=255)                                                                                                     aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACt                                                        >  2:56136/1‑116 (MQ=255)                                                                                                     aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACt                                                        <  1:56136/116‑1 (MQ=255)                                                                                                      aaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACaaa                                <  1:462292/139‑1 (MQ=255)                                                                                                       aTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAcc                                                                                >  1:616484/1‑90 (MQ=255)                                                                                                       aTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAcc                                                                                <  2:616484/90‑1 (MQ=255)                                                                                                                 ccGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAAt                     >  2:428425/1‑139 (MQ=255)                                                                                                                 ccGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAAt                     >  2:562718/1‑139 (MQ=255)                                                                                                                       tATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCa               <  1:681872/139‑1 (MQ=255)                                                                                                                             ttACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCa         <  1:134341/139‑1 (MQ=255)                                                                                                                               aCCAAAGCGACCCAGAGAACCTATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAg                                                                              <  2:58725/68‑1 (MQ=255)                                                                                                                               aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAg                                                                              >  1:58725/1‑68 (MQ=255)                                                                                                                                   aaGCGACCCAGAGAACCCATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaag   >  2:460978/1‑139 (MQ=255)                                                                                                                                   aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaag   >  2:224613/1‑139 (MQ=255)                                                                                                                                    aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaaga  >  1:421022/1‑139 (MQ=255)                                                                                                                                    aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaaga  >  2:535353/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       AAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGA  >  NC_000913/3901630‑3901899

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