Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,045,010	A→C	E271A (GAA→GCA)	yihF →	DUF945 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,045,010	0	A	C	100.0%	165.5 / NA	45	E271A (GAA→GCA)	yihF	DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/19);  total (26/19)

GTTGGTGAAAGCGAGTCGAAATTTAATCAGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTG  >  NC_000913/4044868‑4045153                                                                                                                                               |                                                                                                                                                gTTGGTGAAAGCGAGTCGAAATTTAATCAGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATa                                                                                                                                         <  1:794570/151‑1 (MQ=255)   tGGTGAAAGCGAGTCGAAATTTAATCAGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATa                                                                                                                                         >  1:318461/1‑149 (MQ=255)        aaaGCGAGTCGAAATTTAATCAGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAAc                                                                                                                                  >  1:919442/1‑151 (MQ=255)                   aaaTTTAATCAGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAAt                                                                                                                        >  1:635086/1‑150 (MQ=255)                       ttAATCAGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGc                                                                                                                    >  2:893655/1‑150 (MQ=255)                         aaTCAGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAAc                                                                                                                 >  2:215819/1‑151 (MQ=255)                             aGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCg                                                                                                              <  1:682252/150‑1 (MQ=255)                                     aTTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAgtggtg                                                                                                     <  2:725374/151‑1 (MQ=255)                                        atcaatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAg                                                                                                   <  2:266066/150‑1 (MQ=255)                                         tcaatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTg                                                                                                 >  2:645008/1‑151 (MQ=255)                                           aatcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGt                                                                                               >  2:904566/1‑151 (MQ=255)                                            atcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATAGCAATGTCAATCTCGCCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTgg                                                                                                <  2:384697/149‑1 (MQ=255)                                             tcaCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTc                                                                                              >  1:912272/1‑150 (MQ=255)                                               aCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCatt                                                                                           <  1:144166/151‑1 (MQ=255)                                               aCGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCatt                                                                                           <  1:18783/151‑1 (MQ=255)                                                cGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCattg                                                                                          <  1:833632/151‑1 (MQ=255)                                                  gggAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCattgat                                                                                        <  1:596222/151‑1 (MQ=255)                                                          gTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAAt                                                                                >  2:130132/1‑151 (MQ=255)                                                               ccAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTAtt                                                                           <  1:91434/151‑1 (MQ=255)                                                                  aaTCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGAt                                                                        >  1:358728/1‑151 (MQ=255)                                                                  aaTCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGAt                                                                        >  2:675776/1‑151 (MQ=255)                                                                        tGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTc                                                                   >  1:915417/1‑150 (MQ=255)                                                                         gCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTca                                                                  <  2:696651/150‑1 (MQ=255)                                                                         gCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTc                                                                   >  1:172864/1‑149 (MQ=255)                                                                            ttCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTcagca                                                               <  2:572752/150‑1 (MQ=255)                                                                              cGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTg                                                            >  2:215744/1‑151 (MQ=255)                                                                                  gCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGATAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGcc                                                        >  1:88910/1‑151 (MQ=255)                                                                                         ccAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTat                                                 >  1:536347/1‑151 (MQ=255)                                                                                         ccAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTa                                                  >  2:465676/1‑150 (MQ=255)                                                                                          cAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTatc                                                <  2:707198/151‑1 (MQ=255)                                                                                                ggATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATcc                                            <  2:172730/149‑1 (MQ=255)                                                                                                 gATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGt                                         >  1:379466/1‑151 (MQ=255)                                                                                                         tAAAGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACAtt                                 >  2:217501/1‑151 (MQ=255)                                                                                                          aaaGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTg                                >  2:415207/1‑151 (MQ=255)                                                                                                           aaGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGc                               <  1:906428/151‑1 (MQ=255)                                                                                                            aGACAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCg                              <  2:876418/151‑1 (MQ=255)                                                                                                              aCAAAGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGAt                            <  2:456252/151‑1 (MQ=255)                                                                                                                  aGATTATATCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGc                          >  2:261363/1‑149 (MQ=255)                                                                                                                         atCAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAagcag                  <  1:776710/150‑1 (MQ=255)                                                                                                                           cAATGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCtt               >  2:888382/1‑151 (MQ=255)                                                                                                                             aTGTCAATCTCACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTtgc             <  2:797916/151‑1 (MQ=255)                                                                                                                                   aTCTCACCTACGCACTTGATAAGATAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGCCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATCTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGca         >  2:101103/1‑149 (MQ=255)                                                                                                                                      tcACCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACAccc    >  2:185595/1‑151 (MQ=255)                                                                                                                                        aCCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTg  >  1:408068/1‑151 (MQ=255)                                                                                                                                        aCCTACGCACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTg  >  1:191220/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                GTTGGTGAAAGCGAGTCGAAATTTAATCAGATTAACATTATCAATCACGGGGAAGACGTTGCCCAAATCGATGCTTTCGTTGCAAAAACCAGGCTGGATCGCGTTAAAGACAAAGATTATATCAATGTCAATCTCACCTACGAACTTGATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTG  >  NC_000913/4044868‑4045153

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