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%	143.4 / NA	46	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 (25/21);  total (25/21)

GGAATATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATT  >  NC_000913/1425848‑1426115                                                                                                                                     |                                                                                                                                        ggAATATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                                    >  1:659585/1‑138 (MQ=255) ggAATATACGCCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                                    <  2:659585/138‑1 (MQ=255)   aaTATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAA‑‑AGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCg                                                                                                                               <  1:157936/139‑1 (MQ=255)     tataCGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCg                                                                                                                               <  2:768563/139‑1 (MQ=255)            cAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTa                                                                                                                        >  2:706408/1‑139 (MQ=255)            cAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTa                                                                                                                        >  2:624456/1‑139 (MQ=255)            cAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGGAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTa                                                                                                                        >  2:768107/1‑139 (MQ=255)               aaCGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCg                                                                                                                               >  2:150795/1‑129 (MQ=255)               aaCGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCg                                                                                                                               <  1:150795/129‑1 (MQ=255)                         cgcgCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTc                                                                                                               >  1:482930/1‑135 (MQ=255)                         cgcgCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTc                                                                                                               <  2:482930/135‑1 (MQ=255)                             ctgctTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGgaga                                                                                                       <  1:706408/139‑1 (MQ=255)                                          aTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAg                                                                                          >  1:58107/1‑139 (MQ=255)                                          aTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTAGCGGAGTTCGGGGGAGAAAGGGGGGAAAAg                                                                                          >  1:850287/1‑139 (MQ=255)                                           tCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAgg                                                                                         <  1:624456/139‑1 (MQ=255)                                                               tCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGaaaaa                                                                     >  1:308369/1‑139 (MQ=255)                                                               tCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGaaaaa                                                                     >  2:215370/1‑139 (MQ=255)                                                               tCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGaaaaa                                                                     >  2:758573/1‑139 (MQ=255)                                                                    gCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAg                                                                >  1:901833/1‑139 (MQ=255)                                                                          gcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGtt                                                             <  2:349686/136‑1 (MQ=255)                                                                          gcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGtt                                                             >  1:349686/1‑136 (MQ=255)                                                                          gcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGcc                                                          >  2:167449/1‑139 (MQ=255)                                                                             gccagcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACt                                                                   >  2:266781/1‑127 (MQ=255)                                                                             gccagcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACt                                                                   <  1:266781/127‑1 (MQ=255)                                                                                   aTCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATAcc                                                 >  2:708944/1‑139 (MQ=255)                                                                                    tCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCg                                                <  2:901833/139‑1 (MQ=255)                                                                                       cccAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggt                                                                                                                             >  1:882813/1‑59 (MQ=255)                                                                                       cccAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggt                                                                                                                             <  2:882813/59‑1 (MQ=255)                                                                                       cccAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGaa                                             >  2:122042/1‑139 (MQ=255)                                                                                        ccAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAAc                                            <  1:54231/139‑1 (MQ=255)                                                                                        ccAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAAc                                            <  2:718626/139‑1 (MQ=255)                                                                                                      agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCg                                                                                                                               <  2:609361/42‑1 (MQ=255)                                                                                                      agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCg                                                                                                                               >  1:609361/1‑42 (MQ=255)                                                                                                       gcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCgg                             >  2:510305/1‑139 (MQ=255)                                                                                                        cagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGCAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGt                            >  2:364539/1‑139 (MQ=255)                                                                                                          gAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGGAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGc                          >  1:722372/1‑139 (MQ=255)                                                                                                             gCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCgggg                                                                                                           <  2:989387/55‑1 (MQ=255)                                                                                                             gCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCgggg                                                                                                           >  1:989387/1‑55 (MQ=255)                                                                                                               aGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAg                     <  2:722372/139‑1 (MQ=255)                                                                                                                          aaaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTa          <  2:308369/139‑1 (MQ=255)                                                                                                                           aaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg                  >  2:83916/1‑130 (MQ=255)                                                                                                                           aaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg                  <  1:83916/130‑1 (MQ=255)                                                                                                                           aaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTAc         >  1:744697/1‑139 (MQ=255)                                                                                                                             aaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGc       <  2:434823/139‑1 (MQ=255)                                                                                                                               gCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCt     <  1:758573/139‑1 (MQ=255)                                                                                                                                  aTAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTAtt  <  1:510305/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        GGAATATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATT  >  NC_000913/1425848‑1426115

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