Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,750,092	G→A	*287* (TAG→TAA)  E3K (GAG→AAG)	sgbU → sgbE →	putative L‑xylulose 5‑phosphate 3‑epimerase L‑ribulose‑5‑phosphate 4‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,750,092	0	G	A	100.0%	189.0 / NA	58	*287* (TAG→TAA)  E3K (GAG→AAG)	sgbU sgbE	putative L‑xylulose 5‑phosphate 3‑epimerase L‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/28);  total (30/28)

TGCATAAACTGAACTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTT  >  NC_000913/3749959‑3750215                                                                                                                                      |                                                                                                                            tGCATAAACTGAACTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCaa                                                                                                                        >  1:1643294/1‑139 (MQ=255)    aTAAACTGAACTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTg                                                                                                                     >  2:480342/1‑139 (MQ=255)     tAAACTGAACTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGa                                                                                                                    <  1:836639/139‑1 (MQ=255)                  cgcgGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCt                                                                                                       >  2:675861/1‑139 (MQ=255)                      gTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTggcg                                                                                                   >  2:1288067/1‑139 (MQ=255)                        tCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTggcggc                                                                                                 <  1:675861/139‑1 (MQ=255)                           ttCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGaa                                                                                              >  1:960793/1‑139 (MQ=255)                                    gagaTGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCt                                                                                     >  1:977906/1‑139 (MQ=255)                                         gtgGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCg                                                                                >  1:68493/1‑139 (MQ=255)                                           ggACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCaa                                                                                                                        <  2:447957/97‑1 (MQ=255)                                           ggACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCaa                                                                                                                        >  1:447957/1‑97 (MQ=255)                                            gACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTc                                                                             >  2:578396/1‑139 (MQ=255)                                               cGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCAcc                                                                          <  2:1420133/139‑1 (MQ=255)                                                 aaaaaGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTg                                                                                                         >  1:1286635/1‑106 (MQ=255)                                                 aaaaaGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTg                                                                                                         <  2:1286635/106‑1 (MQ=255)                                                 aaaaaGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCAt                                                                        <  2:68493/139‑1 (MQ=255)                                                           agagCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGt                                                              <  2:1643294/139‑1 (MQ=255)                                                                cGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCAcc                                                         <  2:960793/139‑1 (MQ=255)                                                                  gTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCt                                                        <  1:1232945/138‑1 (MQ=255)                                                                  gTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCt                                                        >  2:1232945/1‑138 (MQ=255)                                                                       ggAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCa                                                  >  2:27887/1‑139 (MQ=255)                                                                             tattCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTaaag                                                                                                                           <  2:1305293/60‑1 (MQ=255)                                                                             tattCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCt                                                                                     >  1:887848/1‑98 (MQ=255)                                                                             tattCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCt                                                                                     <  2:887848/98‑1 (MQ=255)                                                                                  aGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGcgcg                                       <  2:829128/139‑1 (MQ=255)                                                                                  aGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCCCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGcgcg                                       <  1:480342/139‑1 (MQ=255)                                                                                    gcgcGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGt                                                              <  1:716998/114‑1 (MQ=255)                                                                                    gcgcGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGt                                                              >  2:716998/1‑114 (MQ=255)                                                                                     cgcgGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCa                                                                            >  1:214888/1‑99 (MQ=255)                                                                                     cgcgGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCa                                                                            <  2:214888/99‑1 (MQ=255)                                                                                       cggcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGAc                                                                                                            <  2:1674650/65‑1 (MQ=255)                                                                                       cggcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGAc                                                                                                            >  1:1674650/1‑65 (MQ=255)                                                                                         gcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTcc                                                                                  <  2:1066860/89‑1 (MQ=255)                                                                                         gcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTcc                                                                                  >  1:1066860/1‑89 (MQ=255)                                                                                         gcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTc                                                                             >  2:907224/1‑94 (MQ=255)                                                                                         gcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTc                                                                             <  1:907224/94‑1 (MQ=255)                                                                                         gcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACg                                <  1:1288067/139‑1 (MQ=255)                                                                                            ttGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGaaa                             <  1:578396/139‑1 (MQ=255)                                                                                                aTTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGTGGCAATGTCAGCGCGGTAGACGAGAcgcg                         <  2:977906/139‑1 (MQ=255)                                                                                                    aaGCGCGTATGCAGGAGGCTGGATTTATATGTAAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCaa                     >  2:1004946/1‑139 (MQ=255)                                                                                                       cgcgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAgcgc                                        <  2:1419728/117‑1 (MQ=255)                                                                                                       cgcgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAgcgc                                        >  1:1419728/1‑117 (MQ=255)                                                                                                       cgcgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGc                       <  2:85240/134‑1 (MQ=255)                                                                                                       cgcgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGc                       >  1:85240/1‑134 (MQ=255)                                                                                                         cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAgcgc                                        <  2:1482327/115‑1 (MQ=255)                                                                                                         cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAgcgc                                        >  1:1482327/1‑115 (MQ=255)                                                                                                         cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAg                                           <  1:102568/112‑1 (MQ=255)                                                                                                         cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAg                                           >  2:102568/1‑112 (MQ=255)                                                                                                         cgTATGCAGGAGGCTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggat                >  2:1614568/1‑139 (MQ=255)                                                                                                                      cTGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCt   >  1:1390315/1‑139 (MQ=255)                                                                                                                       tGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGgcgc                                                                                      <  1:156500/55‑1 (MQ=255)                                                                                                                       tGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGgcgc                                                                                      >  2:156500/1‑55 (MQ=255)                                                                                                                       tGGATTTATATGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCtt  >  1:1282711/1‑139 (MQ=255)                                                                                                                              atatGTTAAAGCAACTGAAAGCGGACGTGCTGGCGGTGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGcgcg                                       <  2:137962/95‑1 (MQ=255)                                                                                                                              atatGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGcgcg                                       >  1:137962/1‑95 (MQ=255)                                                                                                                               tatGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTgg                                                                    >  1:762793/1‑65 (MQ=255)                                                                                                                               tatGTTAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTgg                                                                    <  2:762793/65‑1 (MQ=255)                                                                                                                                   ttAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAgcgc                                        <  1:445617/89‑1 (MQ=255)                                                                                                                                   ttAAAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAgcgc                                        >  2:445617/1‑89 (MQ=255)                                                                                                                                      |                                                                                                                            TGCATAAACTGAACTATCGCGGTTCTTTCCTGATTGAGATGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTT  >  NC_000913/3749959‑3750215

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