Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	986,671	A→G	R104R (CGT→CGC)	ompF ←	outer membrane porin 1a (Ia,b,F)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	986,671	0	A	G	100.0%	223.3 / NA	64	R104R (CGT→CGC)	ompF	outer membrane porin 1a (Ia,b,F)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/28);  total (36/28)

ATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATG  >  NC_000913/986539‑986806                                                                                                                                     |                                                                                                                                        aTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGtt                                                                                                                                   >  1:1205532/1‑139 (MQ=255)  tCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGttt                                                                                                                                  <  2:1205532/139‑1 (MQ=255)      tGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGtt                                                                                                                              >  1:57814/1‑139 (MQ=255)                    ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGt                                                                                                                >  1:492451/1‑139 (MQ=255)                    ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGt                                                                                                                >  1:404672/1‑139 (MQ=255)                         aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAc                                                                                                                            <  1:157441/122‑1 (MQ=255)                         aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAc                                                                                                                            >  2:157441/1‑122 (MQ=255)                         aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCAGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcg                                                                                                           >  2:258485/1‑139 (MQ=255)                                      atatCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              <  2:583089/139‑1 (MQ=255)                                           ggTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAccc                                                                                         <  2:404672/139‑1 (MQ=255)                                                aaCCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGaa                                                                                    >  2:532456/1‑139 (MQ=255)                                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTc                                                                                                               <  2:253858/105‑1 (MQ=255)                                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTc                                                                                                               >  1:253858/1‑105 (MQ=255)                                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCt                                                                                        >  2:783608/1‑128 (MQ=255)                                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCt                                                                                        <  1:783608/128‑1 (MQ=255)                                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                                             >  1:142130/1‑139 (MQ=255)                                                       tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                                             >  1:23036/1‑139 (MQ=255)                                                        gCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAACGTATTTAAGACCCGCGAATGCCATGCGCGTTTAGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATAtt                                                                            <  2:57814/139‑1 (MQ=255)                                                                 aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   >  2:343943/1‑139 (MQ=255)                                                                 aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   >  2:1115122/1‑139 (MQ=255)                                                                 aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   >  2:1259929/1‑139 (MQ=255)                                                                 aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   <  2:853903/139‑1 (MQ=255)                                                                     cacCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTg                                                                                                                                <  1:720997/74‑1 (MQ=255)                                                                     cacCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTg                                                                                                                                >  2:720997/1‑74 (MQ=255)                                                                      acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGcc                                                                                                         <  2:991801/96‑1 (MQ=255)                                                                      acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGcc                                                                                                         >  1:991801/1‑96 (MQ=255)                                                                      acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTg                                                                     >  2:956199/1‑132 (MQ=255)                                                                      acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTg                                                                     <  1:956199/132‑1 (MQ=255)                                                                         gTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                            >  1:230041/1‑90 (MQ=255)                                                                         gTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                            <  2:230041/90‑1 (MQ=255)                                                                          tAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  <  1:1135050/115‑1 (MQ=255)                                                                          tAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  >  2:1135050/1‑115 (MQ=255)                                                                          tAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCgg                                                          >  2:1285176/1‑139 (MQ=255)                                                                           aGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGt                                                         <  1:1157889/139‑1 (MQ=255)                                                                            gTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCtt                                                                                                       >  1:217391/1‑92 (MQ=255)                                                                            gTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCtt                                                                                                       <  2:217391/92‑1 (MQ=255)                                                                                  gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   <  1:184475/122‑1 (MQ=255)                                                                                  gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   >  2:184475/1‑122 (MQ=255)                                                                                       aaTCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAAtt                                             >  2:305108/1‑139 (MQ=255)                                                                                                aaCCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGt                                    >  1:1210866/1‑139 (MQ=255)                                                                                                   cAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTc                                 <  1:305108/139‑1 (MQ=255)                                                                                                      cGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   <  2:1014151/102‑1 (MQ=255)                                                                                                      cGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   >  1:1014151/1‑102 (MQ=255)                                                                                                          aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGAt                                                    <  2:903424/113‑1 (MQ=255)                                                                                                          aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGAt                                                    >  1:903424/1‑113 (MQ=255)                                                                                                               aTTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAAcc                     >  1:1174653/1‑139 (MQ=255)                                                                                                                   aaGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                 <  1:1259929/139‑1 (MQ=255)                                                                                                                    aGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGAc                <  1:532456/139‑1 (MQ=255)                                                                                                                    aGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGAc                <  2:857547/139‑1 (MQ=255)                                                                                                                         cgcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCa           >  1:1309043/1‑139 (MQ=255)                                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCg                                                                                                                 <  1:81136/36‑1 (MQ=255)                                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCg                                                                                                                 >  2:81136/1‑36 (MQ=255)                                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                                             <  1:488862/40‑1 (MQ=255)                                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                                             >  2:488862/1‑40 (MQ=255)                                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  >  2:630410/1‑67 (MQ=255)                                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  <  1:630410/67‑1 (MQ=255)                                                                                                                          gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCAt          >  1:698045/1‑139 (MQ=255)                                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              <  1:1441374/54‑1 (MQ=255)                                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              >  2:1441374/1‑54 (MQ=255)                                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAc                                                                                           >  1:1329852/1‑57 (MQ=255)                                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAc                                                                                           <  2:1329852/57‑1 (MQ=255)                                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATa         >  1:656404/1‑139 (MQ=255)                                                                                                                           cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTGGTTACCCTGGAAGTTATATTCCCAGTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAATCCAAGACGGGCATa         >  1:939900/1‑139 (MQ=255)                                                                                                                                  cAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg  <  2:698045/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        ATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATG  >  NC_000913/986539‑986806

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