Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,130,088	T→G	pseudogene (56/1101 nt)	yghO ←	pseudogene, IS‑interrupted

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,130,088	0	T	G	100.0%	198.1 / NA	54	pseudogene (56/1101 nt)	yghO	pseudogene, IS‑interrupted
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (23/31);  total (23/31)

CAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCG  >  NC_000913/3129940‑3130235                                                                                                                                                     |                                                                                                                                                    cAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAgcg                                                                                                                                                   >  2:669373/1‑151 (MQ=255)     gTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCgg                                                                                                                                               <  2:245237/151‑1 (MQ=255)      tATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCggg                                                                                                                                              >  2:428535/1‑151 (MQ=255)          tATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATa                                                                                                                                           >  2:764113/1‑150 (MQ=255)           aTTTGCGCGGTTATACGCGCAACTGTCTGCCCCGCTTTTTTTGCCGCCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCCCAGACAAATGCTCGTTGCATTCAATGAATAAAGGGGGTATCCAGTTTTGAGCGGCGGGATaa                                                                                                                                          <  2:368295/150‑1 (MQ=255)               gcgcGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGtg                                                                                                                                      >  2:381149/1‑150 (MQ=255)               gcgcGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGtg                                                                                                                                      <  2:523771/150‑1 (MQ=255)                 gcgGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCGAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGtgatg                                                                                                                                   <  2:510887/151‑1 (MQ=255)                  cgGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGtgatg                                                                                                                                   >  2:435890/1‑150 (MQ=255)                      aaTACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGa                                                                                                                              <  1:884718/151‑1 (MQ=255)                         aCGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGaaaa                                                                                                                           <  2:597301/151‑1 (MQ=255)                          cGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGaaaa                                                                                                                           <  2:749422/150‑1 (MQ=255)                          cGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAg                                                                                                                          <  2:141947/151‑1 (MQ=255)                           gCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGc                                                                                                                         <  1:455740/151‑1 (MQ=255)                              cAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGAt                                                                                                                      <  2:590837/151‑1 (MQ=255)                              cAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGAt                                                                                                                      <  2:100309/151‑1 (MQ=255)                               aaCTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGAt                                                                                                                      >  2:286751/1‑150 (MQ=255)                               aaCTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGAt                                                                                                                      >  1:252007/1‑150 (MQ=255)                                aCTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATa                                                                                                                     >  1:119146/1‑150 (MQ=255)                                  tATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAAt                                                                                                                  >  2:407558/1‑151 (MQ=255)                                   aTCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAAt                                                                                                                  <  1:363683/150‑1 (MQ=255)                                       gCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCt                                                                                                              <  1:210073/150‑1 (MQ=255)                                              tttttttGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTc                                                                                                       >  1:456060/1‑150 (MQ=255)                                                tttttGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCAt                                                                                                     <  1:869739/150‑1 (MQ=255)                                                   ttGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCAtttt                                                                                                  <  2:118889/150‑1 (MQ=255)                                                    tGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCAttttt                                                                                                 <  1:174170/150‑1 (MQ=255)                                                                  gCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTg                                                                                    <  2:533957/149‑1 (MQ=255)                                                                   ccAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGa                                                                                 <  1:372667/151‑1 (MQ=255)                                                                    cAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGa                                                                                 >  2:908747/1‑150 (MQ=255)                                                                       atgatATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACaa                                                                             <  2:458502/151‑1 (MQ=255)                                                                          atatGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTc                                                                          >  1:562378/1‑151 (MQ=255)                                                                            atGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTccc                                                                        >  2:509789/1‑151 (MQ=255)                                                                                     ccccGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                 >  2:208262/1‑149 (MQ=255)                                                                                      cccGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATgaga                                                              <  2:860423/151‑1 (MQ=255)                                                                                       ccGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                 >  1:875364/1‑147 (MQ=255)                                                                                        cGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                 >  2:443738/1‑146 (MQ=255)                                                                                             ttttCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGt                                                        <  1:220563/150‑1 (MQ=255)                                                                                               ttCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGttt                                                      <  2:855252/150‑1 (MQ=255)                                                                                               ttCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGttt                                                      <  1:286777/150‑1 (MQ=255)                                                                                                 cGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTc                                                   >  1:310283/1‑151 (MQ=255)                                                                                                   cAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCAt                                                 <  2:446237/151‑1 (MQ=255)                                                                                                         aaTGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTgg                                             >  2:123571/1‑149 (MQ=255)                                                                                                               cGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAg                                     <  2:750138/151‑1 (MQ=255)                                                                                                                  tGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGaa                                   <  1:697588/150‑1 (MQ=255)                                                                                                                  tGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAAc                                  >  1:487511/1‑151 (MQ=255)                                                                                                                    cGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAAc                                  >  1:766590/1‑149 (MQ=255)                                                                                                                        cAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGt                             <  1:19802/150‑1 (MQ=255)                                                                                                                                aaaGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATg                     >  1:695902/1‑150 (MQ=255)                                                                                                                                           ccAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACtt         <  2:197692/151‑1 (MQ=255)                                                                                                                                             aGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTAc       <  1:446630/151‑1 (MQ=255)                                                                                                                                             aGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTAc       <  2:620269/151‑1 (MQ=255)                                                                                                                                                ttGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGGACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCtt    >  1:796858/1‑151 (MQ=255)                                                                                                                                                ttGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCtt    >  1:665430/1‑151 (MQ=255)                                                                                                                                                  ggAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTTCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCg  <  2:356764/151‑1 (MQ=255)                                                                                                                                                  ggAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCg  >  1:606422/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAAGGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCG  >  NC_000913/3129940‑3130235

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