Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,471,320	A→G	pseudogene (804/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,471,320	0	A	G	100.0%	208.4 / NA	58	pseudogene (804/3497 nt)	ydbA	pseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/32);  total (26/32)

GGCAGCAAATACTCTGGTAGGGATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGG  >  NC_000913/1471172‑1471465                                                                                                                                                     |                                                                                                                                                  ggCAGCAAATACTCTGGTAGGGATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCt                                                                                                                                                 >  1:99099/1‑151 (MQ=255)    agcaAATACTCTGGTAGGGATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTc                                                                                                                                              <  2:563656/151‑1 (MQ=255)    agcaAATACTCTGGTAGGGATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTc                                                                                                                                              >  2:708372/1‑151 (MQ=255)    agcaAATACTCTGGTAGGGATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTc                                                                                                                                              <  1:512368/151‑1 (MQ=255)     gcaAATACTCTGGTAGGGATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCg                                                                                                                                             >  2:345718/1‑151 (MQ=255)           aCTCTGGTAGGGATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCgtgt                                                                                                                                       >  2:738858/1‑151 (MQ=255)                    gggATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCgg                                                                                                                              <  1:356224/151‑1 (MQ=255)                      gATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCggg                                                                                                                             <  1:610192/150‑1 (MQ=255)                           cTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATgg                                                                                                                       >  2:473874/1‑150 (MQ=255)                           cTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATg                                                                                                                        >  2:75021/1‑150 (MQ=255)                            tGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATg                                                                                                                        >  1:269423/1‑149 (MQ=255)                             gTTTATGAGCATGGTACTGCTATCAACGACCAGACGTGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGATTTATAACGATGGCAAAGGAACAATTGTTAACTATGATACAATCTGCGCTTGCGGCGTGTGACAATCGGGGAAt                                                                                                                         >  1:703016/1‑147 (MQ=255)                             gTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGt                                                                                                                     <  2:479965/151‑1 (MQ=255)                              tttATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGt                                                                                                                     >  2:591623/1‑150 (MQ=255)                                   gAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACaa                                                                                                                 >  1:13254/1‑149 (MQ=255)                                     gCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACaataat                                                                                                             <  2:696909/151‑1 (MQ=255)                                      cATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAataata                                                                                                            <  1:177859/151‑1 (MQ=255)                                         ggTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACa                                                                                                          <  2:194366/150‑1 (MQ=255)                                         ggTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAg                                                                                                         <  1:46068/151‑1 (MQ=255)                                         ggTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAg                                                                                                         <  2:706654/151‑1 (MQ=255)                                          gTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAga                                                                                                        <  1:652578/151‑1 (MQ=255)                                           tACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATAc                                                                                                           >  2:367721/1‑147 (MQ=255)                                            aCTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAgat                                                                                                       <  1:318399/150‑1 (MQ=255)                                            aCTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAgat                                                                                                       <  2:147828/150‑1 (MQ=255)                                               gCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAgatgat                                                                                                    <  2:378349/150‑1 (MQ=255)                                                  aTCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTc                                                                                                 <  1:301741/150‑1 (MQ=255)                                                  aTCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCa                                                                                                <  2:624328/151‑1 (MQ=255)                                                    cAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCAc                                                                                               <  1:614293/150‑1 (MQ=255)                                                         aCCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACt                                                                                         <  2:131120/151‑1 (MQ=255)                                                          ccAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTg                                                                                        <  1:697398/151‑1 (MQ=255)                                                          ccAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTg                                                                                        <  2:106733/151‑1 (MQ=255)                                                                                   gTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGAtt                                                                >  1:514842/1‑150 (MQ=255)                                                                                   gTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTa                                                               <  2:477113/151‑1 (MQ=255)                                                                                             gTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCgg                                                     <  1:122733/151‑1 (MQ=255)                                                                                              tCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCgg                                                     >  2:132469/1‑150 (MQ=255)                                                                                               cAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGga                                                    <  2:79491/150‑1 (MQ=255)                                                                                                aGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCgg                                                     >  1:554949/1‑148 (MQ=255)                                                                                                  gCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGaa                                                 >  2:217702/1‑150 (MQ=255)                                                                                                    gTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACt                                              <  2:645822/151‑1 (MQ=255)                                                                                                    gTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACt                                              >  1:7660/1‑151 (MQ=255)                                                                                                       ttATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTg                                             <  2:699678/149‑1 (MQ=255)                                                                                                        tataACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTa                                           >  2:122465/1‑150 (MQ=255)                                                                                                        taaaaCGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTACGGCGTGTGCCAATCGGGGAACTAGTACAATAATACAGATGATGTCAATCCACTGATCTATACCGGTGGCGATACGATTAGACGCAGCAGCGGAACTGTa                                           >  1:70925/4‑150 (MQ=255)                                                                                                           aaCGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACg                                        <  1:707372/150‑1 (MQ=255)                                                                                                             cGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCt                                      >  2:410094/1‑150 (MQ=255)                                                                                                                 ggCACAGGAACAATTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCtaaata                                 <  2:379273/151‑1 (MQ=255)                                                                                                                             aTTGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATctgctgctg                     >  2:81871/1‑151 (MQ=255)                                                                                                                               tGTTAACTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCtgt                    <  2:353931/150‑1 (MQ=255)                                                                                                                                   aaCTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTg               <  1:286697/151‑1 (MQ=255)                                                                                                                                    aCTATGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTg               >  2:509028/1‑150 (MQ=255)                                                                                                                                       aTGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATa            >  2:29228/1‑150 (MQ=255)                                                                                                                                       aTGGTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATa            <  1:586880/150‑1 (MQ=255)                                                                                                                                         ggTACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAgc         <  2:336732/151‑1 (MQ=255)                                                                                                                                           tACAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAgctg       >  1:65212/1‑151 (MQ=255)                                                                                                                                             cAATCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAgctggc     >  2:18741/1‑151 (MQ=255)                                                                                                                                              aaTCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAgctggct    <  1:478857/151‑1 (MQ=255)                                                                                                                                                tCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGctggctgg  >  2:44020/1‑151 (MQ=255)                                                                                                                                                tCTGCGCTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGctggctgg  >  2:648896/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGCAGCAAATACTCTGGTAGGGATGGCTGTTTATGAGCATGGTACTGCTATCAACGACCAGACGGGTGTTATCAATATCAATGTTGGTACTGGTCAGGCGTTTTATAACGATGGCACAGGAACAATTGTTAACTATGGTACAATCTGCACTTTCGGCGTGTGCCAATCGGGGAATGAGTACAATAATACAGATGATTTCACCTCACTGATCTATACCGGTGGCGATACGATTACACGAAGCGGAGAAACTGTAACGCTAAATAAATCTGCTGCTGTGACTGATAAGCTGGCTGG  >  NC_000913/1471172‑1471465

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