Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,787,434	G→A	pseudogene (1038/2253 nt)	ygaQ →	uncharacterized protein, putative enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,787,434	0	G	A	100.0%	204.1 / NA	63	pseudogene (1038/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/31);  total (32/31)

ACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAT  >  NC_000913/2787321‑2787564                                                                                                                  |                                                                                                                                   aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                              >  2:1275626/1‑136 (MQ=255) aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                              <  1:1275626/136‑1 (MQ=255)   cTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcatt                                                                                                         >  1:364842/1‑139 (MQ=255)    tttAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTa                                                                                                        <  2:811372/139‑1 (MQ=255)      tAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTaaa                                                                                                      >  1:1471360/1‑139 (MQ=255)                    ttGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                        >  1:876270/1‑139 (MQ=255)                           tAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGt                                                                                 >  1:920528/1‑139 (MQ=255)                           tAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGt                                                                                 >  1:689210/1‑139 (MQ=255)                             gTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTaa                                                                               >  2:650943/1‑139 (MQ=255)                                   aaaaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                         <  2:1471360/139‑1 (MQ=255)                                        gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                    <  2:364842/139‑1 (MQ=255)                                         aCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   <  2:1159451/139‑1 (MQ=255)                                         aCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   >  1:1419236/1‑139 (MQ=255)                                         aCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   <  2:689210/139‑1 (MQ=255)                                         aCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   <  2:607195/139‑1 (MQ=255)                                          ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                  >  1:724222/1‑139 (MQ=255)                                            tGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                <  1:650943/139‑1 (MQ=255)                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                              >  2:246808/1‑81 (MQ=255)                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                              <  1:246808/81‑1 (MQ=255)                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   <  2:1403350/108‑1 (MQ=255)                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   >  1:1403350/1‑108 (MQ=255)                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                <  2:1505793/127‑1 (MQ=255)                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                >  1:1505793/1‑127 (MQ=255)                                                        ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                    <  2:1073277/139‑1 (MQ=255)                                                              atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATGCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                >  1:935911/1‑121 (MQ=255)                                                              atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                              >  2:1467506/1‑59 (MQ=255)                                                              atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                              <  1:1467506/59‑1 (MQ=255)                                                              atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                <  2:935911/121‑1 (MQ=255)                                                              atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                 <  1:1522135/120‑1 (MQ=255)                                                              atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                 >  2:1522135/1‑120 (MQ=255)                                                                     aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGAc                                       >  1:1315083/1‑139 (MQ=255)                                                                      gCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACt                                      >  2:336664/1‑139 (MQ=255)                                                                        gAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGt                                                                                              >  1:1433136/1‑81 (MQ=255)                                                                        gAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGt                                                                                              <  2:1433136/81‑1 (MQ=255)                                                                         aaGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAg                                                                                                     >  2:434022/1‑73 (MQ=255)                                                                         aaGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAg                                                                                                     <  1:434022/73‑1 (MQ=255)                                                                            aTCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                <  1:806708/139‑1 (MQ=255)                                                                                  aaaGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTa                          <  2:1419236/139‑1 (MQ=255)                                                                                      atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                         <  2:544555/88‑1 (MQ=255)                                                                                      atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                         >  1:544555/1‑88 (MQ=255)                                                                                          aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  <  2:1315083/139‑1 (MQ=255)                                                                                          aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  <  1:336664/139‑1 (MQ=255)                                                                                           gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  <  2:134132/138‑1 (MQ=255)                                                                                           gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  >  1:134132/1‑138 (MQ=255)                                                                                               atggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAtctc             >  2:957751/2‑139 (MQ=255)                                                                                                tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGa                                  >  1:725708/1‑117 (MQ=255)                                                                                                tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  >  2:1516173/1‑133 (MQ=255)                                                                                                tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  <  1:1516173/133‑1 (MQ=255)                                                                                                tagtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGa                                  <  2:725708/115‑1 (MQ=255)                                                                                                 ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                        >  1:774325/1‑62 (MQ=255)                                                                                                 ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                        <  2:774325/62‑1 (MQ=255)                                                                                                  gtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTaa          >  1:1120960/1‑139 (MQ=255)                                                                                                    aTCCCT‑GCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATaa                   <  1:316880/127‑1 (MQ=255)                                                                                                    aTCCCT‑GCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCAAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATaa                   >  2:316880/1‑127 (MQ=255)                                                                                                     tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGt                                                                                              <  1:828707/52‑1 (MQ=255)                                                                                                     tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGt                                                                                              >  2:828707/1‑52 (MQ=255)                                                                                                     tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                  <  2:747743/80‑1 (MQ=255)                                                                                                     tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                  >  1:747743/1‑80 (MQ=255)                                                                                                         tGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGAt                     >  2:1431674/1‑121 (MQ=255)                                                                                                         tGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGAt                     <  1:1431674/121‑1 (MQ=255)                                                                                                          ggCTGTTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCGATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAt  >  2:892947/1‑139 (MQ=255)                                                                                                            cTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                <  2:1054684/107‑1 (MQ=255)                                                                                                            cTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                >  1:1054684/1‑107 (MQ=255)                                                                                                                  |                                                                                                                                   ACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAT  >  NC_000913/2787321‑2787564

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