Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,288,561	(C)7→6	pseudogene (354/2525 nt)	yejO ←	pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,288,555	0	C	.	100.0%	168.6 / NA	40	pseudogene (360/2525 nt)	yejO	pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (22/18);  total (22/18)

TTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGA  >  NC_000913/2288425‑2288682                                                                                                                                   |                                                                                                                                tttGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccccccatc                                                                                                                         >  2:340549/1‑136 (MQ=255)      tCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCttttt                                                                                                                   >  1:566997/1‑139 (MQ=255)               cTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTa                                                                                                          >  1:717323/1‑139 (MQ=255)                           cgTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                              >  1:1017554/1‑139 (MQ=255)                           cgTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                              >  1:1022679/1‑139 (MQ=255)                              ggTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                           <  2:566997/139‑1 (MQ=255)                                   ttttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                                                      <  2:1017554/139‑1 (MQ=255)                                   ttttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                                                      <  1:340549/139‑1 (MQ=255)                                            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                   >  2:1003049/1‑117 (MQ=255)                                            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                   <  1:1003049/117‑1 (MQ=255)                                               cGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCt                                                                          >  2:491123/1‑139 (MQ=255)                                                gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                         <  2:544144/139‑1 (MQ=255)                                                 ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgc                                                                        >  2:877175/1‑139 (MQ=255)                                                            aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCt                                                             >  2:627317/1‑139 (MQ=255)                                                               tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccc                                                                                                                               <  2:797350/71‑1 (MQ=255)                                                               tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccc                                                                                                                               >  1:797350/1‑71 (MQ=255)                                                                aCCGCAAATGCAAATCCCTCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                         <  1:871398/139‑1 (MQ=255)                                                                  cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       <  2:1052112/139‑1 (MQ=255)                                                                    caaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                     >  2:1215729/1‑139 (MQ=255)                                                                    caaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                     <  1:1092225/139‑1 (MQ=255)                                                                                   ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTa                                      >  1:720727/1‑139 (MQ=255)                                                                                    ctctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAg                                     <  1:491123/139‑1 (MQ=255)                                                                                       tttAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCt                                  >  2:1056344/1‑139 (MQ=255)                                                                                               ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       <  2:627430/110‑1 (MQ=255)                                                                                               ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       >  1:627430/1‑110 (MQ=255)                                                                                               ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGcc                          >  1:588640/1‑139 (MQ=255)                                                                                                cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcg                         <  2:8094/139‑1 (MQ=255)                                                                                                 agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcgc                        >  1:1086592/1‑139 (MQ=255)                                                                                                     gtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                    <  2:1086592/139‑1 (MQ=255)                                                                                                      tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                    <  1:420845/138‑1 (MQ=255)                                                                                                      tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                    >  2:420845/1‑138 (MQ=255)                                                                                                        ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       >  2:585349/1‑101 (MQ=255)                                                                                                        ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       <  1:585349/101‑1 (MQ=255)                                                                                                          cccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       <  2:922336/99‑1 (MQ=255)                                                                                                          cccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       >  1:922336/1‑99 (MQ=255)                                                                                                           ccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCAtt              >  2:217142/1‑139 (MQ=255)                                                                                                                  gATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACaa       <  1:217142/139‑1 (MQ=255)                                                                                                                       cGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                         <  2:385664/84‑1 (MQ=255)                                                                                                                       cGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                         >  1:385664/1‑84 (MQ=255)                                                                                                                       cGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGa  >  2:665390/1‑139 (MQ=255)                                                                                                                       cGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGa  >  2:871825/1‑139 (MQ=255)                                                                                                                         aCAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       >  2:290447/1‑84 (MQ=255)                                                                                                                         aCAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                       <  1:290447/84‑1 (MQ=255)                                                                                                                                   |                                                                                                                                TTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGA  >  NC_000913/2288425‑2288682

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