Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,220,971	G→A	pseudogene (1371/2648 nt)	ycgH →	putative transporter component, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,220,971	0	G	A	100.0%	185.5 / NA	52	pseudogene (1371/2648 nt)	ycgH	putative transporter component, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/20);  total (32/20)

AATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCT  >  NC_000913/1220834‑1221100                                                                                                                                          |                                                                                                                                  aaTTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGc                                                                                                                                 <  1:383075/140‑1 (MQ=255)        aCAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGAc                                                                                                                          >  2:716003/1‑140 (MQ=255)               tatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTa                                                                                                                   <  2:177256/140‑1 (MQ=255)                    aTCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCTGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                              <  1:309146/140‑1 (MQ=255)                    aTCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                              <  1:309147/140‑1 (MQ=255)                      cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAGTGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                            >  2:797525/1‑140 (MQ=255)                      cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                            >  2:532854/1‑140 (MQ=255)                      cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                            >  1:205777/1‑140 (MQ=255)                      cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                            >  1:367423/1‑140 (MQ=255)                      cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                            >  1:462824/1‑140 (MQ=255)                          cGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGt                                                                                                        >  1:182853/1‑140 (MQ=255)                                 tgatgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACtt                                                                                                 >  1:543212/1‑140 (MQ=255)                                   atgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCg                                                                                               >  1:3674/1‑140 (MQ=255)                                       tCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGtt                                                                                           <  2:853994/140‑1 (MQ=255)                                         aaaCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCa                                                                                         >  1:553463/1‑140 (MQ=255)                                          aaCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCaa                                                                                        >  1:100108/1‑140 (MQ=255)                                          aaCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCaa                                                                                        >  2:102287/1‑140 (MQ=255)                                          aaCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCaa                                                                                        >  2:452637/1‑140 (MQ=255)                                           aCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAAc                                                                                       <  1:874160/140‑1 (MQ=255)                                            cATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACt                                                                                      <  2:553463/140‑1 (MQ=255)                                                       ggAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTAtgt                                                                           >  1:769546/1‑140 (MQ=255)                                                          aaaaCATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGaa                                                                        <  1:578105/140‑1 (MQ=255)                                                                  aaCCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAAGTACCGTCTATGTGAAAAATGCTg                                                                >  2:639879/1‑140 (MQ=255)                                                                        aaaaGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGaa                                                          >  1:694190/1‑140 (MQ=255)                                                                               aaTTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgc                                                   >  2:581162/1‑140 (MQ=255)                                                                               aaTTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgc                                                   >  2:157868/1‑140 (MQ=255)                                                                                 ttATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgcgg                                                 >  1:637723/1‑140 (MQ=255)                                                                                      tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACg                                            <  2:769546/140‑1 (MQ=255)                                                                                       cGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGc                                           >  1:711872/1‑140 (MQ=255)                                                                                       cGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGc                                           >  1:860763/1‑140 (MQ=255)                                                                                        ggTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCt                                          <  2:543212/140‑1 (MQ=255)                                                                                           aataatGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTaaa                                       <  1:452637/140‑1 (MQ=255)                                                                                            ataatGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAAc                                      >  2:839830/1‑140 (MQ=255)                                                                                                tGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGa                                  <  1:839830/140‑1 (MQ=255)                                                                                                 ggTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGAt                                 <  1:102287/140‑1 (MQ=255)                                                                                                  gTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATc                                <  1:157868/140‑1 (MQ=255)                                                                                                   tCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCg                               >  2:581997/1‑140 (MQ=255)                                                                                                    cAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGa                              <  2:100108/140‑1 (MQ=255)                                                                                                    cAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGa                              <  1:233325/140‑1 (MQ=255)                                                                                                       cTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                           >  2:332978/1‑140 (MQ=255)                                                                                                       cTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACt                           <  2:637723/140‑1 (MQ=255)                                                                                                                cGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACtgt                  >  2:773833/1‑140 (MQ=255)                                                                                                                 gCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACtgtg                 >  1:880396/1‑140 (MQ=255)                                                                                                                 gCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACtgtg                 >  1:809665/1‑140 (MQ=255)                                                                                                                 gCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACtgtg                 >  1:537024/1‑140 (MQ=255)                                                                                                                  cACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGa                <  2:860763/140‑1 (MQ=255)                                                                                                                   aCCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGaa               <  2:694190/140‑1 (MQ=255)                                                                                                                        aTTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCg          >  1:883508/1‑140 (MQ=255)                                                                                                                          tAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGAt        <  1:383254/140‑1 (MQ=255)                                                                                                                               gatgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtc   >  1:205306/1‑140 (MQ=255)                                                                                                                                atgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtct  >  1:721245/1‑140 (MQ=255)                                                                                                                                atgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAAtct  >  1:721246/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                  AATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATCT  >  NC_000913/1220834‑1221100

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