Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,425,980	G→A	pseudogene (351/453 nt)	ydaY →	pseudogene, Rac prophage,Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,425,980	0	G	A	100.0%	105.7 / NA	35	pseudogene (351/453 nt)	ydaY	pseudogene, Rac prophage,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/17);  total (18/17)

GGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATT  >  NC_000913/1425857‑1426115                                                                                                                            |                                                                                                                                        ggCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggtggt                                                                                                                          >  2:221733/1‑139 (MQ=255)      aaCGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCg                                                                                                                     <  2:392808/139‑1 (MQ=255)        cGAGGGATCGCGCTGCTTAAGCAGAATCCTGAGTTGTACTAGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGGAGCAGAAGCAGGATAAGATAAAAAAGCCATAACTGCCCAGCGGGGGTTATCGGa                                                                                                                   >  1:525144/1‑139 (MQ=255)        cGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGAACTGGCAGATGCTTGCAGAAGCAGCCGGCATCGCCCAGGGTAAAGAGCAGTAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGa                                                                                                                   >  2:511927/1‑139 (MQ=255)            ggATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTc                                                                                                               <  1:221733/139‑1 (MQ=255)                       ctTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAg                                                                                                    <  2:430471/139‑1 (MQ=255)                                    cAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCa                                                                                       <  2:679230/139‑1 (MQ=255)                                          tGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCa                                                                                                                                  <  2:783506/90‑1 (MQ=255)                                          tGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCa                                                                                                                                  >  1:783506/1‑90 (MQ=255)                                                 gCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGaa                                                                                             <  2:748971/120‑1 (MQ=255)                                                 gCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGaa                                                                                             >  1:748971/1‑120 (MQ=255)                                                 gCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCgag                                                                          >  2:546916/1‑139 (MQ=255)                                                   aGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGaa                                                                        >  1:463123/1‑139 (MQ=255)                                                           gCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAg                                                                >  2:202823/1‑139 (MQ=255)                                                                 gcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                                    <  1:31676/65‑1 (MQ=255)                                                                 gcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                                    >  2:31676/1‑65 (MQ=255)                                                                    gccagcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCAcc                                                       >  1:170465/1‑139 (MQ=255)                                                                     ccagcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCg                                                      >  2:774387/1‑139 (MQ=255)                                                                                  gggTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAAcc                                           <  1:580116/137‑1 (MQ=255)                                                                                  gggTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAAcc                                           >  2:580116/1‑137 (MQ=255)                                                                                            cagcagAAGCAGGATACGATAAAAAATCCATAACTGCTCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGAccc                               <  1:511927/139‑1 (MQ=255)                                                                                            cagcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGAccc                               <  1:767900/139‑1 (MQ=255)                                                                                               cagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGt                            >  1:470111/1‑139 (MQ=255)                                                                                                      aGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAg                     <  2:463123/139‑1 (MQ=255)                                                                                                      aGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAg                     <  1:202823/139‑1 (MQ=255)                                                                                                         aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGa                                        <  1:690712/117‑1 (MQ=255)                                                                                                         aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGa                                        >  2:690712/1‑117 (MQ=255)                                                                                                         aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg                  >  2:557779/1‑139 (MQ=255)                                                                                                                  aaaaGCCATAACTGCCCAGCGGTGGTTAGCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATGGACCCGGTGCTTAAGGAGTTGTTGTAc         >  1:548414/1‑139 (MQ=255)                                                                                                                       ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCAc                                                        >  1:778204/1‑87 (MQ=255)                                                                                                                       ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCAc                                                        <  2:778204/87‑1 (MQ=255)                                                                                                                       ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAAcc                                           >  2:580812/1‑100 (MQ=255)                                                                                                                       ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAAcc                                           <  1:580812/100‑1 (MQ=255)                                                                                                                       ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTa    <  2:491252/139‑1 (MQ=255)                                                                                                                         aTAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTAtt  <  1:189620/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                        GGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATT  >  NC_000913/1425857‑1426115

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