Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	720,542	A→G	intergenic (‑82/‑41)	speF ← / → ybfK	ornithine decarboxylase, degradative/uncharacterized protein YbfK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	720,542	0	A	G	100.0%	178.9 / NA	48	intergenic (‑82/‑41)	speF/ybfK	ornithine decarboxylase, degradative/uncharacterized protein YbfK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/25);  total (23/25)

TGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGA  >  NC_000913/720409‑720674                                                                                                                                      |                                                                                                                                     tGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGc                                                                                                                                >  2:377408/1‑140 (MQ=255) tGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGc                                                                                                                                <  2:496677/140‑1 (MQ=255)            tCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGAtaaataaa                                                                                                                     >  1:66931/1‑140 (MQ=255)            tCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGAtaaataaa                                                                                                                     >  2:5202/1‑140 (MQ=255)             ccGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATaaataaat                                                                                                                    <  1:80064/140‑1 (MQ=255)             ccGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATaaataaat                                                                                                                    >  2:213166/1‑140 (MQ=255)                         cACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGc                                                                                                        <  2:69962/140‑1 (MQ=255)                                cGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGTAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTc                                                                                                 >  1:499598/1‑140 (MQ=255)                                cGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTc                                                                                                 >  1:64285/1‑140 (MQ=255)                                    aTTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGaa                                                                                             <  2:64285/140‑1 (MQ=255)                                      ttttAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATg                                                                                           <  1:214199/140‑1 (MQ=255)                                             tttGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGTAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCTGAAAc                                                                                    <  2:499598/140‑1 (MQ=255)                                             tttGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAAc                                                                                    <  1:377408/140‑1 (MQ=255)                                               tGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAAccc                                                                                  >  2:45018/1‑140 (MQ=255)                                                       ttCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGt                                                                          >  2:313332/1‑140 (MQ=255)                                                       ttCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGt                                                                          <  1:471648/140‑1 (MQ=255)                                                           tctcTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAgg                                                                      <  1:434256/140‑1 (MQ=255)                                                               tttATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCt                                                                                 <  1:245782/125‑1 (MQ=255)                                                               tttATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCt                                                                                 >  2:245782/1‑125 (MQ=255)                                                                         gTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATTATACa                                                        <  2:359102/140‑1 (MQ=255)                                                                         gTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACa                                                        <  2:40181/140‑1 (MQ=255)                                                                         gTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACa                                                        >  1:359555/1‑140 (MQ=255)                                                                             tAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTgg                                                                   <  2:174591/125‑1 (MQ=255)                                                                             tAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAggg                                                    >  1:140119/1‑140 (MQ=255)                                                                             tAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAggg                                                    >  2:201334/1‑140 (MQ=255)                                                                             tAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAggg                                                    >  2:499572/1‑140 (MQ=255)                                                                             tAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGGAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTgg                                                                   >  1:174591/1‑125 (MQ=255)                                                                              aGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGt                                                   <  2:140119/140‑1 (MQ=255)                                                                                  ggTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAAc                                               >  1:259451/1‑140 (MQ=255)                                                                                  ggTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAAc                                               <  1:213166/140‑1 (MQ=255)                                                                                   gTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACa                                              >  1:416172/1‑140 (MQ=255)                                                                                   gTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACa                                              >  1:414030/1‑140 (MQ=255)                                                                                     ccTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACaga                                            >  2:471711/1‑140 (MQ=255)                                                                                      cTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACagag                                           >  1:245665/1‑140 (MQ=255)                                                                                        cAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGaa                                         >  2:391152/1‑140 (MQ=255)                                                                                            ttCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAg                                     <  1:313332/140‑1 (MQ=255)                                                                                              ctcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCa                                   <  2:236105/140‑1 (MQ=255)                                                                                              ctcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCa                                   <  2:292168/140‑1 (MQ=255)                                                                                              ctcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCa                                   <  1:201334/140‑1 (MQ=255)                                                                                                cGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACt                                 >  1:308508/1‑140 (MQ=255)                                                                                                 gAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTg                                >  2:163098/1‑140 (MQ=255)                                                                                                 gAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTg                                <  2:245665/140‑1 (MQ=255)                                                                                                   aTGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                              <  2:53715/140‑1 (MQ=255)                                                                                                        aaaTCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATa                         >  1:263366/1‑140 (MQ=255)                                                                                                        aaaTCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATa                         <  2:66931/140‑1 (MQ=255)                                                                                                         aaTCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATAc                        <  1:5202/140‑1 (MQ=255)                                                                                                         aaTCGATCGTGCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATAc                        <  1:163098/140‑1 (MQ=255)                                                                                                                   gCCGGAAAAAAGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGc              <  2:180033/140‑1 (MQ=255)                                                                                                                             aGCAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACa    <  2:308508/140‑1 (MQ=255)                                                                                                                               cAAAGTTGCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACaga  <  2:416172/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     TGGTAAAACAGTCCGGGCAAGAATCACTAACCGCAATTTTTAATTTTGACATTTTTCATCTCTTTATTTTAGGTAATAGCAGGTCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGA  >  NC_000913/720409‑720674

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑