Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,995,314	C→T	pseudogene (447/447 nt)	ygeN ←	pseudogene, orgB family, part of T3SS PAI ETT2 remnant

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,995,314	0	C	T	100.0%	105.5 / NA	33	pseudogene (447/447 nt)	ygeN	pseudogene, orgB family, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/16);  total (17/16)

TGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGT‑AAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTAT  >  NC_000913/2995183‑2995447                                                                                                                                     |                                                                                                                                      tgaAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAg                                                                                                                                 >  2:781029/1‑139 (MQ=255) tgaAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAg                                                                                                                                 >  2:23189/1‑139 (MQ=255)          aatGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATAATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAggg                                                                                                                        >  1:697128/1‑139 (MQ=255)              aaaTATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTg                                                                                                                    >  2:384459/1‑139 (MQ=255)              aaaTATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGACAGGGAGTTATAAGTCCTCAGGGCTTg                                                                                                                    >  1:628880/1‑139 (MQ=255)                 tataATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCa                                                                                                                 <  2:434091/139‑1 (MQ=255)                          aaataaaatTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAat                                                                                                        >  2:151864/1‑139 (MQ=255)                                                             aTTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAttttt                                                                     >  2:574152/1‑139 (MQ=255)                                                              ttAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                    <  2:176959/139‑1 (MQ=255)                                                                       cATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAaat                                                           <  1:200626/139‑1 (MQ=255)                                                                       cATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAaat                                                           <  2:628880/139‑1 (MQ=255)                                                                                atTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTgg                                                  >  1:265160/1‑139 (MQ=255)                                                                                          tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCt                                                                                                                             <  1:789814/54‑1 (MQ=255)                                                                                          tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCt                                                                                                                             >  2:789814/1‑54 (MQ=255)                                                                                          tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCata                                                                                                >  2:772747/1‑83 (MQ=255)                                                                                          tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCata                                                                                                <  1:772747/83‑1 (MQ=255)                                                                                          tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCa                                        >  2:774205/1‑139 (MQ=255)                                                                                          tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCa                                        >  2:626199/1‑139 (MQ=255)                                                                                            cAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGa                                      <  1:781029/139‑1 (MQ=255)                                                                                            cAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGa                                      <  1:151864/139‑1 (MQ=255)                                                                                               ttAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGc                                                                                                                   <  1:458945/59‑1 (MQ=255)                                                                                               ttAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGc                                                                                                                   >  2:458945/1‑59 (MQ=255)                                                                                                       tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAcc                           >  2:549044/1‑139 (MQ=255)                                                                                                            aaaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTTGTTAACCTTGTc                      >  2:197240/1‑139 (MQ=255)                                                                                                            aaaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTc                      <  1:384459/139‑1 (MQ=255)                                                                                                            aaaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTc                      <  2:265160/139‑1 (MQ=255)                                                                                                                cATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAggg                  <  2:234966/139‑1 (MQ=255)                                                                                                                    aaaaTTGAGCCAGGGAGTTATAAGTTCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGTAGATTTCGTTAACCTTGTCAGGGgaga              <  2:517368/139‑1 (MQ=255)                                                                                                                          gAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAAtgt                                                                                    >  2:491465/1‑63 (MQ=255)                                                                                                                          gAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAAtgt                                                                                    <  1:491465/63‑1 (MQ=255)                                                                                                                           aGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATg       >  2:285147/1‑139 (MQ=255)                                                                                                                                gggAGTTATTAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGattat  <  1:549044/139‑1 (MQ=255)                                                                                                                                gggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCGAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGattat  <  1:574152/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      TGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGT‑AAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTAT  >  NC_000913/2995183‑2995447

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: ‑