Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,473,811	A→G	pseudogene (3295/3497 nt)	ydbA →	pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,473,811	0	A	G	100.0%	149.7 / NA	46	pseudogene (3295/3497 nt)	ydbA	pseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/25);  total (21/25)

GGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAAC  >  NC_000913/1473666‑1473953                                                                                                                                                  |                                                                                                                                               ggAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGAAAAGACGGCTGGAGAGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGACGGCg                                                                                                                                             >  1:710254/1‑149 (MQ=255)          aaCCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGGCTACGCAGGTAAAGACGGCGGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGGGCTACGCGAAGAGCCGGGGTACGGCAAggc                                                                                                                                  >  2:355618/1‑149 (MQ=255)           aCCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCt                                                                                                                                 <  1:343912/151‑1 (MQ=255)            ccTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTg                                                                                                                                <  1:851569/151‑1 (MQ=255)               gAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTgg                                                                                                                               <  1:288292/149‑1 (MQ=255)                aaTATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGgc                                                                                                                              >  2:607702/1‑149 (MQ=255)                  tatGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGgcaggc                                                                                                                          >  1:598310/1‑151 (MQ=255)                   atGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGgcaggcg                                                                                                                         <  2:807386/151‑1 (MQ=255)                   atGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGgcaggcg                                                                                                                         >  2:117966/1‑151 (MQ=255)                       aCAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGgcaggcgcagg                                                                                                                     >  2:425776/1‑151 (MQ=255)                           tGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCgcaggcag                                                                                                                  >  1:97019/1‑150 (MQ=255)                            ggCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTc                                                                                                                >  2:549174/1‑151 (MQ=255)                             gCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCa                                                                                                               <  2:730173/151‑1 (MQ=255)                              cAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAg                                                                                                              <  1:122991/151‑1 (MQ=255)                                 cGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCa                                                                                                            >  1:720111/1‑150 (MQ=255)                                       ggcggTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAAc                                                                                                     >  1:505837/1‑151 (MQ=255)                                             ggACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGAt                                                                                               >  2:783677/1‑151 (MQ=255)                                             ggACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGAt                                                                                               <  1:410705/151‑1 (MQ=255)                                                 aGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGAc                                                                                            <  1:433668/150‑1 (MQ=255)                                                      cgtcgGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTc                                                                                        >  2:275890/1‑149 (MQ=255)                                                      cgtcgGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCa                                                                                       <  2:505873/150‑1 (MQ=255)                                                      cgtcgGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCa                                                                                       <  2:488224/150‑1 (MQ=255)                                                      cgtcgGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAg                                                                                      <  2:156921/151‑1 (MQ=255)                                                         cgGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAg                                                                                   <  1:815040/151‑1 (MQ=255)                                                            gCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGggc                                                                                >  1:739153/1‑151 (MQ=255)                                                               gAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGggcgg                                                                              <  1:165776/150‑1 (MQ=255)                                                                   cTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGggcggcggc                                                                          >  2:132113/1‑150 (MQ=255)                                                                      gACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATc                                                                       >  1:760751/1‑150 (MQ=255)                                                                      gACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCa                                                                      <  2:139309/151‑1 (MQ=255)                                                                                   aaGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACa                                                           >  1:415574/1‑149 (MQ=255)                                                                                       cGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTc                                                     >  1:755807/1‑151 (MQ=255)                                                                                            ggAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCg                                                 <  1:498963/150‑1 (MQ=255)                                                                                             gAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTg                                               <  2:695035/151‑1 (MQ=255)                                                                                                 gctagcgctaCGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGt                                           >  2:320762/1‑151 (MQ=255)                                                                                                      cgctaCGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTAcagc                                      >  2:560520/1‑151 (MQ=255)                                                                                                              tGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCaa                                 >  1:642214/1‑148 (MQ=255)                                                                                                                    gCGGGCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTc                        <  2:261075/151‑1 (MQ=255)                                                                                                                      gggCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGc                      <  1:415489/151‑1 (MQ=255)                                                                                                                       ggCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGc                      <  1:355595/150‑1 (MQ=255)                                                                                                                       ggCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCt                     <  2:190905/151‑1 (MQ=255)                                                                                                                        gCCGAACCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTg                    <  1:549137/151‑1 (MQ=255)                                                                                                                            aaCCTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATc                <  2:742546/151‑1 (MQ=255)                                                                                                                              ccTGAGCTACGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTg              >  1:103126/1‑151 (MQ=255)                                                                                                                                      aCGCGAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGt       <  2:609545/150‑1 (MQ=255)                                                                                                                                          gAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAAc  <  1:622828/151‑1 (MQ=255)                                                                                                                                          gAAGAGCCGGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAAc  <  2:802425/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               GGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAAC  >  NC_000913/1473666‑1473953

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