Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,631,895	G→A	S232S (AGC→AGT)	yhiJ ←	DUF4049 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,631,895	0	G	A	100.0%	183.8 / NA	48	S232S (AGC→AGT)	yhiJ	DUF4049 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/18);  total (30/18)

TAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGGCTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCAACGGGTAAACAATTGGAG  >  NC_000913/3631754‑3632036                                                                                                                                              |                                                                                                                                              tAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGaa                                                                                                                                       >  2:311092/1‑150 (MQ=255)   aaTAAGGCTGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAAtt                                                                                                                                    >  2:184299/1‑151 (MQ=255)      aaGGCTGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAAtttt                                                                                                                                  <  1:401531/150‑1 (MQ=255)           tGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCtt                                                                                                                            >  2:192636/1‑151 (MQ=255)            gAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTc                                                                                                                           <  1:192635/151‑1 (MQ=255)             aaGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTc                                                                                                                           >  1:401565/1‑150 (MQ=255)             aaGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTc                                                                                                                           >  1:465477/1‑150 (MQ=255)                         aataatGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTccc                                                                                                              >  1:284327/1‑151 (MQ=255)                             ctGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAAt                                                                                                          >  2:478738/2‑151 (MQ=255)                             atGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAAt                                                                                                          >  2:478754/1‑151 (MQ=255)                                    ttctGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGt                                                                                                   >  2:59120/1‑151 (MQ=255)                                     tctGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGtt                                                                                                  >  2:135404/1‑151 (MQ=255)                                      ctGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTg                                                                                                 >  2:257679/1‑151 (MQ=255)                                           cAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTATCAATACGTTGAGAACTTCATCTTCTATATTAAGCAGCTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGCCCCGTTGCTaaa                                                                                            <  1:325715/151‑1 (MQ=255)                                           cAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTCTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGaa                                                                                             >  2:203820/1‑150 (MQ=255)                                               aTACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAgcg                                                                                        >  2:273393/1‑151 (MQ=255)                                                   tGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGa                                                                                     >  1:256527/1‑150 (MQ=255)                                                              ttAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCAc                                                                         >  2:333912/1‑151 (MQ=255)                                                                aaCATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCg                                                                       >  2:306888/1‑151 (MQ=255)                                                                aaCATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCg                                                                       <  2:104180/151‑1 (MQ=255)                                                                      aGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTg                                                                   >  2:226854/1‑149 (MQ=255)                                                                       gtgtTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGc                                                                >  2:593985/1‑151 (MQ=255)                                                                           tAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATc                                                            >  2:334139/1‑151 (MQ=255)                                                                            aaGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATc                                                            <  2:401569/150‑1 (MQ=255)                                                                            aaGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATc                                                            <  1:213338/150‑1 (MQ=255)                                                                                     aTTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCAcc                                                   >  2:543102/1‑150 (MQ=255)                                                                                     aTTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCa                                                  <  1:289341/151‑1 (MQ=255)                                                                                            ttCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAgg                                           >  2:544399/1‑151 (MQ=255)                                                                                             tctcCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCAc                                                    >  2:535370/1‑141 (MQ=255)                                                                                             tctcCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCAc                                                    <  1:535367/141‑1 (MQ=255)                                                                                             tctcCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGt                                          >  2:447645/1‑151 (MQ=255)                                                                                                cccTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTgg                                       >  1:91977/1‑151 (MQ=255)                                                                                                  cTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTc                                     >  1:521326/1‑151 (MQ=255)                                                                                                         aTACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACt                                                                    <  2:615027/113‑1 (MQ=255)                                                                                                         aTACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACt                                                                    >  1:615016/1‑113 (MQ=255)                                                                                                         aTACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACt                              >  1:107522/1‑151 (MQ=255)                                                                                                         aTACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACt                              <  1:369723/151‑1 (MQ=255)                                                                                                         aTACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACt                              >  1:325567/1‑151 (MQ=255)                                                                                                          tACGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACtt                             >  2:520916/1‑151 (MQ=255)                                                                                                           aCGTTGAGAACTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACtt                             <  1:579756/150‑1 (MQ=255)                                                                                                                   aaCTTCATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCa                    <  1:322057/151‑1 (MQ=255)                                                                                                                        cATCTTCTATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCAACggg               <  1:308069/151‑1 (MQ=255)                                                                                                                             tctATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCAACGGGt              <  2:116902/147‑1 (MQ=255)                                                                                                                             tctATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCAACGGGTaa            <  2:325569/149‑1 (MQ=255)                                                                                                                             tctATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCAACGGGTaa            <  2:24911/149‑1 (MQ=255)                                                                                                                             tctATATTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCAACGGGTaa            <  2:230736/149‑1 (MQ=255)                                                                                                                                atatTAAGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCAACGGGTAAACAAt       >  2:28205/1‑151 (MQ=255)                                                                                                                                     aaGCAACTGACTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATATGTTGGTCTGGTACTTGTGGGATCAACGGGTAAACAATTGgag  <  1:226854/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                              TAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCATCTTCTATATTAAGCAACTGGCTAAGGAAATTTTCTCCTTCAGCCATGATTCCCAAATGTCCCGTTGCGAAAAGCGCGAACAAAACTGCACCGACTGAGCCATCGGTATCACCAAAATAGGTTGGTCTGGTACTTGTGGGATCAACGGGTAAACAATTGGAG  >  NC_000913/3631754‑3632036

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