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%	101.0 / NA	32	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/15);  total (17/15)

AATGAAATATAATTAAAAATAAAATTTTTGCGT‑AAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTAT  >  NC_000913/2995192‑2995447                                                                                                                            |                                                                                                                                      aatGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAggg                                                                                                                        <  1:789143/139‑1 (MQ=255)                                 taaaaaaaaTACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAATTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAattatt                                                                                                    >  2:600212/3‑127 (MQ=255)                                 taaaaaaaaTACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAATTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAattatt                                                                                                    <  1:600212/125‑1 (MQ=255)                                                    aTTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGa                                                                                        <  2:381850/120‑1 (MQ=255)                                                    aTTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGa                                                                                        >  1:381850/1‑120 (MQ=255)                                                       aaaaaaTCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCa                                                                  >  2:946849/1‑139 (MQ=255)                                                            aTCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTaa                                                             >  2:1463003/1‑139 (MQ=255)                                                                       atTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGTCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTgg                                                  >  1:1144037/1‑139 (MQ=255)                                                                       atTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTgg                                                  >  2:179825/1‑139 (MQ=255)                                                                               atTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCg                                          <  1:1463003/139‑1 (MQ=255)                                                                                 tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCa                                        >  1:1287208/1‑139 (MQ=255)                                                                                           ttATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTa                              >  1:827235/1‑139 (MQ=255)                                                                                             atatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                            <  2:827235/139‑1 (MQ=255)                                                                                             atatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                            <  1:946849/139‑1 (MQ=255)                                                                                                   aaaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                    >  2:30616/1‑93 (MQ=255)                                                                                                   aaaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                    <  1:30616/93‑1 (MQ=255)                                                                                                     aTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCa                                                                              <  1:456506/81‑1 (MQ=255)                                                                                                     aTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCa                                                                              >  2:456506/1‑81 (MQ=255)                                                                                                     aTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAaat                                                           >  1:293578/1‑100 (MQ=255)                                                                                                     aTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAaat                                                           <  2:293578/100‑1 (MQ=255)                                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                    >  1:1428844/1‑88 (MQ=255)                                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                    <  2:1428844/88‑1 (MQ=255)                                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttt                                                                      <  1:499599/86‑1 (MQ=255)                                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttt                                                                      >  2:499599/1‑86 (MQ=255)                                                                                                        aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAgggg                 <  1:179825/139‑1 (MQ=255)                                                                                                          tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGTTTGTTCGCAGATTTCGTTAAcc                           >  2:1199508/1‑127 (MQ=255)                                                                                                          tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGTTTGTTCGCAGATTTCGTTAAcc                           <  1:1199508/127‑1 (MQ=255)                                                                                                                 gAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAAt        >  2:151096/1‑139 (MQ=255)                                                                                                                  aGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATg       >  2:133158/1‑139 (MQ=255)                                                                                                                       gggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGattat  <  2:292008/139‑1 (MQ=255)                                                                                                                        ggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGt                       <  2:1221469/117‑1 (MQ=255)                                                                                                                        ggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGt                       >  1:1221469/1‑117 (MQ=255)                                                                                                                            |                                                                                                                                      AATGAAATATAATTAAAAATAAAATTTTTGCGT‑AAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTAT  >  NC_000913/2995192‑2995447

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