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_000913239,3780GA98.4% 198.8 / ‑6.6 63pseudogene (273/273 nt)yafFpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/32);  minor base C (0/1);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGA  >  NC_000913/239256‑239493
                                                                                                                          |                                                                                                                   
aTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACtctc                                                                                                     >  1:884769/1‑139 (MQ=255)
 ttGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAataata                                                                                                                <  2:303908/127‑1 (MQ=25)
 ttGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAataata                                                                                                                >  1:303908/1‑127 (MQ=25)
       ttAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATAcccc                                                                                                            >  2:679955/1‑125 (MQ=255)
       ttAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATAcccc                                                                                                            <  1:679955/125‑1 (MQ=255)
         aaTATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCtt                                                                                            >  2:221582/1‑139 (MQ=255)
              tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAaca                                                                                       >  1:183891/1‑139 (MQ=255)
              tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGAAGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAaca                                                                                       >  2:915031/1‑139 (MQ=255)
               ttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAacac                                                                                      <  2:474070/139‑1 (MQ=255)
                tGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                 >  1:97214/1‑127 (MQ=255)
                tGACGAATGAGAAGGTAGTCAAGGCAGGGTTAAGACGTAAGATGCGTAAAGCAGCCATGGGCAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                 <  2:97214/127‑1 (MQ=255)
                    gAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAccc                                                                                 <  1:420593/139‑1 (MQ=255)
                    gAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAccc                                                                                 <  1:915031/139‑1 (MQ=255)
                       tGAGAAGGTATTTAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCATTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   <  2:1095950/118‑1 (MQ=255)
                       tGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   >  1:1095950/1‑118 (MQ=255)
                         agaAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCGCCCCGTCCTTAAACACAACCCCCACt                                                                            >  2:823535/1‑139 (MQ=255)
                              gTATTCAAGTCAGGGTTAAGACGTAAGATGCGAACAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                       <  1:830514/139‑1 (MQ=255)
                               tATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAATCTACCTGGCGTCAGTCCTTGCGGCGAGCGGGCTTTCGTACTAATACCCCGCCTCTCCCCGTCCTTAAACACAACCCCCACTcaccac                                                                      <  1:823535/139‑1 (MQ=255)
                                aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACcaca                                                                     >  2:1027512/1‑139 (MQ=255)
                                  tCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAAc                                                                   <  1:1027512/139‑1 (MQ=255)
                                     aGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAc                                                                                   >  1:209886/1‑120 (MQ=255)
                                     aGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAc                                                                                   <  2:209886/120‑1 (MQ=255)
                                      ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   >  2:503590/1‑103 (MQ=255)
                                      ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   <  1:503590/103‑1 (MQ=255)
                                         aGGGTTAAGACGTAAGATGCGAAAAGCAGTCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACt                                                            >  1:1216183/1‑139 (MQ=255)
                                         aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCCGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                 >  1:514797/1‑102 (MQ=255)
                                         aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                 <  2:514797/102‑1 (MQ=255)
                                         aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAac                                                                                        >  2:628251/1‑111 (MQ=255)
                                         aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAac                                                                                        <  1:628251/111‑1 (MQ=255)
                                         aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACt                                                            <  2:1265905/139‑1 (MQ=255)
                                           ggTTACGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGAGGGGAGCGGGATTTCGTAATAATACCCCGATTCTCCTCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCa                                                          >  1:162569/1‑139 (MQ=21)
                                           ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTcccc                                                                                                  >  1:1038387/1‑99 (MQ=255)
                                           ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTcccc                                                                                                  <  2:1038387/99‑1 (MQ=255)
                                           ggTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCa                                                          >  1:957237/1‑139 (MQ=255)
                                            gTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCAt                                                         <  2:183891/139‑1 (MQ=255)
                                                aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTcc                                                                                                    >  2:386691/1‑92 (MQ=255)
                                                aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTcc                                                                                                    <  1:386691/92‑1 (MQ=255)
                                                aGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCGACTCACCACAACCTAAACTCATCCGc                                                     >  2:822740/1‑139 (MQ=255)
                                                      aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   >  2:203055/1‑87 (MQ=255)
                                                      aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   <  1:203055/87‑1 (MQ=255)
                                                      aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTg                                               >  2:1288868/1‑139 (MQ=255)
                                                       aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGtaataa                                                                                                                  >  1:424839/1‑71 (MQ=12)
                                                       aGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGtaataa                                                                                                                  <  2:424839/71‑1 (MQ=9)
                                                           gCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                       <  1:226394/110‑1 (MQ=255)
                                                           gCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacca                                                                       >  2:226394/1‑110 (MQ=255)
                                                            cGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATc                                                        >  2:1026887/1‑124 (MQ=255)
                                                            cGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATc                                                        <  1:1026887/124‑1 (MQ=255)
                                                                     gcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTGCTTAAACACAACCCCGACTCACCACAGCCTAAACTCATCCGCATCCTGCCATGCCGGAAACtt                                <  1:822740/139‑1 (MQ=255)
                                                                      ccATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACttt                               <  1:221582/139‑1 (MQ=255)
                                                                       cATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                 <  1:1241866/72‑1 (MQ=255)
                                                                       cATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                 >  2:1241866/1‑72 (MQ=255)
                                                                        aTGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTtc                             <  1:1288868/139‑1 (MQ=255)
                                                                                aaaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   <  1:503173/61‑1 (MQ=255)
                                                                                aaaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   >  2:503173/1‑61 (MQ=255)
                                                                                 aaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTc                    <  2:957237/139‑1 (MQ=255)
                                                                                   cTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTccc                  >  2:733807/1‑139 (MQ=255)
                                                                                        tGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCg                                       <  1:15296/113‑1 (MQ=255)
                                                                                        tGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCg                                       >  2:15296/1‑113 (MQ=255)
                                                                                         ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTACCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg            <  2:1216183/139‑1 (MQ=255)
                                                                                          gCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCa                                                                              <  2:19288/72‑1 (MQ=255)
                                                                                          gCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCa                                                                              >  1:19288/1‑72 (MQ=255)
                                                                                                   cTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGa  <  2:884769/139‑1 (MQ=255)
                                                                                                   cTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGa  <  1:733807/139‑1 (MQ=255)
                                                                                                      gCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   <  2:491477/39‑1 (MQ=255)
                                                                                                      gCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTccc                                                                                                   >  1:491477/1‑39 (MQ=255)
                                                                                                             gCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAccc                                                                                 >  2:1008038/1‑50 (MQ=255)
                                                                                                             gCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAccc                                                                                 <  1:1008038/50‑1 (MQ=255)
                                                                                                                          |                                                                                                                   
ATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGA  >  NC_000913/239256‑239493

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: