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%	167.0 / NA	52	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 (24/28);  total (24/28)

TTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATGG  >  NC_000913/3187461‑3187719                                                                                                                                |                                                                                                                                    ttGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCt                                                                                                                          >  1:1301350/1‑139 (MQ=255) ttGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCt                                                                                                                          >  1:89482/1‑139 (MQ=255)  tGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATTCACCtt                                                                                                                         <  2:1210382/139‑1 (MQ=255)    aaTTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGa                                                                                                                       <  1:1206842/139‑1 (MQ=255)     aTTTAACAAAGATTTAATCGAAGCCCAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAACAATAAGAGGACACCAATCCACCTTGAc                                                                                                                      <  2:1041848/139‑1 (MQ=255)          acaaAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCa                                                                                                                 <  2:663237/139‑1 (MQ=255)           caaAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                <  2:1088335/139‑1 (MQ=255)                   tAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTAtt                                                                                                        <  1:1450865/139‑1 (MQ=255)                   tAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTAtt                                                                                                        <  1:220220/139‑1 (MQ=255)                                      tGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCa                                                                                     >  1:130925/1‑139 (MQ=255)                                             gTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCg                                                                              <  2:1374042/139‑1 (MQ=255)                                              ttAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                             <  2:130925/139‑1 (MQ=255)                                              ttAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCATCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                             <  2:1138078/139‑1 (MQ=255)                                                      tCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCag                                                                     >  2:21752/1‑139 (MQ=255)                                                       cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCag                                                                     <  2:968684/138‑1 (MQ=255)                                                       cacAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACTTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCag                                                                     >  1:968684/3‑138 (MQ=255)                                                               gAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTTATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTaa                                                            <  2:1301350/139‑1 (MQ=255)                                                                      aTGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattatt                                                     <  2:495422/139‑1 (MQ=255)                                                                      aTGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattatt                                                     <  1:21752/139‑1 (MQ=255)                                                                           cAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGAt                                                <  1:702109/139‑1 (MQ=255)                                                                              ttACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTa                                             >  2:810825/1‑139 (MQ=255)                                                                                     tCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGCCCTCCa                                                                                                                 >  2:729379/1‑64 (MQ=255)                                                                                     tCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCa                                                                                                                 <  1:729379/64‑1 (MQ=255)                                                                                      cGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                             >  1:707572/1‑99 (MQ=255)                                                                                      cGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                             <  2:707572/99‑1 (MQ=255)                                                                                        gCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACatta                                                       >  2:770260/1‑119 (MQ=255)                                                                                        gCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACatta                                                       <  1:770260/119‑1 (MQ=255)                                                                                                  aCTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGaaaaa                         >  2:1466217/1‑139 (MQ=255)                                                                                                      aGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTc                     >  1:496406/1‑139 (MQ=255)                                                                                                       gCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                >  1:114089/1‑47 (MQ=255)                                                                                                       gCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                <  2:114089/47‑1 (MQ=255)                                                                                                         aCTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTaa                                                            >  1:1205715/1‑97 (MQ=255)                                                                                                         aCTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTaa                                                            <  2:1205715/97‑1 (MQ=255)                                                                                                           tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAAcc                                                                                             >  1:711585/1‑62 (MQ=255)                                                                                                           tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAAcc                                                                                             <  2:711585/62‑1 (MQ=255)                                                                                                              gATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                        >  1:1586917/1‑80 (MQ=255)                                                                                                              gATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                        <  2:1586917/80‑1 (MQ=255)                                                                                                                   ataataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCa                                                                                     <  2:1584020/62‑1 (MQ=255)                                                                                                                   ataataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCa                                                                                     >  1:1584020/1‑62 (MQ=255)                                                                                                                   ataataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGAt        >  2:258886/1‑139 (MQ=255)                                                                                                                     aataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCg      >  1:342736/1‑139 (MQ=255)                                                                                                                      ataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGa     <  1:1466217/139‑1 (MQ=255)                                                                                                                       taaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGAt    <  1:258886/139‑1 (MQ=255)                                                                                                                         agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                             >  2:93288/1‑64 (MQ=255)                                                                                                                         agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                             <  1:93288/64‑1 (MQ=255)                                                                                                                         agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAAc                                                           >  1:1265613/1‑82 (MQ=255)                                                                                                                         agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAAc                                                           <  2:1265613/82‑1 (MQ=255)                                                                                                                         agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATgg  >  2:459860/1‑139 (MQ=255)                                                                                                                         agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATgg  >  1:759671/1‑139 (MQ=255)                                                                                                                         agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACATGATCGATgg  >  1:370326/1‑139 (MQ=255)                                                                                                                             gACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCa                                                                      <  2:660050/67‑1 (MQ=255)                                                                                                                             gACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCa                                                                      >  1:660050/1‑67 (MQ=255)                                                                                                                                |                                                                                                                                    TTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATGG  >  NC_000913/3187461‑3187719

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