Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	239,378	G→A	pseudogene (273/273 nt)	yafF →	pseudogene, H repeat‑associated protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	239,378	0	G	A	100.0%	204.2 / NA	63	pseudogene (273/273 nt)	yafF	pseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/31);  total (32/31)

GGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCG  >  NC_000913/239244‑239513                                                                                                                                       |                                                                                                                                        gggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTaataat                                                                                                                                     >  1:52854/1‑139 (MQ=25)    aTACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATAcc                                                                                                                                  <  2:339895/139‑1 (MQ=37)        ggCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGa                                                                                                                              <  2:927891/139‑1 (MQ=255)            caTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGActct                                                                                                                          >  1:719325/1‑139 (MQ=255)               tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGa                                                                                                                              <  2:509827/132‑1 (MQ=255)               tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGa                                                                                                                              >  1:509827/1‑132 (MQ=255)                gCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTcccc                                                                                                                      <  2:719325/139‑1 (MQ=255)                    ttAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTcc                                                                                                                  >  1:149805/1‑139 (MQ=255)                           tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAaca                                                                                                           >  1:989428/1‑139 (MQ=255)                            ttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAacac                                                                                                          >  1:741643/1‑139 (MQ=255)                              gACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACaa                                                                                                        >  2:607679/1‑139 (MQ=255)                                   aTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAccccc                                                                                                   >  2:76459/1‑139 (MQ=255)                                            tATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcaccac                                                                                          <  1:76459/139‑1 (MQ=255)                                             aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGAct                                                                                                                            >  1:1120901/1‑104 (MQ=255)                                             aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGAct                                                                                                                            <  2:1120901/104‑1 (MQ=255)                                             aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTa                                                                                                               >  2:198224/1‑117 (MQ=255)                                             aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACGCTCCCCGTCCTTa                                                                                                               <  1:198224/117‑1 (MQ=255)                                              ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacc                                                                                            <  1:312779/135‑1 (MQ=255)                                              ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacc                                                                                            >  2:312779/1‑135 (MQ=255)                                              ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACaa                                                                                        >  1:920807/1‑139 (MQ=255)                                              ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACaa                                                                                        >  1:1147909/1‑139 (MQ=255)                                              ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACaa                                                                                        >  1:672857/1‑139 (MQ=255)                                                  aGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                                       <  2:1313873/104‑1 (MQ=255)                                                  aGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                                       >  1:1313873/1‑104 (MQ=255)                                                     cAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGTGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAAc                                                                                 >  2:711021/1‑139 (MQ=255)                                                      aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCt                                                                                     <  2:655056/134‑1 (MQ=255)                                                      aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCt                                                                                     >  1:655056/1‑134 (MQ=255)                                                          ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATc                                                                            <  1:723256/139‑1 (MQ=255)                                                              gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                           <  2:381410/120‑1 (MQ=255)                                                              gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                           <  2:725620/120‑1 (MQ=255)                                                              gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                           >  1:381410/1‑120 (MQ=255)                                                              gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                           >  1:725620/1‑120 (MQ=255)                                                                  tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAacac                                                                                                          <  1:1290922/101‑1 (MQ=255)                                                                  tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAacac                                                                                                          >  2:1290922/1‑101 (MQ=255)                                                                  tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCt                                                                    >  1:375238/1‑139 (MQ=255)                                                                  tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCt                                                                    >  1:234850/1‑139 (MQ=255)                                                                   aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                                       <  2:106789/87‑1 (MQ=255)                                                                   aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                                       >  1:106789/1‑87 (MQ=255)                                                                             aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGa                                                         <  2:52854/139‑1 (MQ=255)                                                                             aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGa                                                         <  2:989428/139‑1 (MQ=255)                                                                                      tGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTtct                                                <  2:1173083/139‑1 (MQ=255)                                                                                             aaaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCATAATAATACCCCGACTCTCGCCGTCCTTAAACACAACCCCCACTGACCACAACCTAAACTCATCCGCATCCTGCCTTGCCGGAAACTTTTCTCTATAtt                                         >  2:866195/1‑139 (MQ=39)                                                                                                      ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg                                <  2:672857/139‑1 (MQ=255)                                                                                                      ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg                                <  2:920807/139‑1 (MQ=255)                                                                                                      ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg                                <  2:385398/139‑1 (MQ=255)                                                                                                       gCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg                                <  1:564655/138‑1 (MQ=255)                                                                                                       gCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg                                >  2:564655/1‑138 (MQ=255)                                                                                                        cGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGa                                                                                                                              <  2:1097454/43‑1 (MQ=255)                                                                                                        cGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGa                                                                                                                              >  1:1097454/1‑43 (MQ=255)                                                                                                        cGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtgct                              >  1:598413/1‑139 (MQ=255)                                                                                                            aGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACt                                                                                                >  1:800146/1‑69 (MQ=255)                                                                                                            aGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACt                                                                                                <  2:800146/69‑1 (MQ=255)                                                                                                            aGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAAc                                                      <  1:123492/111‑1 (MQ=255)                                                                                                            aGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAAc                                                      >  2:123492/1‑111 (MQ=255)                                                                                                                cTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGa                      <  2:375238/139‑1 (MQ=255)                                                                                                                  tGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTtc                                                 <  2:306733/110‑1 (MQ=255)                                                                                                                  tGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTtc                                                 >  1:306733/1‑110 (MQ=255)                                                                                                                   gCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAg                   <  2:234850/139‑1 (MQ=255)                                                                                                                        gAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATc                                                                      <  2:752993/83‑1 (MQ=255)                                                                                                                        gAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATc                                                                      >  1:752993/1‑83 (MQ=255)                                                                                                                        gAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCg              >  1:797020/1‑139 (MQ=255)                                                                                                                        gAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCg              <  2:149805/139‑1 (MQ=255)                                                                                                                                   cGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTctct                                              >  1:665936/1‑96 (MQ=255)                                                                                                                                   cGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTctct                                              <  2:665936/96‑1 (MQ=255)                                                                                                                                    gtaataATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGcgcg  <  2:797020/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        GGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGCG  >  NC_000913/239244‑239513

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