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%	190.6 / NA	55	G152G (GGC→GGT)	bglH	carbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/27);  total (28/27)

CCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGA  >  NC_000913/3901637‑3901896                                                                                                                                 |                                                                                                                                    ccGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAAc                                                                                                                           >  1:358129/1‑139 (MQ=255)           aCTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCa                                                                                                                 <  1:512718/139‑1 (MQ=255)            ctgctgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCAt                                                                                                                >  2:227157/1‑139 (MQ=255)               ctgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAt                                                                                                             >  2:1043760/1‑139 (MQ=255)                  aCCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTa                                                                                                          >  1:816439/1‑139 (MQ=255)                  aCCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTa                                                                                                          >  2:588559/1‑139 (MQ=255)                  aCCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTa                                                                                                          >  1:604135/1‑139 (MQ=255)                           aCCATCCATCATGACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCa                                                                                                 >  1:631576/1‑139 (MQ=255)                               tccatcATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATgag                                                                                             <  1:44489/139‑1 (MQ=255)                                   tcatcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGc                                                                                         >  2:942801/1‑139 (MQ=255)                                    catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                        <  2:824613/139‑1 (MQ=255)                                    catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                        <  2:293834/139‑1 (MQ=255)                                    catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                        <  2:677548/139‑1 (MQ=255)                                    catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                        <  2:753088/139‑1 (MQ=255)                                         caACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTccc                                                                                   <  1:761334/139‑1 (MQ=255)                                          aaCGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                   <  1:595448/106‑1 (MQ=255)                                          aaCGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                   >  2:595448/1‑106 (MQ=255)                                                aTCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCa                                                                            >  1:436404/1‑139 (MQ=255)                                                aTCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCa                                                                            >  1:122522/1‑139 (MQ=255)                                                 tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAg                                                                           <  2:693877/139‑1 (MQ=255)                                                  cAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGa                                                                          <  2:358129/139‑1 (MQ=255)                                                   aaCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                         <  2:604135/139‑1 (MQ=255)                                                   aaCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                         <  1:188413/139‑1 (MQ=255)                                                            gCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAg                                                                <  2:694673/139‑1 (MQ=255)                                                             ccGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGc                                                               >  2:527367/1‑139 (MQ=255)                                                                ttttCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGt                                                            <  2:816439/139‑1 (MQ=255)                                                                                gttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGTCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTa                                                    <  1:364171/131‑1 (MQ=255)                                                                                gttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGTCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTa                                                    >  2:364171/1‑131 (MQ=255)                                                                                 ttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAGCTAAAGCCAAtt                                           >  2:356519/1‑139 (MQ=255)                                                                                   gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTc                                         >  2:595863/1‑139 (MQ=255)                                                                                          cTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTc                                  >  1:846275/1‑139 (MQ=255)                                                                                            gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                       >  2:649110/1‑52 (MQ=255)                                                                                            gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                       <  1:649110/52‑1 (MQ=255)                                                                                            gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                   <  1:1070919/56‑1 (MQ=255)                                                                                            gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                   >  2:1070919/1‑56 (MQ=255)                                                                                            gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                   <  2:655464/56‑1 (MQ=255)                                                                                            gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                   >  1:655464/1‑56 (MQ=255)                                                                                              aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACaa                              >  1:716772/1‑139 (MQ=255)                                                                                                    aaCCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTAc                        <  2:436404/139‑1 (MQ=255)                                                                                                       cAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGa                     >  2:723043/1‑139 (MQ=255)                                                                                                         gCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaaa                   >  2:532876/1‑139 (MQ=255)                                                                                                          ccGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                       >  1:718764/1‑54 (MQ=255)                                                                                                          ccGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                       <  2:718764/54‑1 (MQ=255)                                                                                                          ccGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAAt                  >  2:549525/1‑139 (MQ=255)                                                                                                           cGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATc                 >  1:695541/1‑139 (MQ=255)                                                                                                                 aTTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCag           <  2:695541/139‑1 (MQ=255)                                                                                                                   tCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCagag         <  1:588559/139‑1 (MQ=255)                                                                                                                        aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg    <  1:723043/139‑1 (MQ=255)                                                                                                                        aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg    <  1:549525/139‑1 (MQ=255)                                                                                                                         ccAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATgg   <  2:716772/139‑1 (MQ=255)                                                                                                                          cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGga  >  2:339428/1‑139 (MQ=255)                                                                                                                            aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                         <  1:111904/66‑1 (MQ=255)                                                                                                                            aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                         >  2:111904/1‑66 (MQ=255)                                                                                                                             aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAg                                                                <  1:983710/74‑1 (MQ=255)                                                                                                                             aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAg                                                                >  2:983710/1‑74 (MQ=255)                                                                                                                                 |                                                                                                                                    CCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGA  >  NC_000913/3901637‑3901896

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