Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,381,789	G→A	pseudogene (13/126 nt)	ycjV →	pseudogene, putative ATP‑binding component of a transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,381,789	0	G	A	100.0%	202.0 / NA	62	pseudogene (13/126 nt)	ycjV	pseudogene, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/29);  total (33/29)

AATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGC  >  NC_000913/1381666‑1381917                                                                                                                            |                                                                                                                                 aaTACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTg                                                                                                                      <  2:80759/136‑1 (MQ=255) aaTACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTg                                                                                                                      >  1:80759/1‑136 (MQ=255)  aTACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATg                                                                                                                                >  1:1602767/1‑125 (MQ=255)    aCGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGt                                                                                                                >  1:1331067/1‑139 (MQ=255)        ccGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTaaa                                                                                                            <  1:789689/139‑1 (MQ=255)         cGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAAt                                                                                                           <  1:631865/139‑1 (MQ=255)                     aTCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCa                                                                                               <  1:838705/139‑1 (MQ=255)                      tCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                              >  2:308262/1‑107 (MQ=255)                      tCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                              <  1:308262/107‑1 (MQ=255)                      tCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCg                                                                                                                           <  2:800830/110‑1 (MQ=255)                      tCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCg                                                                                                                           >  1:800830/1‑110 (MQ=255)                      tCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAg                                                                                              >  1:662765/1‑139 (MQ=255)                               gcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAata                                                                                     >  2:1593462/1‑139 (MQ=255)                                cAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAatat                                                                                    <  2:1331067/139‑1 (MQ=255)                                         aaTAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                 >  2:795555/1‑117 (MQ=255)                                         aaTAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                 <  1:795555/117‑1 (MQ=255)                                           tAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCAtt                                                                         >  2:1298476/1‑139 (MQ=255)                                                tttCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                              >  2:513500/1‑81 (MQ=255)                                                tttCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                              <  1:513500/81‑1 (MQ=255)                                                tttCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGat                                                                    <  2:606662/139‑1 (MQ=255)                                                  tCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTgg                                                                                                                     <  2:521369/88‑1 (MQ=255)                                                  tCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTgg                                                                                                                     >  1:521369/1‑88 (MQ=255)                                                   ccGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                            >  2:1103464/1‑112 (MQ=255)                                                   ccGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                            <  1:1103464/112‑1 (MQ=255)                                                   ccGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaaa                                                                                       <  1:427290/117‑1 (MQ=255)                                                   ccGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaaa                                                                                       >  2:427290/1‑117 (MQ=255)                                                             tAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTc                                                       <  1:380637/139‑1 (MQ=255)                                                                     cAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTa                                                                                                      >  2:336711/1‑84 (MQ=255)                                                                     cAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTa                                                                                                      <  1:336711/84‑1 (MQ=255)                                                                               cGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGaaa                                     >  2:1139419/1‑139 (MQ=255)                                                                                  tGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCAt                                                     >  2:1506557/1‑120 (MQ=255)                                                                                  tGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCAt                                                     <  1:1506557/120‑1 (MQ=255)                                                                                    cGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                              <  2:110621/45‑1 (MQ=255)                                                                                    cGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                              >  1:110621/1‑45 (MQ=255)                                                                                    cGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGa                                       <  1:1381947/132‑1 (MQ=255)                                                                                    cGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGa                                       >  2:1381947/1‑132 (MQ=255)                                                                                    cGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGaa                                >  2:1273061/1‑139 (MQ=255)                                                                                            aTGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt                        <  1:236419/139‑1 (MQ=255)                                                                                                 ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTg                                                                                                                      <  1:330204/40‑1 (MQ=255)                                                                                                 ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTg                                                                                                                      >  2:330204/1‑40 (MQ=255)                                                                                                 ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACg                                                              >  1:611187/1‑96 (MQ=255)                                                                                                 ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACg                                                              <  2:611187/96‑1 (MQ=255)                                                                                                  tACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCATAAACGGAAATATCAATTCGCTa                  >  2:1461988/1‑139 (MQ=255)                                                                                                  tACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTa                  <  1:1139419/139‑1 (MQ=255)                                                                                                   acacCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaaa                                                                                       >  2:1086608/1‑69 (MQ=255)                                                                                                   acacCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaaa                                                                                       <  1:1086608/69‑1 (MQ=255)                                                                                                    caccacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                 <  2:782758/58‑1 (MQ=255)                                                                                                    caccacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                 >  1:782758/1‑58 (MQ=255)                                                                                                    caccacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTaaa                >  2:406377/1‑139 (MQ=255)                                                                                                         cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt                        <  2:846798/126‑1 (MQ=255)                                                                                                         cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt                        >  1:846798/1‑126 (MQ=255)                                                                                                         cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCt                   <  2:76886/131‑1 (MQ=255)                                                                                                         cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCt                   >  1:76886/1‑131 (MQ=255)                                                                                                         cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAg           >  2:871474/1‑139 (MQ=255)                                                                                                         cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAg           >  2:1467215/1‑139 (MQ=255)                                                                                                          ggTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAgg          >  2:584077/1‑139 (MQ=255)                                                                                                            ttGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCAt                                                     <  2:144264/94‑1 (MQ=255)                                                                                                            ttGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCAt                                                     >  1:144264/1‑94 (MQ=255)                                                                                                            ttGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAgggg        >  1:307018/1‑139 (MQ=255)                                                                                                              gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGat                                                                    <  1:1483865/77‑1 (MQ=255)                                                                                                              gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGat                                                                    >  2:1483865/1‑77 (MQ=255)                                                                                                                  gCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGc  <  1:871474/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                 AATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGC  >  NC_000913/1381666‑1381917

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