Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,770,087 C→A pseudogene (55/114 nt) yibU → pseudogene, H repeat‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,770,0870CA98.7% 263.5 / NA 77pseudogene (55/114 nt)yibUpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (38/38);  total (38/39)
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

GAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAG  >  NC_000913/3769991‑3770181
                                                                                                |                                                                                              
gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                             <  2:1183573/100‑1 (MQ=255)
gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCGCATCCTGCGGAAGa                                                                                             >  2:1074771/1‑100 (MQ=255)
   aaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACa                                                                                          <  1:363382/100‑1 (MQ=255)
    aGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACaa                                                                                         >  1:654902/1‑100 (MQ=255)
     gCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAg                                                                                        <  1:1219136/100‑1 (MQ=255)
          ccTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGa                                                                                   <  1:112355/100‑1 (MQ=255)
          ccTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGa                                                                                   >  1:358712/1‑100 (MQ=255)
           cTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGAt                                                                                  >  1:459201/1‑100 (MQ=255)
             gCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaa                                                                                >  2:194714/1‑100 (MQ=255)
                 ttttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGCTGTATATTGTCATCCTGCGGAAGAACAAGACGCATaaa                                                                               >  2:1231858/1‑97 (MQ=255)
                  tttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAAt                                                                           <  1:1388663/100‑1 (MQ=255)
                   ttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATc                                                                          >  2:717360/1‑100 (MQ=255)
                   ttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATGGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATc                                                                          <  2:1010482/100‑1 (MQ=255)
                      tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAGAATCCaa                                                                       <  1:830837/100‑1 (MQ=255)
                      tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                       >  1:1049759/1‑100 (MQ=255)
                      tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                       <  2:415354/100‑1 (MQ=255)
                         ttAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATga                                                                    <  2:749939/100‑1 (MQ=255)
                          tAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                   >  2:1343554/1‑100 (MQ=255)
                          tAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                   <  1:761047/100‑1 (MQ=255)
                          tAGAGACGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCAACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                   <  2:887228/100‑1 (MQ=255)
                           agagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgaga                                                                  >  1:937253/1‑100 (MQ=255)
                            gagCCGTCGCAGGAGGAACAACTCATATCATCATTTCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGagag                                                                 <  1:856876/100‑1 (MQ=255)
                            gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGagag                                                                 >  2:815441/1‑100 (MQ=255)
                            gagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCGCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgag                                                                   >  2:395991/1‑98 (MQ=255)
                               cagTCGCAGGAGGAACAACACATATCATCATTGAAAACTCAACCATGTACAGTCTCATCATACGGAAGAACAAAACCGATAAAAAATCCGATgagat                                                                 <  2:1524119/95‑2 (MQ=255)
                                   cGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTc                                                          >  1:490817/1‑100 (MQ=255)
                                    gCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCa                                                         >  1:1354958/1‑100 (MQ=255)
                                    gCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCa                                                         >  2:631656/1‑100 (MQ=255)
                                    gCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCa                                                         <  1:1098036/100‑1 (MQ=255)
                                     cAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGCGCACACTCAt                                                        >  2:221532/1‑100 (MQ=255)
                                      aggaggAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATACAAAATCCAATGAGAGCACACTCATc                                                       >  1:94114/1‑100 (MQ=255)
                                      aggaggAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATc                                                       <  1:1150068/100‑1 (MQ=255)
                                          ggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGa                                                   >  1:331708/1‑100 (MQ=255)
                                            aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                 >  1:473936/1‑100 (MQ=255)
                                            aacaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAAc                                                 <  1:450811/100‑1 (MQ=255)
                                                 cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                            >  1:639242/1‑100 (MQ=255)
                                                 cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                            <  1:225541/100‑1 (MQ=255)
                                                  tCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGTGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatta                                           <  2:1010368/100‑1 (MQ=255)
                                                  tCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatta                                           <  2:209818/100‑1 (MQ=255)
                                                   cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                          <  2:1447798/100‑1 (MQ=255)
                                                   cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                          <  2:967780/100‑1 (MQ=255)
                                                   cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                          >  1:1387822/1‑100 (MQ=255)
                                                    atatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCTAGAACTTATTAAc                                         <  1:1074777/100‑1 (MQ=255)
                                                     tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                        <  1:1565154/100‑1 (MQ=255)
                                                     tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                        >  2:227143/1‑100 (MQ=255)
                                                      atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                       >  2:315769/1‑100 (MQ=255)
                                                      atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                       >  1:251631/1‑100 (MQ=255)
                                                         atcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAAt                                    <  2:297765/100‑1 (MQ=255)
                                                          tcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATg                                   <  2:937248/100‑1 (MQ=255)
                                                          tcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATg                                   <  1:1286905/100‑1 (MQ=255)
                                                            atTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGaa                                 >  1:1356026/1‑100 (MQ=255)
                                                                cGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATAtt                             >  2:185146/1‑100 (MQ=255)
                                                                 gATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTa                            >  1:799596/1‑100 (MQ=255)
                                                                 gATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTa                            >  2:1069641/1‑100 (MQ=255)
                                                                   tctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaaa                          >  2:379/1‑100 (MQ=255)
                                                                    ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                         >  1:564854/1‑100 (MQ=255)
                                                                      cgacgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaa                       >  2:1007300/1‑100 (MQ=255)
                                                                        acgaTGTATAGTCTCATACTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACACATCCAGAACTTATTAACATAATGAATATTAAAg                         <  2:1345469/96‑1 (MQ=255)
                                                                          gaTCTATAGTCTCATCCTGCCGAAGAACAATACGGATAAAAAAACCAATTAGAGCACACTCATCCACAACTTATTAACATAATGAATATTaaa                          <  2:460055/93‑1 (MQ=255)
                                                                           aTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATCTTAAAGAAAAatc                  <  1:1074184/100‑1 (MQ=255)
                                                                            tGTATAGTCTCATCCTGCGGAAGAACAAGACGGCTAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatca                 >  1:1294007/1‑100 (MQ=255)
                                                                            tGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatca                 <  1:156430/100‑1 (MQ=255)
                                                                            tGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatca                 <  1:460355/100‑1 (MQ=255)
                                                                              tataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatc               >  1:793500/1‑100 (MQ=255)
                                                                                  ggctcATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg           <  1:815446/98‑1 (MQ=255)
                                                                                  gTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg           >  1:840707/1‑100 (MQ=255)
                                                                                  gTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAg           <  2:530448/100‑1 (MQ=255)
                                                                                   tctcATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGc          <  1:1384146/100‑1 (MQ=255)
                                                                                   tctcATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGc          <  1:1384104/100‑1 (MQ=255)
                                                                                    ctcATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGc          <  1:1104332/99‑1 (MQ=255)
                                                                                     tcATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTa        <  2:720360/100‑1 (MQ=255)
                                                                                       aTCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATACTGAATATTAAAGAAAAATCATCGAAGCTAAt      >  1:1339460/1‑100 (MQ=255)
                                                                                       aTCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAAt      >  1:206339/1‑100 (MQ=255)
                                                                                         ccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGa    <  1:517099/100‑1 (MQ=255)
                                                                                         ccTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGa    <  1:624902/100‑1 (MQ=255)
                                                                                          cTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGaa   >  1:1041078/1‑100 (MQ=255)
                                                                                           tGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAg  >  1:743825/1‑100 (MQ=255)
                                                                                                |                                                                                              
GAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAG  >  NC_000913/3769991‑3770181

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: