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%	186.1 / NA	57	pseudogene (273/273 nt)	yafF	pseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/34);  total (23/34)

AGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCA  >  NC_000913/239243‑239506                                                                                                                                        |                                                                                                                                 aGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGtaataa                                                                                                                               >  2:1533617/1‑139 (MQ=25) aGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGtaataa                                                                                                                               >  1:1474222/1‑139 (MQ=18) aGGGATACGGCACATGTCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGtaataa                                                                                                                               <  1:1546222/139‑1 (MQ=255)      tACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATAccc                                                                                                                          <  2:1474222/139‑1 (MQ=255)           cacaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGAct                                                                                                                     >  2:1566777/1‑139 (MQ=255)                tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                                <  1:821636/139‑1 (MQ=255)                tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                                <  2:291599/139‑1 (MQ=255)                 gCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTcccc                                                                                                               >  1:291085/1‑139 (MQ=255)                  cTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                              <  1:1533617/139‑1 (MQ=255)                   tATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                              <  2:359835/138‑1 (MQ=255)                   tATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                              >  1:359835/1‑138 (MQ=255)                       aaTATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCtt                                                                                                         <  1:1252067/139‑1 (MQ=255)                         tatTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTaa                                                                                                       >  2:226464/1‑139 (MQ=255)                                aCGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAc                                                                                                <  2:485589/139‑1 (MQ=255)                                          aGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcac                                                                                      >  2:264245/1‑139 (MQ=255)                                           ggTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacc                                                                                     <  2:291085/139‑1 (MQ=255)                                            gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                    <  1:478378/139‑1 (MQ=255)                                            gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                    <  2:188284/139‑1 (MQ=255)                                            gTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCCTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                    <  2:1183605/139‑1 (MQ=255)                                             tATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcaccac                                                                                   <  1:354456/139‑1 (MQ=255)                                                      cAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAAc                                                                          >  2:1499261/1‑139 (MQ=255)                                                       aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcaccac                                                                                   <  1:377735/129‑1 (MQ=255)                                                       aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcaccac                                                                                   >  2:377735/1‑129 (MQ=255)                                                       aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACt                                                                         <  1:485390/139‑1 (MQ=255)                                                         ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCa                                                                       >  2:122396/1‑139 (MQ=255)                                                              aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGc                                                                  >  1:1218094/1‑139 (MQ=255)                                                               gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCa                                                                 <  1:264245/139‑1 (MQ=255)                                                                   tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCt                                                             >  1:40225/1‑139 (MQ=255)                                                                   tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCt                                                             >  1:63688/1‑139 (MQ=255)                                                                   tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCt                                                             >  2:910817/1‑139 (MQ=255)                                                                     aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTcccc                                                                                                               <  1:1116462/87‑1 (MQ=255)                                                                     aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTcccc                                                                                                               >  2:1116462/1‑87 (MQ=255)                                                                     aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                              >  2:294751/1‑88 (MQ=255)                                                                     aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                              <  1:294751/88‑1 (MQ=255)                                                                     aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGc                                                           >  2:233205/1‑139 (MQ=255)                                                                              aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGa                                                  <  2:1218094/139‑1 (MQ=255)                                                                                gcagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGaaa                                                <  1:910817/139‑1 (MQ=255)                                                                                gcagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGaaa                                                <  2:63688/139‑1 (MQ=255)                                                                                gcagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGaaa                                                <  2:40225/139‑1 (MQ=255)                                                                                     cATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACtttt                                           <  1:1566777/139‑1 (MQ=255)                                                                                            aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                    >  2:181092/1‑91 (MQ=255)                                                                                            aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                                    <  1:181092/91‑1 (MQ=255)                                                                                             gAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTatat                                   <  2:1266078/139‑1 (MQ=255)                                                                                              aaaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATAtt                                  >  1:479658/1‑139 (MQ=255)                                                                                                   aCCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGc                             >  2:679021/1‑139 (MQ=255)                                                                                                    ccTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTg                                                            >  1:502264/1‑107 (MQ=255)                                                                                                    ccTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTg                                                            <  2:502264/107‑1 (MQ=255)                                                                                                       ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg                         <  1:233205/139‑1 (MQ=255)                                                                                                          gTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtgctg                      <  1:226464/139‑1 (MQ=255)                                                                                                                  ttGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGAc              <  1:1499261/139‑1 (MQ=255)                                                                                                                  ttGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGAc              <  1:122396/139‑1 (MQ=255)                                                                                                                   tGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTctct                                       >  1:1326586/1‑113 (MQ=255)                                                                                                                   tGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTctct                                       <  2:1326586/113‑1 (MQ=255)                                                                                                                            cGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGcc                    <  1:518807/123‑1 (MQ=255)                                                                                                                            cGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGcc                    >  2:518807/1‑123 (MQ=255)                                                                                                                            cGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATc                 <  1:382049/126‑1 (MQ=255)                                                                                                                            cGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATc                 >  2:382049/1‑126 (MQ=255)                                                                                                                              ggCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCa  <  2:235252/139‑1 (MQ=255)                                                                                                                                    cGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcaccac                                                                                   <  1:1073709/52‑1 (MQ=255)                                                                                                                                    cGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcaccac                                                                                   >  2:1073709/1‑52 (MQ=255)                                                                                                                                        |                                                                                                                                 AGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCA  >  NC_000913/239243‑239506

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