Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,655,360	(C)5→4	coding (314/327 nt)	ynfA ←	UPF0060 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,655,356	0	C	.	100.0%	141.9 / NA	35	coding (318/327 nt)	ynfA	UPF0060 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (21/14);  total (21/14)

GAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATAA  >  NC_000913/1655218‑1655488                                                                                                                                           |                                                                                                                                     gAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTgcgc                                                                                                                                      >  2:1100991/1‑139 (MQ=255)  aaaaaaGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcc                                                                                                                                     <  2:652144/139‑1 (MQ=255)   aaaaaGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccc                                                                                                                                    <  1:876014/139‑1 (MQ=255)     aaaGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                                  <  2:89508/139‑2 (MQ=255)     aaaGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGcccca                                                                                                                                  <  1:1106001/139‑2 (MQ=255)        gCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCGccccagcc                                                                                                                               <  2:156444/139‑5 (MQ=255)             aGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                                         <  2:154244/139‑1 (MQ=255)             aGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                                         >  1:357498/1‑139 (MQ=255)                      aaaaaaaGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCa                                                                                                                <  2:910808/139‑1 (MQ=255)                           aaGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATg                                                                                                           >  1:10989/1‑139 (MQ=255)                                        gcgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCaa                                                                                              >  2:875133/1‑139 (MQ=255)                                        gcgcTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCaa                                                                                              >  1:81382/1‑139 (MQ=255)                                                  tGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCAcc                                                                                    <  1:352284/139‑1 (MQ=255)                                                        aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                              >  1:80396/1‑139 (MQ=255)                                                        aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                              >  1:145802/1‑139 (MQ=255)                                                        aCAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCa                                                                              >  1:158613/1‑139 (MQ=255)                                                                    aTAGAATACCTTTTGTCAGCGGACACTTCCTCTTATCTTATTGATAAAATGGAGTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGCCCATTCATAAAGAGt                                                                  <  2:37902/139‑1 (MQ=255)                                                                        aaTACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGt                                                              >  2:1227293/1‑139 (MQ=255)                                                                          tACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGttt                                                            >  2:361010/1‑139 (MQ=255)                                                                              ttttGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACg                                                        <  1:1227293/139‑1 (MQ=255)                                                                                  gTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCAt                                                    >  1:807520/1‑139 (MQ=255)                                                                                       cTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCcaca                                               >  1:597104/1‑139 (MQ=255)                                                                                       cTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCcaca                                               >  1:768910/1‑139 (MQ=255)                                                                                          acacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCa                                                          <  1:1001014/125‑1 (MQ=255)                                                                                          acacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCa                                                          >  2:1001014/1‑125 (MQ=255)                                                                                          acacTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcg                                            <  1:1100991/139‑1 (MQ=255)                                                                                            acTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAAcgcg                                          <  1:302111/139‑1 (MQ=255)                                                                                              ttCCTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCa                                        <  2:81382/139‑1 (MQ=255)                                                                                                ccTCTTATCTTATTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGc                                      >  1:669403/1‑139 (MQ=255)                                                                                                     tatcttatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACAt                                 <  2:80396/139‑1 (MQ=255)                                                                                                      atcttatTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATc                                >  2:727027/1‑139 (MQ=255)                                                                                                           atTGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGc                           <  2:145802/139‑1 (MQ=255)                                                                                                            ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                                         <  1:159235/44‑1 (MQ=255)                                                                                                            ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGcaa                                                                                                                         >  2:159235/1‑44 (MQ=255)                                                                                                            ttGATAAAATGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGcc                          >  1:954931/1‑139 (MQ=255)                                                                                                                     tGGATTTATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACa                 >  2:735709/1‑139 (MQ=255)                                                                                                                         tttATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATTAACGCCGTGCAGACATaaa             <  1:875133/139‑1 (MQ=255)                                                                                                                          ttATATTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAAc            >  2:626852/1‑139 (MQ=255)                                                                                                                           tATGTTCTTACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACg           <  1:361010/139‑1 (MQ=255)                                                                                                                                 cttACGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGccacca     >  2:301744/3‑139 (MQ=255)                                                                                                                                    aCGTGCG‑CCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATaa  >  1:313375/1‑139 (MQ=255)                                                                                                                                           |                                                                                                                                     GAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTTTGTCAGCTGACACTTCCTCTTATCTTATTGATAAAATGGATTTATGTTCCTACGTGCGCCCCCAGCCCGCAACAATGATCAACATGCCGCAAAGCGCAATCAACGCACCCGTCCAGTCATAAAGAGTCAGTTTCACGCCATCCACAACGCGCAGCCACATCAACGCCGTGCAGACATAAACGCCACCATAA  >  NC_000913/1655218‑1655488

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: ‑