Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,424,386	G→A	*229* (TAG→TAA)	ulaF →	L‑ribulose 5‑phosphate 4‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,424,386	0	G	A	100.0%	215.8 / NA	66	*229* (TAG→TAA)	ulaF	L‑ribulose 5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/32);  total (34/32)

GATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAG  >  NC_000913/4424254‑4424519                                                                                                                                     |                                                                                                                                      gATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACAcg                                                                                                                                 <  2:147019/139‑1 (MQ=255)       agaagTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGc                                                                                                                           >  1:1381514/1‑139 (MQ=255)         aagTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGTTGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGTTGCgg                                                                                                                         >  1:1029689/1‑139 (MQ=255)           gTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGaa                                                                                                                       >  1:440622/1‑139 (MQ=255)               cGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAAttt                                                                                                                   >  1:75475/1‑139 (MQ=255)                  aaaTGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATttcct                                                                                                                <  2:75475/139‑1 (MQ=255)                           ggATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                       <  2:1598242/139‑1 (MQ=255)                                 cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATa                                                                                                   <  1:845373/137‑1 (MQ=255)                                 cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATa                                                                                                   >  2:845373/1‑137 (MQ=255)                                 cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 <  1:902922/139‑1 (MQ=255)                                 cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 >  1:382058/1‑139 (MQ=255)                                 cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 <  2:1347057/139‑1 (MQ=255)                                              cccACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCg                                                                                    >  1:1326797/1‑139 (MQ=255)                                                  caACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 <  1:1466704/122‑1 (MQ=255)                                                  caACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 >  2:1466704/1‑122 (MQ=255)                                                  caACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGc                                                                                >  1:800052/1‑139 (MQ=255)                                                         aTCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGttttt                                                                         >  2:286937/1‑139 (MQ=255)                                                          tCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTc                                                                        >  2:121105/1‑139 (MQ=255)                                                             caTCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaa                                                                     >  2:1564412/1‑139 (MQ=255)                                                                     gCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGcc                                                             <  1:286937/139‑1 (MQ=255)                                                                        tCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 >  1:153500/1‑100 (MQ=255)                                                                        tCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 <  2:153500/100‑1 (MQ=255)                                                                         ccTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTACGCAGCCGTTTTTCAAAAAAAAGCCCCCt                                                         <  1:1106427/139‑1 (MQ=255)                                                                                aaTAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAgg                                                  >  1:1512320/1‑139 (MQ=255)                                                                                aaTAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAgg                                                  >  1:379711/1‑139 (MQ=255)                                                                                  tAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAgggg                                                >  2:1102016/1‑139 (MQ=255)                                                                                   aaaCACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAggggg                                               <  2:379711/139‑1 (MQ=255)                                                                                         ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACAcg                                                                                                                                 >  1:1564745/1‑51 (MQ=255)                                                                                         ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACAcg                                                                                                                                 <  2:1564745/51‑1 (MQ=255)                                                                                         ttCATCCGTAACCACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACTCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTC‑AAAAAAAGCCCCCTGCGAAGGGGgcaaagca                                        <  2:1626446/139‑1 (MQ=255)                                                                                          tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca                                        >  1:1094444/1‑139 (MQ=255)                                                                                          tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca                                        <  1:1564412/139‑1 (MQ=255)                                                                                          tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca                                        <  1:1368349/139‑1 (MQ=255)                                                                                           cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATttcc                                                                                                                 >  1:747565/1‑65 (MQ=255)                                                                                           cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATttcc                                                                                                                 <  2:747565/65‑1 (MQ=255)                                                                                             tGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATttcc                                                                                                                 >  2:578938/1‑63 (MQ=255)                                                                                             tGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATttcc                                                                                                                 <  1:578938/63‑1 (MQ=255)                                                                                                        ggTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 <  1:1335693/68‑1 (MQ=255)                                                                                                        ggTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 >  2:1335693/1‑68 (MQ=255)                                                                                                          tCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAAtt                                                                                       <  2:805992/76‑1 (MQ=255)                                                                                                          tCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAAtt                                                                                       >  1:805992/1‑76 (MQ=255)                                                                                                          tCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGttt                        >  1:81652/1‑139 (MQ=255)                                                                                                           ccTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTc                       <  2:800052/139‑1 (MQ=255)                                                                                                            cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                      <  2:1136752/139‑1 (MQ=255)                                                                                                            cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                      <  2:1326797/139‑1 (MQ=255)                                                                                                            cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                      <  2:1512320/139‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATttcct                                                                                                                >  1:1555086/1‑41 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATttcct                                                                                                                <  2:1555086/41‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 >  2:1314036/1‑56 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                                 <  1:1314036/56‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTgg                                                                                              <  2:1413599/59‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTgg                                                                                              >  1:1413599/1‑59 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGc                               <  1:339556/122‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGc                               >  2:339556/1‑122 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                      >  2:1312638/1‑131 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                      <  1:1312638/131‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat              >  2:1319532/1‑139 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat              >  1:1134966/1‑139 (MQ=255)                                                                                                                        aCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCt          <  1:739814/139‑1 (MQ=255)                                                                                                                          gggCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGTTAACTGGTCTAATTCCACAGCCGTTTTTCAAAAAAAAGCACCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTgg        >  1:1420933/1‑139 (MQ=255)                                                                                                                           ggCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt       >  1:209428/1‑139 (MQ=255)                                                                                                                              aGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTgtt    <  1:1102016/139‑3 (MQ=255)                                                                                                                               gAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAACTTATGGCAAGGTTTCGTTGGTTATACCTGGTGTTa   >  1:980021/1‑139 (MQ=255)                                                                                                                                aaGTAAAACACGCGCTGCGGAAATTTCCTTCCTCgg                                                                                                         <  1:1342596/36‑1 (MQ=255)                                                                                                                                aaGTAAAACACGCGCTGCGGAAATTTCCTTCCTCgg                                                                                                         >  2:1342596/1‑36 (MQ=255)                                                                                                                                aaGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAg  <  1:1134366/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      GATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAG  >  NC_000913/4424254‑4424519

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