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%	179.9 / NA	57	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 (29/28);  total (29/28)

ATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTAC  >  NC_000913/1425853‑1426108                                                                                                                                |                                                                                                                                 ataCGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCgg                                                                                                                       <  1:741843/139‑1 (MQ=255)     ggCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggtggt                                                                                                                   >  2:450403/1‑139 (MQ=255)          aaCGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCg                                                                                                              >  1:406291/1‑139 (MQ=255)               gggATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGtt                                                                                                         <  1:450403/139‑1 (MQ=255)                        ctgctTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGTAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGgaga                                                                                                <  1:888282/139‑1 (MQ=255)                        ctgctTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGgaga                                                                                                <  2:37323/139‑1 (MQ=255)                                    aaTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGgaga                                                                                                <  2:145890/127‑1 (MQ=255)                                    aaTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGgaga                                                                                                >  1:145890/1‑127 (MQ=255)                                     aTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                             >  2:661662/1‑97 (MQ=255)                                     aTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                             >  2:606800/1‑97 (MQ=255)                                     aTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                             <  1:661662/97‑1 (MQ=255)                                     aTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                             <  1:606800/97‑1 (MQ=255)                                         agagTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCac                                                                               >  2:1084642/1‑138 (MQ=255)                                                  cTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGc                                                                      >  2:1169361/1‑139 (MQ=255)                                                  cTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGAAACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGc                                                                      <  2:635573/139‑1 (MQ=255)                                                     gCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGc                                                                      <  1:1472843/136‑1 (MQ=255)                                                     gCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGc                                                                      >  2:1472843/1‑136 (MQ=255)                                                       aGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGaa                                                                 >  1:296751/1‑139 (MQ=255)                                                          tCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATg                                                                           <  1:386496/126‑1 (MQ=255)                                                          tCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATg                                                                           >  2:386496/1‑126 (MQ=255)                                                          tCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGaaaaa                                                              >  1:459596/1‑139 (MQ=255)                                                                     gcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGcc                                                   >  1:624293/1‑139 (MQ=255)                                                                        gccagcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGTGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCAcc                                                >  1:955332/1‑139 (MQ=255)                                                                        gccagcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCAcc                                                >  1:613834/1‑139 (MQ=255)                                                                        gccagcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCAcc                                                >  2:451854/1‑139 (MQ=255)                                                                             cATCGCCCAGGGAAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGCGTTCGGGGGTGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATAc                                           >  2:193484/1‑139 (MQ=255)                                                                                 gCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGa                                       <  1:175791/139‑1 (MQ=255)                                                                                    cAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAAcc                                    <  2:514981/139‑1 (MQ=255)                                                                                            agagCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAg                            >  1:1310028/1‑139 (MQ=255)                                                                                                cagcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGAccc                        <  1:1169361/139‑1 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggt                                                                                                                      <  2:1179856/44‑1 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggt                                                                                                                      >  1:1179856/1‑44 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCg                                                                                                              <  2:1103797/52‑1 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCg                                                                                                              >  1:1103797/1‑52 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGa                                                                                                            >  1:1441362/1‑54 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGa                                                                                                            <  2:1441362/54‑1 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCATGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCg                       >  2:986562/1‑139 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGc                                                                      <  2:1373665/92‑1 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGc                                                                      >  1:1373665/1‑92 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCgg                                        >  1:479247/1‑122 (MQ=255)                                                                                                 agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCgg                                        <  2:479247/122‑1 (MQ=255)                                                                                                        gCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACt                                                            >  2:1438107/1‑95 (MQ=255)                                                                                                        gCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACt                                                            <  1:1438107/95‑1 (MQ=255)                                                                                                        gCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTa                <  2:624293/139‑1 (MQ=255)                                                                                                          aGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGTTGCTTAAg              <  2:955332/139‑1 (MQ=255)                                                                                                          aGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAg              <  2:296751/139‑1 (MQ=255)                                                                                                             aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg           >  2:1121329/1‑139 (MQ=255)                                                                                                              tACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGaa                                                                                      >  2:877864/1‑63 (MQ=255)                                                                                                              tACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGaa                                                                                      <  1:877864/63‑1 (MQ=255)                                                                                                              tACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAgt          <  2:1310028/139‑1 (MQ=255)                                                                                                                      aaaaGCCATAACTGCCCAGCGGTGGTTATCTGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATg                                                                           <  1:123563/66‑1 (MQ=255)                                                                                                                      aaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATg                                                                           >  2:123563/1‑66 (MQ=255)                                                                                                                      aaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTAc  <  1:1226949/139‑1 (MQ=255)                                                                                                                        aaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAgt          >  1:1293642/1‑129 (MQ=255)                                                                                                                        aaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAgt          <  2:1293642/129‑1 (MQ=255)                                                                                                                           ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGtt                                                      >  1:1171011/1‑82 (MQ=255)                                                                                                                           ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGtt                                                      <  2:1171011/82‑1 (MQ=255)                                                                                                                                |                                                                                                                                 ATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTAC  >  NC_000913/1425853‑1426108

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