Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,187,588	G→A	pseudogene (162/2439 nt)	yqiG →	pseudogene, fimbrial export usher family,putative membrane, Not classified, putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,187,588	0	G	A	100.0%	179.5 / NA	56	pseudogene (162/2439 nt)	yqiG	pseudogene, fimbrial export usher family,putative membrane, Not classified, putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/31);  total (25/31)

ATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGAT  >  NC_000913/3187465‑3187717                                                                                                                            |                                                                                                                                  aTTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGAc                                                                                                                    <  2:459437/139‑1 (MQ=255)       caaAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              <  2:162135/139‑1 (MQ=255)       caaAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              <  1:173906/139‑1 (MQ=255)       caaAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              <  2:492163/139‑1 (MQ=255)        aaaGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCTATCCACCTTGACCTCCAAt                                                                                                             <  1:1022858/139‑1 (MQ=255)         aaGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATg                                                                                                            <  1:107560/139‑1 (MQ=255)         aaGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATg                                                                                                            <  1:1029591/139‑1 (MQ=255)          aGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATgg                                                                                                           >  1:1088996/1‑139 (MQ=255)           gATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATggg                                                                                                          <  1:988014/139‑1 (MQ=255)              ttAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTAt                                                                                                       <  2:167829/139‑1 (MQ=255)                      aGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAAt                                                                                                             >  2:717836/1‑125 (MQ=255)                      aGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAAt                                                                                                             <  1:717836/125‑1 (MQ=255)                      aGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGAt                                                                                               >  1:119291/1‑139 (MQ=255)                      aGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGAt                                                                                               >  1:351234/1‑139 (MQ=255)                           aaGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCa                                                                                          <  2:1088996/139‑1 (MQ=255)                                  tGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCa                                                                                   >  2:443831/1‑139 (MQ=255)                                    aaaaCGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATTCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGt                                                                                 <  2:1085964/139‑1 (MQ=255)                                        cGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              >  2:1223448/1‑106 (MQ=255)                                        cGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              <  1:1223448/106‑1 (MQ=255)                                          ttAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                           <  2:351234/139‑1 (MQ=255)                                           tAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                           <  2:1306361/138‑1 (MQ=255)                                           tAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                           >  1:1306361/1‑138 (MQ=255)                                               cTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                      <  1:1357485/139‑1 (MQ=255)                                               cTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                      <  2:119291/139‑1 (MQ=255)                                                  tCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              >  2:963847/1‑96 (MQ=255)                                                  tCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              <  1:963847/96‑1 (MQ=255)                                                   cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAAtt                                                                                                  >  2:598377/1‑107 (MQ=255)                                                   cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAAtt                                                                                                  <  1:598377/107‑1 (MQ=255)                                                      aaTTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCa                                                               <  2:716606/139‑1 (MQ=255)                                                         ttGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAAtt                                                            <  1:443831/139‑1 (MQ=255)                                                                  aTGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattat                                                    >  2:1386320/1‑138 (MQ=255)                                                                  aTGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattat                                                    <  1:1386320/138‑1 (MQ=255)                                                                              ccGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACtga                                       <  1:926480/139‑1 (MQ=255)                                                                                  cGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGaa                                   >  2:1531017/1‑139 (MQ=255)                                                                                            tCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGaaa                         >  2:552500/1‑139 (MQ=255)                                                                                                  aGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCGCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTc                   >  2:728327/1‑139 (MQ=255)                                                                                                       tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCaa              >  1:353557/1‑139 (MQ=255)                                                                                                       tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCaa              >  2:522914/1‑139 (MQ=255)                                                                                                       tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCaa              >  2:1252520/1‑139 (MQ=255)                                                                                                         tGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              >  1:1009607/1‑41 (MQ=255)                                                                                                         tGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                              <  2:1009607/41‑1 (MQ=255)                                                                                                          gATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGCTGAATTTATTGAAAAAACTCAGCAAAAc           >  1:919533/1‑139 (MQ=255)                                                                                                            ttAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCt         >  1:377819/1‑139 (MQ=255)                                                                                                              aataataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCag                                                                   <  1:1178147/79‑1 (MQ=255)                                                                                                              aataataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCag                                                                   >  2:1178147/1‑79 (MQ=255)                                                                                                                 aataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAAt                                                             >  1:778388/1‑82 (MQ=255)                                                                                                                 aataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAAt                                                             <  2:778388/82‑1 (MQ=255)                                                                                                                   taaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCg                                                                            >  1:295822/1‑65 (MQ=255)                                                                                                                   taaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCg                                                                            <  2:295822/65‑1 (MQ=255)                                                                                                                   taaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGAt  <  2:353557/139‑1 (MQ=255)                                                                                                                   taaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGAt  >  2:604971/1‑139 (MQ=255)                                                                                                                   taaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAAGTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGAt  <  2:919533/139‑1 (MQ=255)                                                                                                                     agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGt                                                                                 >  1:1357818/1‑58 (MQ=255)                                                                                                                     agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGt                                                                                 <  2:1357818/58‑1 (MQ=255)                                                                                                                         gACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                           >  2:435576/1‑60 (MQ=255)                                                                                                                         gACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCa                                                                    <  2:1535909/67‑1 (MQ=255)                                                                                                                            |                                                                                                                                  ATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGAT  >  NC_000913/3187465‑3187717

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